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PREFACE 


This report on the feasibility of utilizing SKYLAB spacecraft 
photography to provide control for small scale mapping' opera- 
tions was prepared so as to enhance the comprehension of those 
readers not familiar with the principles of analytic aerotri- 
angulation procedures and the SKYLAB Earth Resources Experiment 
mission. 

All of the work involved in this study was performed in the 
offices of, and on equipment operated by, the National Oceanic 
and Atmospheric Administration, National Ocean Survey/Coastal 
Mapping Division, which is located in the Washington Science 
Center, Rockville, Maryland. Computer processing was performed 
on a CDC 6600 computer operated by NOAA and located at Suitland 
Maryland. 

The author wishes to express his sincere appreciation to Mr. D. 
Norman and Mr. I. Raborn of the Aerotriangulation Section who 
performed the photocoordinate measurements, assembled the data, 
and processed t.he material through the analytic aerotriangula- 
tion system of computer programs. Thanks are due to Commander 
W. V. Hull, Chief of the Coastal Mapping Division, and to Mr. 
e. Slama, Chief of the Photogrammetric Research Branch, for 
allowing the author to provide the time and effort needed to 
bring this study to a successful conclusion. A final vote of 
thanks is due Mrs. M. Taglieri , secretary for the Photogram- 
metric Research Branch, for her patience and diligence in 
preparing this report for delivery to the National Aeronautics 
and Space Administration. 
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ABSTRACT 


The objective of this study was to investigate the feasibility 
of utilizing SKYLAB spacecraft Earth Terrain Camera (S-190B) 
1:946,000 scale photography in analytic aerotriangulation 
procedures to provide low-order, high-density- control suitable 
for small-scale mapping operations . 

The long range application is the employment of this technique 
for coastal zone mapping at medium and small scales, surveys 
in remote areas, forest and range management, yarious planning 
activities, and route location for highways, pipelines, trans- 
mission lines, and canals. 

The National Oceanic and Atmospheric Administration, National 
Ocean Survey, (NOAA/NOS)., office-identified the locations of 
29 photo control points of known position and elevation on a 
strip of 12 photographs ranging along a 350-mile track from 
Charlotte , North Carolina , • to the Rappahannock River in 
Virginia. The coordinates of pertinent images on each photo- 
graph were then processed through an established analytic 
aerotriangulation system of computer programs. 

The inherent errors in using nonmetric SKYLAB photography and 
office-identified photo control made it necessary to perform^ 
numerous block adjustment solutions involving different combin 
ations of control and weights. The final block adjustment was 
executed holding to 14 of the office-identified photo control 
points. The accuracy of the solution was evaluated by 
comparing the analytically computed ground positions of the 
15 withheld photo control points with their known ground 
positions and also by determining the standard errors of these 
point's from the variance values. A horizontal position RMS 
error of 15 meters was attained. The maximum observed error 
in position at a control point was 25 meters. 
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BACKGROUND 

A basic framework of horizontal and vertical geodetic control 
is essential for coordinating surveys and the mapping of large 
areas. In the United States, the first- and second-order 
horizontal and vertical control surveys conducted by the 
National Ocean Survey of the National Oceanic and Atmospheric 
Administration provide this basic framework of geodetic control. 
Additional control surveys of third-order accuracy by various 
federal and state agencies then subdivide or extend the- basic • 
network by triarigulation, traverse, and leveling methods in 
order to bring the control into the areas to be mapped. The 
control stations established by these surveys are usually 
monumented for future use. 

Photogrammetry is a system of measuring and interpreting data 
recorded on photographs and is applicable to all sciences that 
depend on reliable geometric measurements of physical quantities 
occurring in a fixed or transitory state. The widest application 
of' the photogramme trie art has been in the topographic mapping 
of the earth surface, where it provides an alternative and/or 
supplement to conventional ground methods for establishing 
geodetic control and mapping geographical features. 

Photogrammetric mapping procedures offer certain ’worthwhile 
advantages over ground methods, such as: 1. Detailed mapping 

can be performed more accurately, completely, efficiently, and 
economically than by ground methods. This advantage increases 
with the complexity of the details, as, for example, in city 
and harbor areas and along irregular coasts j 2. otherwise 
inaccessible areas can be more easily mapped; 3. maps can 
generally be produced from photographs with less ground control • 
than is necessary for ground methods; and 4. the workload is 
transferred to the office, where operations are independent 
of weather and daylight. 

Photogrammetry can provide the following primary services: 

Provide three-dimensional stereoscopic models of the 
terrain that can be set in stereoscopic plotting instru- 
ments, so that planimetric and topographic details can 
be compiled from these models. 

2 . Extend the basic control network directly into the area 
of photography by using aerotriangulation methods to., 
bridge between the high-order arcs of existing control. 
This procedure yields primarily fourth-order nonmonumented 
control and minimizes the need for field work to establish 
the photo control required to properly orient the. stereo- 
scopic mod'els oh the .plotting instruments. 
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Aerotriangulation is a photogrammetric techni.que for deriving 
the ground" co^dinates of objects from a set of overlapping 
aerial photographs that show images of these objects and also 
of a relatively sparse distribution of other objects whose 
coordinates are known from classical measurements on the ground. 
The two principal methods employed to determine the desired three- 
dimensional ground coordinates for the objects are stereotriangu- 
lation and analytic aerotriangulation. Stereotriangulation 
depends on measurements made on a sequential series of overlapping 
stereoscopic models formed on a high-precision photogrammetric 
plotting instrument. ■ Analytic aerotriangulation is a digital - ■ 
solution based on observed coordinates of the images created by 
pertinent objects appearing on each of the photographs covering 
the area. The analytic solution possesses a remarkably high 
accuracy potential as compared to stereotriangulation, because 
of the advantages accruing from automation, digital accuracy, 
least- squares adjustment, and freedom from the mechanical 
discrepancies contributed by the stereoscopic plotting instruments 
In addition, the systematic errors such as camera- lens distortion, 
film shrinkage, atmospheric refraction distortion, etc., can 
be more effectively eliminated by analytic* methods than in 
stereotriangulation procedures. A disadvantage of the analytical 
solution is that the computations are complicated and- require 
a large-size electronic computer to contain and process the 
large volume of data with economy and speed. 

THE MATHEMATICAL BASIS OF ANALYTIC AEROTRIANGULATION 

Several different variations in analytical aerotriangulation 
techniques have evolved. However, all of the methods basically 
consist of writing equations which relate the unknown elements ^ ' 
of exterior orientation of each photograph to camera constants 
and refined x and y image coordinates observed on a comparator. 

The equations are solved for the 'unknown camera orientation 
parameters and the ground coordinates for each object creating 
observed images on the photographs. Since more observational 
information is normally available than is re'quired for a unique 
solution, the method of least squares is used to obtain the most 
probable values of the unknown parameters in such a fashion 
that the sum of the squares of the residual observational dis- 
crepancies is a minimum. 

The observation equations must- be linear with respect to the 
unknown independent parameters; otherwise, a direct solution of 
the equations becomes difficult. If the mathematical model is 
nonlinear, as most photogrammetric problems are, a Taylor's 
expansion series is usually employed to linearize the equations. 
The' computation requires initial approximations of the unknown 
parameters and is iterative because the second and higher 
degree terms of- the Taylor's series are neglected to -simpl-ify- 
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the mathematics. The least squares solution provides corrections 
to the approximate values of the parameters. If the initial 
approximations are coarse, the corrections are added to them, 
■giving fresh and improved approximations for a 'new solution. 

Least squares is used again to provide another set of corrections, 
and the procedure is repeated until some criterion of convergence 
is satisfied. 

The most commonly used analytic methods are designed to enforce ■ 
one of two conditions: coplanarity or collinearity . ^ Coplanarity 

is the condition that the two perspective centers of an over- 
lapping pair of photographs, any object point, and its correspond- 
ing image points on the two pictures all lie ’in a common plane ; 
i.e., the rays passing through the two camera stations should 
intersect at a single object point. The purpose of the computa- 
tion is to minimize the distance between the two rays at the 
object location. The observation equation utilized in each 
object space angle thus contains the, four residual errors, v^, v^. 

on first photo and ,v on second photo, involved in measuring 
the X and y image coordinates on each picture- This causes the 
solution to become cumbersome and difficult to solve properly 
when the object occurs on several or more' photographs. 

Collinearity is the condition that every object, its photographic 
image, and the camera exposure station must lie on a common 
straight line, as defined by the method of least squares in which . 
the sum of the squares of the residual errors of image coordinate 
measurement is minimized. Two observation equations are written 
for every image and, except for the few control station images,' 
•contain only one residual error, . v^ or v^, in each equation. 

This simplifies the application of least squares so that any 
number of photographs can be routinely accomraodate-d^ 

The principle of collinearity provides the basis for the NOS 
-method of analytic aerotriangulation. The condition is utilized 
in an iterative manner to determine incremental corrections to 
initial approximations for the unknowns, which are reasonably 
close to the correct values. 

COLLINEARITY CONDITION FORMULATION 

The well known equations of collineation comprising the projective 
transformation are : 

_ (X-Xo) a,, + (y-7o) + ( 2 -Z ^ ) 

T (X-X,) ^31 +iY-Yo)a>, + (Z-2o) 

(X-Xq) g„ (y-n) 033. + 

Z (X-Xo) Qai + (.Y-Yo ) Qaa + C33 
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where the a- terms are the nine elements of the rotation matrix 
relating the x:, z image coordinate system to the X, Y, Z 
ground coordinate system of the objects-, and Xo , Yo,; Zo are the- 
coordinates of the camera station expressed in the ground 
coordinate system. 
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A Taylor’s expansion series. is applied to the transcendental 
collinearity equations to obtain linearized observation equations 
which can then be solved for the linear independent unknowns by 
the application of least squares. The complete - form of the 
linearized observation equations is: 

where the nine terms dw through ^ are incremental corrections 
to be applied to initxal approximations of the unknowns.. These 
two equations occur for each image on each photograph. Block 
adjustment requires the presence of all nine terms, whereas the 
space resection computations use only the first' six terms dw 
through dZo- . ' Sufficient photographs , and images are needed to 
provide at least as many equations as there are unknowns to be - 
computed. The solution is iterative and terminates when the 
incremental corrections to the angular parameters are smaller 
than the observed precision. 

THE NATIONAL OCEAN SURVEY ANALYTIC AEROTRIANGULATION SYSTEM 
Aerial Photography ■ 

The Coastal Mapping Division of NOS utilizes well calibrated 
precision aerial cameras to secure the near vertical aerial 
photography required for its photo gramme trie operations-. These 
cameras include the Wild RC-8, 6-inch focal length camera and 
the Wild RC-10 camera equipped with interchangeable 6-inch and 
3.5-inch focal length cones. 

The cameras are mounted in two aircraft operated by the Division; 
one being a' DeHavilland Buffalo, while the other is an Aero 
Commander 69 OA. The Buffalo is a twin-engined turboprop craft, 
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having a cruising speed between 120 and ^180 knots at altitudes 
UP to 32,000 feet. Its cruising range is 10 hours. ihe cabin 
is unpressurized, thereby requiring the crew to use oxygen above 
10 000 feet. The Buffalo has been modified to simultaneously 
accommodate "three aerial cameras -mounted in three hatches. The. 
plane is used on nearly all photographic missions conducued by 

the Division. 


The Aer o Commander 69 OA is a leased aircraft, employed to obtain 
photography for airport surveys. It is a twin-engined turbo- 
prop craft, having a cruising speed between 100 and 270 knots 
at altitudes up xo 32,000 feet. Its cruising range is five hours. 
The pressurized cabin permits the crew to operate without relying 
on oxygen. "However, this requires the photography he taken 
through an optically, flat glass window which covers, the single 
hatch. A backup oxygen supply is available ^ that permits aeria _ 
photography to be taken without the window in place over the hatch. 


Photo Control Points 

In order to implement analytic aerotriangulation methods, it is 
necessary to establish sufficient photo control to properly 
orient the aerial photography. Photo control refers to the 
establishment of horizontal positions and/or elevations with 
respect to the basic framework of geodetic control monuments of 
carefully chosen ground objects, which create sharp, distinct, 
and easily identifiable point images on the photographs. The 
positions and/or elevations are determined from the monumented 
control network by third- and fourth-order triangulation, _ traverse , 
and leveling methods. The photo control points selected in the 
field are usually prominent natural or cultural features providing 
sharp imagery on the overlapping photographs. Such examples are 
road intersections, fence-line intersections, lone trees, corners 
of buildings and wharves, and smaller stacks or towers. 

Por the more precise photo gramme trie surveys, the placement of 
'specially prepared targets on the control points, prior- to 
photography, is desirable for facilitating accurate office 
identification of the photo control points. The targets are 
symmetrical in design, centered on the photo control points, and 
are of sufficient size to show on the photography. The targets 
are usually in the shape of .a Y or a cross . 


Pass Points 

' ’ X 

In analytical control extension, pass points are established in 
the nine standard relative orientation ^ locations on each photo- 
graph, *The pass points may be easily identifiable point images 
similar to those selected as photo control points. Usually, 
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however j . the pass points are siniply holes drilled into the ^ 
photographic emulsion, with the Wild PUG-2 stereoscopic point 
transfer device in areas providing an optimum stereoscopic- 
perception. On 60 percent overlap photography, one ground 
object can appear as a pass point image on three successive 
photographs of a strip. 

NOS programs permit two pass points to be used in each relatxve^ 
orientation location, even -though only one pass point is .suffi- 
cient to provide all of the data needed for .the analytic computations 
Then, if one of the pass points in a relative orientation location 
should exhibit an excessively large residual discrepancy during the 
solution, it can be discarded and its companion pass point 
substituted in its place. ' • 

Marking and 'Photocoordinate Measurement 

A stereoscopic point transfer device, such as the Wild PUG-2, is 
first used to select, mark, and transfer suitable photo control 
and pass point images to adjacent photographs. The coordinates 
of each image on the photograph transparency are then measured 
to micron accuracy by either a digitized Mann monocomparator 
or by a stereo comparator, such as the Wild STK or the Zeiss PSK. 

The use of a mo no c omp ar a t or means that the’ images on only one 
plate can be measured at a time; It, thus, is necessary to _use 
a stereoscopic point transfer device to mark and transfer all 
images to every photograph. Stereocomparators , however, can 
simultaneously measure the photocoordinates of corresponding , 
images bn a stereoscopic pair of photographs. As a consequence, 
it is necessary to mark (drill) all of the images on only one 
photograph in each strip in order-’to specifically identify the 
image for measurement purposes.- The' drilled .image can then be 
stereoscopically' transferred to the second photograph and measured 
with the stereocomparator without actually ever drilling the 
image on the second photo. • Since the images are drilled on only 
one plate , it is not necessary to use a stereoscopic point 
transfer device to mark the images . In practice , the PUG device 
is needed only as a means for transferring the images to any 
adjoining strips. 

In essence, a stereocomparator consists of two mechanically united 
monocomparators. A monocomparator is, therefore, more accurate 
-than a stereocomparator. In practice, however, the need ^ for 
additional PUG operations with a monocomparator results- in the 
accuracy of the monocomparator-PUG combination, being about equal 
to that of the stereocomparator. 

After the photocoordinates have been measured, the data is 
submitted for processing through the series of computer programs 
comprising the analytic aerotriangulation system. 
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Computer Processing 

The analytic aerotriangulation system developed at. NOS consists 
of five programs: Cl) Image coordinate refinement and three- 

photo orientation, (2) strip adjustment to ground control., 

C3) secant plane coordinate transformation, (M-) block- adjustment, 
and C5) accuracy analysis-. 

(1) Image Coordinate Refinement and Three-PhotO' •Orientataon : 

To obtain the highest possible degree of accuracy in analytical 
solutions, measured photocoordinates must be corrected for system- 
atic errors which cause distortions in the image positions. The 
first program, therefore, begins with a refinement of the raw x 
and y image coordinates measured on each photograph on the compar-- 
ator. The popularity of making positive prints of the photo- 
graphs on glass plates (diapositi'ves) has declined in recent 
years. Today polyester plastic bases are used because they 
provide a dimensionally stable base film that is much less 
sensitive to humidity, temperature, and laboratory processing. 

The photocoordinates observed on- the photograph transparency 
are corrected for the systematic distortions introduced by the 
comparator, film shrinkage, camera lens, and atmospheric refrac- 
tion. The problem of earth curvature is recognized in the third 
program, in which the ground coordinates of all objects are 
expressed in a geocentric, three-dimensional secant ■ plane ■ system 
that takes earth curvature into account. " ■ 

The refined image coordinates are -punched out to. serve as input 
to the block adjustment program (4). The refined image coordinates 
theoretically should be nearly all free of systematic error 
and contain only residual observational discrepancies in them. 

The program then proceeds to the three-photo camera orientation 
phase, which comprises an interrelated geometric fitting of the 
photographs based only on the refined -image coordinates ‘and is 
entirely independent from any ground control data. The computa- - 
tion is iterative and derives the orientation of each photograph 
relative to the previous two in the strip. It also determines 
the positions of all pertinent objects in a three-dimensional 
coordinate system at the scale of the photography. The collinea- 
tion principle is imposed in a least squares solution that 
minimizes the discrepancies in the observed .image coordinates. 

The residual errors are analyzed. by the computer, which discards, 
those images exhibiting excessively large discrepancies. The 
removal of these blunders provides "clean" image coordinate 
data for all subsequent compu'tations , 

(2) Strip Adjustment to Ground Control : The analysis of three 

photographs at a time automates the joining of the separate 
triplets into a continuous strip and develops a set of model 
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coordinates, which are analogous to the product obtained from 
conventional stereotriangulation on stereoscopic plotting 
instruments. The horizontal and vertical strip adjustment program 
transforms the model coordinate data into the prevailing ground 
control coordinate system by fitting to control stations through 
the application of polynomial equations and least squares. Any 
large residual discrepancies appearing in the resulting adjustment 
are corrected in .order to obtain provisional ground position data 
that are free of blunders prior to entering the block adjustment 
computation. 

The analytic computations may be terminated after strip adjustment 
or may continue through block adjustment, depending on the desired 
accuracy. While block adjustment can be performed without 
actually using the three-photo orientation and strip adjustment 
programs, these preliminary programs are employed in practice 
to furnish improved and complete data for the block adjustment 
in an effort to reduce the time and cost of "debugging” and 
computer operations. 

(3) Secant Plane Coordinate Transformation : If maximum accuracy 

is desired, the provisional ground coordinates are first trans- 
formed into a geocentric and then into a special secant plane 
system that takes earth curvature into account. The block adjust- 
ment solution is performed using these secant plane coordinates 
for the objects, together with the previously obtained refined 
image coordinates. The secant plane transformation program is 
designed to operate in its inverse mode so that given secant 
plane coordinates can be transformed back into the prevailing 
ground coordinate system. 

(4) Block Adjustment : This program permits the simultaneous 

solution of the absolute orientation (three linear elements of 
position and three angular elements of orientation) of all 
photographs in a block of overlapping strips of photography. 

Only the pass points and control station objects contribute 
equations and thus influence the least square orientation solution. 
Their finalized ground coordinates are computed simultaneously, 
along with the absolute orientation of all the photographs in the 
block. Those objects that are not pass points or control stations 
do not contribute equations and thus do not influence the orien- 
tation solution. The finalized ■ camera parameters from the 
orientation solution and the refined image coordinates are used 

to compute the final ground X, Y, Z coordinates for these other 
objects by intersection. After the block adjustment is completed, 
the adjusted secant plane coordinates are transformed back into 
the original ground coordinate system by applying the secant 
plane transformation in its inverse mode. 
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The major task of the block adjustment program is the solution 
of the large number of simultaneous equations in a least square 
manner that efficiently utilizes the memory capacity of the 
computer. The largest program written at NOS can accommodate as 
many as 600 photographs in a single simultaneous least square 
adjustment. Some 36,000 observation equations containing about 
10,000 unknowns may be generated in developing the normal 
equations. The number of objects whose final ground positions 
can be computed in the block adjustment solution is unlimited. 

All of the analytic programs have been written to operate on the 
CDC 6600 computer. To date, the largest block adjustment problem 
processed through the computer contained 180 photograhs. The 
solution involved over 15,000 observation equations and about 
4,500 unknowns. The CDC computer running time for the least 
square block adjustment was less than five minutes. 

(5) Accuracy Analysis : In order to appreciate fully the 

accuracy potential of the system, and the error values at test 
points, a final computer program is used to develop the inverse 
of the matrix of normal equations, the variances, and the 
standard errors in centimeters in X, Y, and Z at all the points 
used throughout the area. The error E = Qe at any point is 
composed of two components where ^ is the variance at the point 
as derived from the inverse, and ^ is the standard error of unit 
weight for the problem based on program (4). The quantity Q is 
affected by the geometry of the system, including the amount and 
distribution of control points, and e is related to the precision 
of the steps of the system including image resolution. 

Accuracy of the NOS Analytic Aerotriangulation System 

A significant increase in accuracy results when analytic aerotri- 
angulation computations are continued through block adjustment. 
Studies conducted at NOS have yielded the following accuracy 
results : 

The horizontal position root-mean-square error in. meters when 
using film cameras is SlO-5 where ^ is the denominator of the 
photography scale' fraction. If a glass plate camera is 
used, and the pass points are premarked, then the expected RMS 
error in meters is about 1/5 of SlO-5 . 

The vertical and horizontal errors are equal for a 3. 5 -inch 
focal length camera (base-height ratio = 1) , whereas the 
vertical errors may be about 1.5 times larger than the 
horizontal errors for a 6- inch focal length camera (base- 
height ratio =0.6). 
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This accuracy is achieved when: (a) 60 percent forward and 

side overlap exists between the photographs in the block; 

(b) a strong network of horizontal and vertical photo control 
exists around the perimeter of the area, along with a few 
interior vertical photo control stations; Cc) all of these 
stations are premarked prior to photography; (d) the block 
consists of at least three strips of photography. 

Note: The rms errors are considered to be .essentially equal 

to the standard error of unit weight. The rms value therefore 
has a 68 percent reliability. The 90 percent reliability is 
about 1.6 times larger, and the 93 percent value is about 
‘2.6 times lareer. 
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SKYLAB 

The SKYLAB mission was planned and implemenired to determine the 
capability of man and spacecraft to conduct medical, solar physics 
stellar and solar astronomy, and Earth observational programs. 

On May 14, 1973, the National Aeronautics and Space Administra- 
tion successfully launched the SKYLAB manned orbital facility 
into a nearly circular 234— nautical mile (435 km) orbit above 
the earth (SL-l). The first three-man team of astronauts 
manned the laboratory for 23 days, beginning on May 25, 1973, 
CSL-2); the second team occupied the facility for 60 days, 
starting on July 28, 1973 (SL-3); and the third team followed on 
November ^16, 1973, for an 85-day mission (SL-4). The 50-degree 
inclination of the orbit permitted the astronauts to view 75 
percent of the earth's surface — the area between 50 degrees North 
and 50 degrees South — and to pass over a given point once every 
five days. 

The 100-ton SKYLAB spacecraft is actually a hollowed-out third 
stage of a Saturn rocket, originally assigned to the U.S. moon 
program, which has been converted to provide living and working 
space for the astronauts. Within the SKYLAB space station are 
complex scientific and technical instruments that will enable them 
to conduct investigations directed toward the accomplishment of 
medical experiments, solar astronomy experiments, technical 
experiments, and earth resources experiments as follows: 

a. To study man: Medical experiments will determine 

physiology conditioning and performance capability in • 
real time, in zero-gravity environment, for long- 
duration space flight. 

b. To study the sun: Solar astronomy experiments will 

provide a synoptic survey and study of special phenomena 
on the solar disk in X-ray, ultraviolet (uv), and visible 
spectral wavelengths. 

c. To study space technology: Technical experiments will 

evaluate coating degradation, spacecraft contamination, 
manufacturing and repair techniques , and manned- 
maneuvering units . 

d. To study the Earth: Earth resources experiments will 

provide a synoptic survey of selected areas on the earth 
in visible, infrared (IR), and microwave spectral wave- 
lengths . 
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EARTH RESOURCES EXPERIMENTS 

cirviAR -investisations, the Earth Resources 
Among the various are concerned directly and 

exclusively with vLious^plants , ground scenes, 

energy spectral distributions, not only 

and infrared and microwave portions of the 

in the visible but in the t^rr ’’signatures" can be 

electromagnetic spectrum, These P ^ “ multiband camera 

detected by utilizing Six such 

to infrared spectrometers and micro systems for observing 

electronic and photographic remot ?<h<=r,-nr^ces Experiment 

the earth have been combined into the |agy|^ 

Package (EREP) and mounted on ^^^ard the m stabilize 


cfsKlLAf Is r sriS.po"grlne;; 

^SoecSit.^S musf^ maLuvIrid iSto an earth-oriented mode xn 
order to use the EREP sensors. 

The EREP is designed J I^wide^rSe^of applica- 

space for use by a vane y mp-ni- * The EREP sensors can be 

tions to earth resources ‘ ^ dL on the 

operated singly or in various co such as weather and/or 

scientific requirements or other f ^^%ape and film 
vehicle capability limits. ^ 'hr’in^ back to earth the data 

so that each team of astronauts initial 

recorded during its stay on the spacecra * -^exas, the 

processing at the Johnson Space ^®5^®'"oK?L2rSiniipal investi- 
data are distributed to more than 200 SKYLAB pnncipa 

gators , ' ■ 

The data acquired fY^j^^^^olt^larth^Resouroes disoi- 

for studlee SbsSrvSuoS Sn be applied to 

plines. For example, these _ ecology, geology, geography, 
research in agriculture, ±o ^irio-'ST- bvdrography , oceanography, 

meteorology , the purpose of identifying ■ 

cartography, and similar ii _ rrr,oM+h rates’ assessing crop 
agricultural species i measuring g determining land surface 

vigor and stress; classifying land use^ Assessing water 

composition and structure, ^ shorelines, and estii- 

^rSns^TefiouS-- Toiditionl, and similar projects . 

The SKXLAB. earth resources ^as ^een structured into ten 

major disciplines as outlined in Table 1 beio 
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TABLE 1. 

EREP PROGRAM STRUCTURE 


100 AGRICULTURE/RANGE /FORESTRY 600 

Crop inventory 
Insect infestation 
Soil type' 

Soil moisture 
Range inventory 
Forest inventory 
Forest insect damage 

700 

200 GEOLOGICAL APPLICATIONS 
Mapping 

Metals exploration 
Hydrocarbon exploration 
Rock types 

Volcanoes 800 

Earth movements 

300 CONTINENTAL WATER RESOURCES 

Ground water 
Snow mapping ' 

Drainage basins 
Water quality 

.400.. OCEAN INVESTIGATIONS 900 


COASTAL ZONES, -SHOALS, AND BAYS 

Circulation and pollution in bays 
Underwater topography and 
sedimentation 
Bathymetry 
Coastal circulation 
Wetlands ecology 

REMOTE SENSING TECHNIQUES DEVELOPMENT] 

Pattern recognition 
Microwave signatures 
Data processing 
Sensor performance evaluation 

REGIONAL PLANNING AND DEVELOPMENT 

Land use classification techniques 
Environmental impacts - special 
topics 

State and foreign resources 
Urban applications 
Coastal/plains applications 
■ Mountain/desert applications 

CARTOGRAPHY 


Sea state 
Sea/Lake ice 
Currents 
Temperature 
Geodesy 

Living marine resources 

500 ATMOSPHERIC INVESTIGATIONS 

Storms j fronts , and clouds 
Radiant energy balance 
Air quality 
Atmospheric effects 


Photomapping 
Map revision 
Map accuracy 
Thematic mapping 

000 USER AGENCY TASKS 

Department of the Army 
Department of the Interior 


Each SKYLAB investigation has been given a three-digit task number, 
according to the subdisciplines in which work is done. In addition, 
every EREP study site has been given a three-digit designation • 
that defines its geographic location. 
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DESCRIPTION OF EARTH RESOURCES EXPERIMENT PACKAGE CEREP) SENSORS 

An extensive description of the photographic remote sensing system 
is given below, along with a brief description of the other 
EREP sensors. 

Multispectral Photographic Facility (S-190A and S-190B) 

The experiment objective is to photograph the earth ’ s ■ surf ace in 
a spectral range that includes visible light and extends into 
the near-infrared, with sufficient resolution and spectral defini- 
tion to allow detailed analysis and interpretation by specialists 
in a variety of earth resources' disciplines. 

The facility is arranged in two parts. S-190A consists of an 
array of six 70-inm film cameras, precisely matched and boresighted, 
so that photographs from all six cameras will be accurately in 
register. Thus, all of the features seen in one photograph can 
be simultaneously aligned with the same features in the photo- 
graphs from the other cameras. A combination of black-and-white 
and color films is used in conjunction with selective filters for 
spectral analysis, allowing comparison with imagery obtained with 
the IR spectrometer (S-191) and multispectral scanner (S-192) and 
with the Earth Resources Technology Satellite (ERTS). The 
camera array is mounted behind an optical glass window, just for- 
ward of the radial docking hatch in the Multiple Docking Adapter 
(MDA) . 

The second part, S-190B Earth Terrain Camera, consists of a single 
camera that is located behind an optical glass window in the 
-Scientific Airlock (SAL) on the antisolar side of the Orbital 
Workshop (OWS). This camera is an adaptation of the Lunar Topo- 
graphic Camera carried on the Apollo 14 mission. 

Controls for the six-camera array are integrated with the controls 
for the other EREP sensors located in the MDA. However, the Earth 
Terrain Camera controls are mounted on the side of the camera 
housing and are independent of other EREP sensors. 

For earth resources operations, SKYLAB departs from its normal 
solar orientation to an orbital mode that provides for continuous 
pointing of the cameras and other sensors at the ground directly 
below. The crewmen load film, install filters, set up the camera 
controls, remove the covers. from the camera ports, uncover the 
window, install the Earth 'Terrain Camera in the Scientific Airlock, 
and make other preparations for camera operations . 

The exposed film is the primary data returned at the end of each 
SKYLAB mission, for processing and analysis on the ground. 
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S-190A 

The EREP multispectral photographic camera consists of six high- 
precision 70-mm film cameras with matched distortion and focal 
length. The f/2.8 lenses have a focal length of six inches. 

The camera has a field-of-view of 21.2 degrees across the 
flats based on the photographic format size of 2.25 inches 
square and provides ground coverage of 16 3 k'm square per frame 
at a scale of 1:2,850,000. The system is designed for the 
following wavelength/ film combinations; 


0.5 

to 

0.6 

um 

PAN X BSW 

0.6 

to 

0.7 

um 

PAN X BSW 

0.7 

to 

0.8 

um 

IR B£W 

0.8 

to 

0.3 

um 

IR BSW 

0.5 

'to 

0.88 

um 

IR COLOR 

0.4 

to 

0.7 

um 

HI-RES COLOR 


S~39 0B 

The body of the Earth Terrain Camera (ETC) is an extensively 
modified Hycon KA-74 reconnaissance camera body with a bidirec- 
tional focal-plane shutter and vacuum film flattening.- The ETC 
is equipped with an f/4 lens having a focal length of 460 ram, 
(18 inches), color correction, and a maximum radial distortion 
of 10 urn. 'Forward image-motion compensation is provided by 
rocking the entire camera in its mount during the exposure. 

The ETC has a limited field-of-view of 14 degrees across the 
flats, based on the photographic format size of 4.5 inches 
square and provides ground coverage of 109 km square per 
frame at a scale of 1:946,000. The system is designed to 
utilize the following film types: 


ESTIMATED GROUND 
RESOLUTION 

TYPE DESCRIPTION WAVELENGTH, urn' Cat low ' contrast ) 


SO-242 
EK 3414 
EK 3443 


Aerial color, high 
resolution 
High-definition - 
aerial B£W 
Aerochrome IR, 
color 


0.4 to 0.7 
0,5 to 0.7 
0.5 to 0 ; 8 8 


70 ft. 
5 5 ft. 
100 ft. 


on ground 
on ground 
on ground 


The ETC is not a metric camera in the photogrammet-ric sense. 
Because the image frame is a part of the removable film magazine 
and because of the use of a focal-plane shutter,, the geometric 
quality of the photographs is limited. The shutter motion is in 
the flight direction for one exposure, and opposite the flight 
direction for the next exposure. This causes a slight scale 
compression or stretching in the flight direction, depending 
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principal point cannot be precisely 
therefore analytxcal applications are limite d, 
camera xs operated for 60 percent overlap » the base- 
hexght ratio IS only 0.10; thus, any stereoscopic height measure- 
ments from the photographs have especially limited accuracy. 

In spite ^ of these limitations, the ETC represents a significant 

systems for earth resources observations from 
photography having a ground resolution 
from 3 to 20 tunes better than that of any other space photographic 

of the ETC i ^ the primary objective 

™ high-resolution stereoscopic photography 

sensors by aiding in the interpretation 

of data gathered by them. 

/ S°f® discussion of the ETC is given in a paper 

_ s!i90B EASH"raSllNlAMi^-!|erippe“nd?^ a!"^“ 

Infrared Spectrometer (S-191 ) 

objective of this experiment is to mahe a fundamental 
evaluation of the applicability and usefulness of sensing earth 
resources from orbital altitudes in the visible through near- 
SKYfAR'^QnSoJ far- inf raped spectral regions. Correlation of 

gathered by ground-based and 

its ^^sure that the radiance from the target and 

accurately established. The extent to 
hich the effects of the atmosphere can be removed from the data 

aLurs^JJ r^n ^ importance to all remote sensors, and this 

will be quantitatively tested-. In addition, the param- 

cSlec^ed? atmosphere at the time of acquisition will be 

spectrometer has a 1-milliradian field-of-view, and 
Its spectral range coverage is from 0.4 to 2.4 and 6.2 to 15.5 

spectrometer has a pointing and tracking cipabil- 
degrees forward, 10 degrees aft, and 20 degrees to the 
tJackP^ ground track. The astronaut uses the view-finder/ 

which Le in'^hii'"f-;^?a target sites during data acquisition, 

than a minute. The primary 
'''' magnetic tape and are returned with each crew 

Hultispectral 'Scanner (S-192 ) 

Itl experiment is to assess the feasibil- 

ity of multispectral techniques for remote sensing of earth re- 

simaturf^idpSf??' jci^lly, attempts will be made at spectral 

^ mapping of ground test sites in 

g 1 ulture, forestry, geology, hydrology, and oceanography. 
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The basic instrument design is that of an optical mechanical 
scanner using .an image plane scanning mirror, with a folded 
reflecting telescope used as a radiation collector. The scanner 
operates in 13 spectral intervals of the visible, near-infrared, ‘ 
and thermal-infrared regions of the spectrum ranging from 0.41 ■ 
to 12.5 urn. The primary data are recorded on magnetic tape and 
are returned with each crew rotation. 

The spectral range covered by the scanner overlaps the range of 
the multispectral cameras (S-190) and the IR spectrometer (S-191), 
permitting a cross-check of results deduced from these three . .=■ 

systems- In addition, the IR spectrometer may provide atmospheric- 
density profiles useful for correcting the primary causes of 
atmospheric attenuation of the scanner data. 

Microwave Radiometer /Scatterometer and Altimeter (S-193 ) 

The objectives of this experiment are simultaneous measurement 
of the radar differential backscattering cross section and passive 
microwave thermal emission of the land and ocean on a global 
scale, and engineering data for use in designing radar altimeters. 

The microwave radiometer/ scatterometer experiment is a combination 
of an* active radar scatterometer and passi.ve radiometer. The 
radar backscattering cross section measurement gives a measure 
of the combined effect of the dielectric properties, roughness, 
and brightness temperature of the terrestrial surface. Informa- 
tion over test sites is obtained by the NASA earth resources 
aircraft for validation and extrapolation of spaceborne measure- 
ments. All data are recorded on magnetic tape. 

L-Band Radiometer (5-194 ) 

The experiment objective is to obtain measurements of the bright- 
ness temperature of the earth’s surface along the spacecraft 
track. The L-band radiometer has basically the same operating 
principle as the radiometer part of the microwave radiometer/ 
scatterometer experiment, except that the operating frequency is 
changed from 13.9 GHz to 1.42 GHz. A function of the experiment 
is to supplement the measurement results of Experiment S-193 by 
taking into consideration the effect of clouds on radiometric, 
measurements. By using two frequencies (S-193 at 13.9 GHz and 
S“194 at 1.42 GHz) simultaneously in measurements, corrections 
can be made on radiometric data to include the cloud effects.' 

All data are recorded on magnetic tape'. 

A summary of the SKYLAB EREP sensor characteristics is given in 
Table 2 ' 



OEIGINAL PAGE IS 
OF-pOOE'QUAIZOa 


Table 2, 

SKYLAB 

EREP SENSOR CHARACTERISTICS 


SENSOR 

DESCRIPTION 

SPECTRAL 

COVERAGE 

SPECTRAL 

RESOLUTION 

GROUND 

COVERAGE 

SPATIAL 

RESOLUTION 

S-!90(A) 

MULTtSPECTRAL 

PHOTOGRAPHIC 

CAMERA 

SIX 70mm CAMERA 
MATCHED DISTORTION AND 
FOCAL length (t5,2cm5 
12 METERS REGISTRATION 
18 FILTERS 
21® FOV 

MICROMETERS 
.5 - .6 PANX BiW 
.6- .7 PANX B&W 
.7 - .8 IR B&W 
.8 - .9 IR 8&W 
.5 - .88 IR COLOR 
.4 - .7 HR COLOR 

0.1 MICROMETERS 

163 X 163 Km 

APPROX. 

24m TO 63nri • 

S-t?0(3) 

BARTH TERRAIN 
CAMERA 



460 mm FOCAL LENGTH 
lU mm FILM FORMAT 
3 FILTERS 

0.4 TO 0.7 H. R. 

AERIAL COLOR 
0.5 TO 0.38 IR COLOR 
0.5 TO 0.7 HIGH 
DEFINITION AERIAL B&W 

0.1 MICROMETERS 

• 

109 X 109 Km ' 

APPROX, 
lOm TO 33m * 

S-I9I 

INFRARED 

SPECTROMETER 

POINTED BV CREW 
FILTER WHEEL 
ONE SEC , SCAN RATE 
ONE mRAO FOV 
CRYOGENIC COOLER 
1 6mm CAMERA 

0.4 TO 2.4 AND 
6,2 TO 15,5 micro- 
meter 

TO 4% 

0-45° FV/D 
0-20° SIDE 
0-l0°REAft 

0.44 Km 
SPOT 

S-192 

MULTISPECTRAL 

SCANNER 

IMAGE PLANE SCANNER 
6000RPM SCAN MIRROR 
CRYOGENIC COOLER 
HgCdTe DETECTORS 
^(13 USED) ■ ■ 

'0.186 mRAD FOV 

0.4 TO 2.35 AND 
10.2 TO 12.5 
MICROMETERS 

13 BANDS: 

0,04 TO 0.1 
MICROMETERS 

AS Km 
SWATH 

80 X 80 m 

SPOT 

S-193 

MICROvmE 

RADIOMETER/SCATTEROMETER 
AND altimeter 

l.lm PARABOLIC ANTENNA 
TWO AXIS GIMBAL (0-40<> 
IN FIVE STEPS) 

1.5® FOV 

DUAL POLARIZATION 
ALTIMETER NADIR SEEKER 

13.8 TO 14.0 GHz 
(13.9 GHz CENTER 
FREQUENCY) 

SCAT RECEIVER: 

FIRST IF. 500 MHz 
SECOND IFs SO MHz 
RAD RECEIVER: 

SINGLE FREQUENCY 

0-48 FWO 
0-48 SIDE 

n X 11 Km 
SPOT 

S-194 

L-8AND 

RADIOMETER 

1 m PHASED ARRAY 
(6x8 ELEMENTS) 
COLD and hot P.EF. 

1 .400 TO ) ,427 GHz 

18 MHz FROM 

CENTER FREQUENCY 

111 Km CIRCLE 

Ml Km 
SPOT 


* - DOES NOT INCLUDE LOSS DUE TO ATMOSPHERE EFFEaS OR FILM PROCESSING, 

Table 2. 


NASA hQ ML7I-S75J 
REV. 2-10-73 
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. IDENTIFICATION OF THE NOAA/NOS INVESTIGATION 

The NOAA/NOS proposal to perform analytic aerotriangulation 
utilizing SKYLAB photography has been designated as SKYLAB EREP 
INVESTIGATION NO. 459 and was performed under NASA PURCHASE ORDER 
T-4110B. The task-site identification number for the study is 
9 31651 in which 9 31 classifies the task as being in the CARTOGRAPHY 
MAP ACCURACY discipline, and 651 identifies the test site as the 
CARETS AREA, which runs from Charlotte, North Carolina, northeast 
to Delaware Bay. 

The official NASA description of this investigation is as follows; 

931 Investigate the feasibility of utilizing spacecraft 

(S-190B) -imagery for analytic aerotriangulation methods 
to provide low-order, high-density control network suit- 
able for small-scale mapping applications . 

Employ this technique for coastal zone mapping at medium 
and small scales, surveys- in remote areas, forest and 
range management, various planning activities and route 
location for highways, pipelines, transmission lines, 
and canals. 

SKYLAB S-190B EARTH TERRAIN CAMERA (ETC) PHOTOGRAPHY 

SKY^B S-190B photography was secured over the test' site, during 
orbit 36, on September 12, 1973 (SL-3). The film used in the ETC 
was Aerial-Color, High-Resolution SO-242. Second generation 
transparencies, positive in tone and direction when viewed on the’ 
emulsion side, were made from the original film by printing emulsion- 
to-emulsion in contact. These 1:946,000 scale, 4.5 x 4.5-inch 
transparencies and 3 .x 9 -inch contact paper prints were provided to 
the Coastal Mapping Division of the National Ocean Survey for 
processing through the analytic aerotriangulation system. 

The photography consisted of a strip of 19 photographs ranging along 
a 500-mile track from Charlotte, North Carolina, (frame 86-288) 
to Atlantic City, New Jersey, (frame 86-3'06). A break in the 
required 60 percent overlap reduced 'the usable strip of photography 
to 12 photographs along a 350-mile track from Charlotte, North 
Carolina, (frame 86-288) to the Rappahannock River in Virginia 
( frame 86-299). 

f^lthough the photography provided sharp high-resolution imagery, 
the selection of pertinent images for measurement on the Wild STK 
stereocomparator was hampered by the extensive .cloud cover occurring 
Dn most of the 11 stereoscopic models comprising the analytic strip. 



PHOTO CONTROL 
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Despite the cloud cover,- 29 photo control ■ points of known position 
and elevation were office-identified on the SKYLAB photography. 

Road intersections were located by stereoscopically examining 
the photographs and comparing them with 1:24,000 scale USGS 
quadrangles covering the area. The Geographic Positions were 
scaled from the quadrangles by linear interpolation between the 
2 ' 30” intervals shown on each quadrangle. The scaling was performed 
five times and a mean Geographic Position computed. In addition, 
aeronautical aids to navigation and airport z’unway ends were 
identified on the pictures, and their positions and elevations 
determined from data secured by the Coastal Mapping Division under 
its Airport Obstruction Chart Survey program. The office-identified 
road intersections and aeronautical aids provided images slightly 
superior in quality to that of the office-identified airport 
runway ends. The locations of these photo control points or 
stations on the ETC photographs is shown in Figure 1. Table 3 
describes the stations and their approximate accuracy. 

All of the photo control stations were at least 1/4-inch in from the 
sides of the 1:946,000, 4.5 x 4.5-inch transparencies. Twenty-'five 
of the stations appeared on lonly two consecutive overlapping photo- 
graphs in the strip, while :^our stations appeared on three consecu- 
tive photographs. The location of the camera clock within the 
photographic format prevented two control stations from creating 
imagery on three consecutive overlapping pictures. 

PASS POINTS 

Two pass points were established in the nine conventional relative 
orientation locations on each photograph. The pass points were 
drilled into the transparency emulsion with the Wild PUG-2 stereo- 
scopic point transfer device in areas providing an optimum stereo- 
scopic perception. As a consequence of the intrusion of the camera 
clock into the photographic format , it was necessary to set the 
pass points along the bottom easterly edge of the strip at least 
1/2-inch in from the sides of the transparencies. 

The two pass points were used in the preliminary analytic aerotri- 
angulation programs consisting of the Image Coordinate Refinement 
and Three-Photo Orientation program, the Secant Plane Coordinate 
Transformation program, and the Strip Adjustment to Ground Control. 
However, only one of the two pass points in each relative orientation 
location was used in the block adjustment of the strip. 

MARKING AND PHOTOCOORDINATE MEASUREMENT 

The Wild PUG-2 stereoscopic point transfer device was used to select 
and mark only the pass points by drilling holes into the photographic 
emulsion at these images using a 60-micron diameter diamond-tipped 
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TABLE 3 . ' 

ACCURACY OF OFFICE IDENTIFIED CONTROL 

USED IN 

BLOCK ADJUSTMENT 


CONTROL 

APPROXIMATE 

APPROXIMATE 

STATION 

HORIZONTAL . 

VERTICAL 

NUMBER 

288100 
288101 ■ 

ACCURACY 

(METERS) 

■ ACCURACY 
(METERS) 

293100 
294102 ■ 
235100 

15 / 

5.0 • 

299100 



288110- 

288111 

290110 

296111 

288201 

288202 

290201 

290111 

5 

3.0 

292110 

292111 

293110 

296201 

296110 

298110 

299110 

299111 

5 

uo 

• 

0 

288120 



291120- 

- 


291121 
293120 
293121 
297120 ’ 
297121 

1 . 

0.3 


DESCRIPTION 


AERONAUTICAL AIDS 
Horizontal and vertical 
positions from Airport 
Surveys, N0S« 


ROAD INTERSECTIONS 
Horizontal and vertical 
positions from 1:24,000 
scale USGS quadrangles. 


ROAD INTERSECTION SPOT 
ELEVATIONS 

Horizontal and vertical 
positions from 1:24,000 
scale USGS quadrangles. 


CENTERLINE RUNWAY ENDS' 
Horizontal and vertical 
positions from Airport 
Surveys 5 NOS. 
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drill. The photo control point images were not' drilled in order 
to preserve the sharpness of the imagery. 

The measurement of the x and photocoordinates for the pass points 
and control stations was performed on a Wild STK stereocomparator. 
Use of the stereocomparator allowed the operator to drill the pass 
point images down the center of each photograph only and to then 
stereoscopically transfer the drilled image to the overlapping 
photograph for measurement, without actually .ever drilling the 
image on the second photograph. The stereocomparator measuring 
mark consisted of a 165-micron diameter black circle having a 
20->micron black dot at its center. The dot was centered in the 
60-micron diameter drilled pass point image when observing the 
photocoordinates for the point, 

FIDUCIAL MARKS 

In the time interval occurring between the film exposure, 
development, and the subsequent printing of the glass plate diapos- 
itive or transparency, the aerial film undergoes a random enlarge- 
ment and shrinkage change. Since the accuracy of analytic photo- 
grammetric computations depend on the use of a true central perspec- 
tive, it is necessary to compensate for the film distortion and 
thereby mathematically return the film to the physical format 
present at the instant of exposure. This can be achieyed by 
comparing the positions of the images created by the fiducial marks 
on each photograph with the true positions of these marks in the 
camera focal plane. The photograph is then mathematically stretched 
so as to place the fiducial marks back into their true positions. 

In metric mapping cameras, the fiducial marks are located in the; 
corners of the camera focal plane. Some cameras have additional 
fiducial marks at the midpoints of the sides of the focal plane. 

The marks are normally a part of the lens cone and thus remain in 
a fixed position relative to the camera lens. The intersection 
of the diagonals joining the corner fiducial marks should represent 
the principal point of the photograph, i.e., the foot of the 
perpendicular from the focal plane to the nodal point of the camera 
lens. The corner fiducial marks on the Wild aerial cameras owned 
by NOS consist of an interrupted cross having a 100-micron diameter 
dot at the center. 

The ETC- is not a metric camera in the photogrammetric sense because 
the fiducial marks and the image frame are a part of the removable 
film magazine and hence are not in a fixed position relative to the 
camera lens. Fortunately, the resulting lack of precision in 
locating .the principal point on the photography is of minor conse- 
quence in narrow angle cameras, such as the ETC. 
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The ETC has a series of holes drilled around the perimeter of the 
image frame. These holes created photographic images having an 
approximate diameter of 330 microns. The images were of poor 
quality and rather ragged around the edges. For this reason, the 
stereocomparator operator centered the 165-micron circle of the 
measuring mark in the center of four holes selected to serve as 
fiducial marks. 

NOS normally employs flash plates to provide a photographic record 
of the true relative positions of the camera fiducial marks. The 
flash plates are made in the laboratory by exposing a diapositive 
mounted in the camera so that its emulsion lies in the camera focal 
plane. Because the emulsion is secured on a stable glass base, the 
coordinates of the fiducial marks can be measured on a comparator 
later, with no concern for film shrinkage distortion. 

Since no flash plate was available for the ETC, a nominal set of true 
fiducial coordinates was obtained by mounting each of the transpar- 
encies in turn on the comparator and then reading the photocoordi- 
nates for the four selected fiducial holes on each of the photos. 

The data were entered into -a flash-plate-reduction program to accom- 
plish the following tasks: 1. Correct the observed photocoordinates 

of the fiducial holes for comparator systematic errors, 2. determine 
by least square methods a meaned set of nominal true photocoordinates 
for the fiducial holes in a coordinate system, having its origin at th 
principal point (intersection of diagonals joining fiducials 1—3 and 
2-4) and oriented so that = y^ = y 2 * See Figure 2 below. As a 

consequence of this computation, the direction of flight becomes the 
x-axis of the photoooordinate system. 


Figure 2, Transparency emulsion- side up 


comparator 

operator 



/ 


flight 

direction 


' 2 . 




comparator x-axis 
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The meaned set of nominal true photocoordinates for the fiducial 
holes provided' by the flash-plate-rcduction program are: 


FIDUCIAL HOLE 


X (microns) y 


1 

2 

3 

4 


59227 . 74 
59112.61 
-59202.38 
-59112.77 


50154.18 

-50132.70 

-50132.70 

50132,84 


The data were then prepared for processing through the HOAA/NOS 
analytic aerotriangulation system of computer programs. All of 
the computations were performed on the CDC 6600 computer. 

COMPUTER PROCESSING 


The 12-photo strip extended over the three states of South 
Carolina, North Carolina, and Virginia. The computational 
processing requires all of the ground positional data to be 
expressed in a common three-dimensional coordinate system. In 
order to attain this condition, and also to compensate for the 
presence of earth curvature in the data, the initial computation 
was to develop secant plane ^coordinates for each of the 2 9 office- 
identified control stations.’ Accordingly, the Geographic Posi- 
tions and elevations of these points were processed through the 
Secant Plane Coordinate Transformation program. 

Secant' Plane Coordinate Transformation : The elevations of the 

control points obtained from the 1:24,000 scale USGS quadrangles 
and airport surveys are based on sea level and thus do not recog- 
nize the existence of earth curvature. The program computations 
begin with a conversion of the Geographic Positions and elevations 
to an orthogonal geocentric coordinate system-.having its origin 
at the center of the earth as defined by the Clarke 1866 Spheroid. 
The- geocentric coordinates are then transformed into a secant plane 
coordinate system in which the secant plane intercepts the earth's 
surface near the edges of the area to be mapped so that most of 
the terrain objects will possess a positive Z (elevation) coordi- 
nate. The origin of the secant plane system is placed near the 
center of the project area. The secant plane origin selected 
for the SKYLAB study was Latitude 36® 20' 00", Longitude -78® 45' 00 
The Z-axis is the extension of the normal to -the ellipsoid which, 
because of the earth's ellipsoidal nature, does not pass through 
the center of the earth. The X-axis points towards the East and 
the Y-axis points towards the North. , • 

Table 4 shows the Geographic Position and elevation input to 
the program and the resulting secant plane coordinate output 
from this program. 
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CONVERSION OF CONTROL STATION GEOGRAPHIC POSITIONS 
TO SECANT PLANE COORDINATES 


SECANT 

PLANE ORIGIN 





LATITUDE 36 20 00 LONGITUDE • 

-078 45 00 


GEOGRAPHIC POSITION INPUT 












Elevation ( 

288100 

34 

59 

20,200 

-080 

57 

18.000 

650.0 

288110 

35 

15 

25,015 

-081 

01 

37.254 

710.0 

288201 

35 

15 

05,A.24 

-081 

01 

43.630 

633.0 

288111 

34 

44 

15.094 

-030 

40 

58.811 

565.000 

288202 

34 

43 

34,559 

-080 

42 

15.428 

598.0 

288120 

35 

12 

31.865 

-080 

57 

00,375 

707.0 

288101 

35 

13 

13,390 

-080 

56 

18.073 

740.00 

290110 

35 

44 

30.052 

-080 

21 

43.729 

760.00 

290201 

35 

44 

34.185 

-080 

21 

48.169 

753.0 

290111 

35 

03 

03 .-25 8 

-079 

50 

57.400 

456.0 

291120 

36 

05 

57.511' - 

-079 

57 

09.510 

926.0 

291121 

36 

05 

13.059 

-079 

56 

15.037 

900.0 

292110 

36 

08 

32.741 

-079 

47 

26. 181 

843.0 

292111 

35 

36 

54.282 

-079 

02 

53.990 

185.0 

293100' 

35 

52 

20.520 

-078 

47 

01.000 

430.00 

293120 

35 

51 

52.243 

“078 

47 

51.245 

398.0- 

293121 

35 

51 

57.031 

-078 

46 

49.346 

401,0 

293110 

36 

05 

18.889 

-079 

03 

46.083 

516.00 

294102 

36 

40 

29.600 

-079 

00 

53.200 

530.00 

295100 

36 

49 

04-400 

-077 

54 

11 .700 

350.00 

296111 

37 

13 

22.840 

-077 

59 

24. 112 

322.00 

296201 

37 

1'3 

10.548 

-0 77 

59 

37.971 

339.0 

296110 

36 

39 

51.186 

-077 

34 

17.309 

156.0 

297120 

37 

29 

44-. 350 

-077 

18 

25.908 

161.0 

297121 

37 

30 

14.556 

-077 

18 

38,599 

160.0 

298110 

37 

12 

17.161 

-076 

59 

13 .914 

112 .0 

299100 

37 

26 

54.700 

-076 

42 

42. 100 

10. ,00 

299110 

37 

54 

13.346 

-076 

51 

44.738 

031 .0 

299111 

37 

26 

06.654 

-076 

20 

11.559 

013.0 

SECANT 

PLANE 

OUTPUT in METERS 







X 

T 

Z 

288100 

-201278.692 

-146855.736 

-3874,837 

288110 

-207169.956 

-116985,134 

-3424,735 

288201 

-20 7344.107 

“117584.396 

-3464.941 

288111 

-176997.192 

-175263,639 

-3901.313 

288202 

-178969.743 

-176473.122 

-3979.759 

288120 

-200294.019 

-122480.950 

-3309.711 

288101 

-199196,924 

-121226.254 

-3241,238 

290110 

“145819.163 

-64433,577 

-962. '699 

290201 

“145928.545 

“64304.337 

-966.024 

290111 

-100287,221 

-141707,209 

-1430.977 

291120 

■ “108297.906 

-25295.089 

110.910 

291121 

-106952.115 

-26681.925 

119.998 

2921-10 

“93656.711 

-20679.874 

333.841 

292111 

-27029.441 

-79651,423 

297.541 

293100 

“3035.511 

-51147.971 

721.944 

293120 

“4296.418 

-52018.862 

704.400 

293121 

-2743.370 

-51871.929 

707.371 

293110 

-28172.927 

-27112.861 

834.668 

294102 

"23669.361 

37934.083 

80 1 . 864 

295100 

75550.344 

54100.713 

226.922 

296111 

674A8.740 

98990.886 

■' -231.39,4 

296201 

67110.169 

98609.407 

-216.716 

296110 

105357.883 

37358,925 

-134.052' 

2 97120' 

127578.566 

129930.492 

-1756.009 

297121 

127252,667 

130856.770 

-1768.812 

298110 

156468.136 

98126.586 

-1843.117 

299100 

180328.210 

125649.347 

-2988.252 

299110 

165981.068 

175874.338 

-3784. 177 

299111 

213537.475 

124932.816 

-3997.981 
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Image Coordinate Refinement and Three-Plioto Orientation : The 

X and y photocoor-dinates observed on the stereooomparator for 
the fiducial and nonfiducial images on each transparency were 
then processed through the image coordinate refinement and 
three-photo orientation program. All of the photocoordinates 
were first corrected for comparator calibration errors. The 
program then performed a least squares fit of the fiducial hole 
photoGoordinates to the nominal true photocoordinates for 
these fiducial holes, as previously obtained from the fla'sh- 
plate-reduction program. This operation of placing the fiducial 
holes back into their true positions serves to correct all of 
the data for film shrinkage distortion and to express the photo- 
coordinates in a two-dimensional photocoordinate system having 
its origin at the principal point and oriented so that the 
x-axis is the direction of flight. 

The systematic errors still remaining in the photocoordinate 
data are those due to the distortions introduced by the aerial 
camera lens and atmospheric refraction. The Addendum to LEG/ 

ASD Technical Memo No. TM 73-002 - April 1973 issued on July 11, 
1974, indicated the camera lens symmetrical distortion and the 
asymmetrical distortion-* caused by lens decentration to be 
probably insignificant from zero. In addition, distortion due 
to atmospheric refraction at camera altitudes above 40 miles 
is relatively negligible. For this reason, no attempt was 
made during the image coordinate refinement phase to compensate 
for camera lens and atmospheric refraction distortion. 

The refined image coordinates provided by the program theoreti- 
cally should be nearly all free of systematic error and contain 
'only residual observational discrepancies in them. These 
refined coordinates were then punched out to serve as input to 
the block adjustment program. 

The program then proceeded to the three-photo camera orienta- 
tion phase, which comprises an interrelated geometric fitting 
of the photographs, based only on the refined photocoordinates 
and is entirely independent from any ground ^control data. The 
computation is iterative and derives the orientation of each 
photograph relative to the previous two in the strip. It also 
determines the positions of all pertinent objects in a three- 
dimensional coordinate system at the scale of the photography. 
The collineation principle is imposed in a least _ squares solu- 
tion that minimizes the residual observational discrepancies 
in the image coordinates. The residual discrepancies are ^ 
analyzed by the computer, which discards those images exhib- 
iting excessively large discrepancies. The removal of these 
blunders provides "clean" photocoordinate data for all sub- 
sequent computations. 
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StTrip Adjustment to Groimd Control : The analysis of three photo- 

graphs at a time automates" the joining of the separate triplets 
into a continuous strip and develops a set of mode'l coordinates 
that are analogous" to the product obtained from conventional 
-stereotriangulation on stereoscopic plotting instruments. The 
horizontal and vertical strip adjustment transforms the model 
coordinate data into the prevailing ground control coordinate 
system, which is a secant plans coordinate- system for this study, 
by fitting to .the control stations uhrough the application of 
polynomial equations and least squares. Any large residual 
discrepancies appearing in the resulting adjustment are cor- 
rected in order to obtain blunder-free provisional ground 
position data prior to entering the block adjustment computation. 

The strip adjustment of the SKYLAB photography was performed 
holding to the IM- photo control stations identified on Figure 1 
by a A . Note: These same 14 -stations were employed later to 

control the block adjustment solution. The strip adjustment 
was performed twice — going from frame 86-28.8 to frame 86-29 9 and 
then going from frame 86-295 to frame 86-2 88,. In both adjust- 
ments, the resulting discrepancies at the- 14 held photo control 
stations and 15 withheld stations increased in magnitude from 
about 25 meters at the beginning of the strip to approximately 
100 meters at the tail end of the. bridge. Results of this 
nature -do not occur on normal photogrammetric m.apping operations 
conducted by NOS,. The appearance ‘ of these apparently, systematic 
strip adjustment discrepancies on the SKYLAB bridge is assumed 
to be attributable to the failure to completely compensate for . 
the systematic errors introduced, into the data by the nonmetr.ic 
characteristics of the Earth Terrain Camera, 

I 

Block Adjustment: - In order to maximize the accuracy of the 
analytic aerotriangulation , the block adjustment program was 
applied, using the previously obtained refined photocoordinates 
and the provisional object coordinates. The program permits a 
simultaneous solution of the absolute orientation of all the 
photographs, together with a determination of the finalized 
coordinates for each object. This office has developed three 
simultaneous analytical aerotriangulation block adjustment 
programs for operation on the CDC 6600 computer. 

1, 25-Photo Block Adjustment: This program was designed to 

service smaller organizations not having access to large-size 
computers and will accommodate blocks up to 25 photos in 
size. All input /output is on cards, and the program requires 
less than 50,000 words of computer core storage. The logic 
of the solution. is similar to that of the 185-photo block 
adjustment program. - ‘ • 
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2. 185-Photo Block Adjustment: The block can contain as many 

as 185 photographs in a single simultaneous least squares 
solution. All input /output is on tape. Approximately one 
million words of storage are required, and therefore auxiliary 
disk storage and extended core storage is used to augment 

the CDC 6600 core memory. 

The area to be block adjusted may be of any shape. The strips 
of photography can be of variable length and may have any 
overlap with other' strips in the block. Thus, diagonal cross- 
flights may be included, if desired. The photographs can be 
entered into the solution in any order. Photographs taken 
by aerial cameras having different focal lengths may be used 
simulxaneously in the solution. 

All of the pass points and control stations contribute equa- 
tions to the normal equation matrix and thus influence the 
least squares orientation solution. Corrections to the provi- 
sional coordinates for these objects are computed simultane- 
ously with the determination of the absolute orientation of 
all the photographs in the block. The unweighted control . 
stations perform as if they are pass points. The weighted 
control stations can be computed as if they are pass points 
by using their refined | photocoordinates and the finalized 
■camera parameters from' the orientation solution to determine 
‘their ground coordinates by intersection. 

3. 600-Photo Block Adjustment: Blocks up to 600 photographs in 

size can be accommodated in this version. All input /output 
is on tape. Auxiliary disk storage is necessary because the 
program requires nearly one million words of memory. 

This version is not as flexible as the 185-photo program in 
terms of data input organization, but the arithmetic approach 
employed results in a more efficient computation and a much 
shorter computer running time. The area should be square or 
rectangular in order 'to simplify the arrangement of input data. 
Photographs must be entered into the solution in the exact 
order in which they were taken and may not overlap each other 
by more than 60 percent. No strip can have more than 20 
pictures in it. The program permits the mixing of photographs 
'taken by aerial cameras of different focal lengths. 

All of the pass points and control stations contribute equa- 
tions to the normal equation matrix and thus influence the 
least squares orientation solution. The finalized camera 
parameters from the orientation solution and the refined 
photocoordinates of the pass points and control stations are 
used to compute the final ground coordinates "for these 
objects by intersection. 
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The 600-photo block adjustment program was .used for the SKYLAB 
analytic aerotriangulation study. A Fortran listing of this 
program is given in Appendix B. 

Thirty-six pass points (one in each relative orientation loca- 
tion) and all of the 29 office-identified photo control stations 
were permitted to contribute observation equations to the normal 
equation matrix and thereby influence the least squares orienta- 
tion solution of the 12 photographs comprising the block. The 
provisional coordinates for these objects’ should be reasonably 
close to their true values in order to minimize the number of 
iterations required of the block adjustment solution. For this 
reason, the initial ground (secant plane) coordinates of the pass 
points consisted of the data furnished by the first half of the 
strip adjustment going from frame 86-288 to frame 86-299 — and by 
the first half of the strip adjustment going from frame 86-299 
to frame 86-288. The known true ground (secant plane) coordinates 
were used as the initial coordinates for the 29 photo control 
points. 

Table 5 is a listing of the initial provisional ground (secant 
plane) coordinates for the pass points and the photo control 
stations. A listing of the refined photocoordinates for these 
pass points and photo control stations, as previously punched 
out by the Image Coordinate Refinement and Three-Photo Orienta- 
tion program, is given in Table 6 . 
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292291310 

292291320 

292291330 

292292310 

292292320 

292292330 

292293310 

292293320 

292293330 

292292110 

292291120 

292291121 

292293100 

292293110 

292293120 

292293121 

292292111 

293 

293292310 

293292320 

293292330 

293293310 

293293320 

293293330 

293294310 

293294320 

293294330 

293293100 

293293110 

293293120 

293293121 

293294102 

294 

294293310 
294293320 
294293330 
294294310 
2942 94320 
294294330 
294295310 
294295320 
294294102 
294295330 

295 

295294310 

295294320 

295294330 

295295310 

295295320 

295295330 

295296310 

295296320 

295296330 

295296110 

295296111 

295295100 

295296201 

296 

296295310 

296295320 

296295330 

296296310 

296296320 

296296330 

296297310 

296297320 

296297330 

296296110 

296296111 

296295100 

296296201 


-4.8126509E-02 
-5.2868542E-02 
■-5.1129797E-02 
1 .8499420E-03 
4.0858396E-03 
-9.3830666E-04 
5.1925419E-02 
5.0031725E-02 
4.7535343E-02 
-1.9237138E-02 
-3.4734189E-02 
-3.4488084E-02 
3. 8201346 E-02 
3.2417993 E-02 
3.6554509E-02 
3. 7984083 E~02 
-7.2191183E-04 

-4.6876203E-02 
-4.4675254E-02 
-4.9701058E-02 
3.2954850E-03 
1.3644901E-03 
-1 .1609361E-03 
4.9291420E-02 
5.0331744E-02 
5.0443036E-02 
■1 .0472546E-02 
'1.6304848E-02 
■1 .2119203E-02 
-1.06S4280E-02 
2 .9769123 E-02 

-4.5408740E-02 
-4.7413524E-02 
-4.9922832E-02 
6.4136652E-04 
1 .6050855E-03 
1.7489357E-03 
4.8336748E-02 
4.5921897E-02 
-1.9005334E-02 
5.0752810E-02 

-4.8065487E-02 
■4,7116111E“02 
■4.6971411E-02 
•2.8582516E-04 
•2.7288240E-03 
2.1121530E-03 
4.8350200E-02 
4. 8859741 E-02 
4. 3942 351 E-02 
4.1888299E-02 
4.9557802E-02 
2.7383244E-02 
4.9025404E-02 

•4.8912550E-02 

■5.1417141E-02 

■4.6558160E-02 

■2.1747625E-04 

2.4754045E-04 

■4.6724768E-03 

5.1814312E-02 

5.I438535E-02 

4.9087131E-02 

■6.6781630E-03 

9.5107825E-04 

•2.1251919E-02 

4.1528880E-04 


-3.0527771E-02 
4. 7874533 E-02 
3.5506342E-03 
“3.0333683 E-02 
4.8972518E-02 
7.85324026-03 
-2.8029944E-02 
5.0058976E-02 
4.5751374E-03 
4.3086281E-02 
4.B658235E-02 
4.6586076E-02 
-4.1701161E-02 
-4.9433929E-03 
-4.1634943E-02 
-4.2517091E-02 
-5.0470942E-02 

-3.'1929435E-02 
4.7409829E-02 
6.2704290E-03 
-2.9647102E-02 
4.8547056E-02 
3. 0063077 E-03 
~2.6523827E-02 
3.7760516E-02 
4.6554986E-03 
-4.332a930E-02 
-6.5265715E-03 
-4.3255058E-02 
-4.4135354E-02 
4.6062768E-02 

-3.1276435E-02 
4.6910495E-02 
1 .3757361E-03 
-2.8164236E-02 
3.6179394E-02 
3.0448055E-03 
-2.7948864E-02 
3.1223507E-02 
4.4462580E-02 
7.0171463E-03 

-2.9755754E-02 
3.4612072E-02 
1.4643224E-03 
-2.9563473E-02 
2.9684546E-02 
5 .4490202E-03 
-3.0367873E-02 
4.8204153E-02 
9.7309548E-03 
-4.1150655E-02 
3.5889388E-02 
-7.6541133E-03 
3.5791637E-02 

-3.1175464E-02 
2.8070525E-02 
3.8367575E-03 
-3 .2005718E-02 
4.6660238E-02 
8.1401176E-03 
-2.1334985E-02 
4.2232885E-02 
5.1457554E-03 
-4.2801992E-02 
3.4330292E-02 
-9.2674902E-03 
3.4232754E-02 
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297296310 

297296320 

297296330 

297297310 

297297320 

297297330 

297298310 

297298320 

297298330 

297296111 

297297120 

297297121 

297298110 

297296201 

298 

2982 97310 
298297320 
298297330 
298298310 
298298320 
298298330 
298299310 
298299320 
298299330 
298299110 
298297120 
298297121 
298298110 
298299100 
298299111 

299 

299298310 

299298320 

299298330 

299299310 

299299320 

299299330 

299000301 

299000302 

299000303 

299299110 

299299100 

299299111 


-^,88846838-02 
-4. 847143 9E-02 
-5.3337673E-02 
3.2292043E-03 
2.8155666E-03 
.4.7887551E-04 
4,9364170E-02 
4.8950281E-02 
5.3156090E-02 
-4.7770507E-02 
2.3627514E-02 
2.3945413E-02 
2.7641953E-02 
-4.830775 9E-02 

-4.5428380E-02 
-4.591988aE-02 
-4.8222649E-02 
8.0468845E-04 
2.9431491E-04 
4.5390987E-03 
4. 7350772 E-02 
4.4515357E-02 
5.1835518E-02 
3.7430048E-02 
-2.5071781E-02 
-2.4755306E-02 
-2.0994424E-02 
1.7160147E-02 
4.4993911E-02 
1 1 

-4.7827856E-02 
-4. 8340 1 6 7E -02 
-4.4103047E-02 
-1.1510283E-03 
-4.0500174E-03 
3.2847107E-03 
-1.9679956'E-ro 
-1.9679956E-10 
-1.9679956E-10 
-1.1156657E-02 
-3.1453155E-02 
-3.5411964E-03 

End of !Pabls 6, 


-3. 354083 8E-02 
4.5150300E-02 
6.6214000E-03 
-2.2914424E-02 
4.0743480E-02 
3.6049954E-03 
-3.1794614E-02 
4.5322814E-02 
1.1645544E-02 
3.2823481E-02 
2.0265087E-02 
2.1263032E-02 
-2.557631 2E-02 
3.2723607E-02 

-2. 44704 11 E-02 
3.9161207E-02 
2.0387070E-03 
-3.3377516E-02 
4.3796670E-02 
1.0100176E-02 
-5. 1665726E-02 
4.8243494E-02 
1.2330690E-02 
3.3063014E-02 
1.8706945E-02 
1.9706280E-02 
-2.7139170E-02 
-1.9058386E-02 
-4.1062900E-02 

-3.4851692E-02 

4.2357455E-02 

8.6372888E-03 

-5.3232239E-02 

4.6813308E-02 

1.0854518E~02 

’-’3.8'6237'08E-09 

-3.8623708E-09 

-3.8623708E-09 

3.1610337E-02 

-2.0549835E-02 

-4.2610240E-02 
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Weighting the B'lock Adjustment Solution : The weighting of the 

block adjustment solution is performed by applying image quality 
and control station weights to the data during the computations. 
These weights can be defined as follows: 

.Image Quality Weights: It is logical to weight a block adjust- 

ment in favor of those observation equations provided by the 
better quality images because their equations are more reliable. 
Image quality is primarily influenced by lens resolution and 
■ the type of ground object creating the image. The weighting 
is accomplished by multiplying each observation equation by a 
number expressing its relative reliability. 

Control Station Weights : The observation equations are written 

for each image on every photograph in the block created by the 
pass points, control stations, and other objects which are used 
to influence the least square orientation solution. When 
written for image created by the control stations, it is neces- 
sary to recognize that their initial provisional X, Y, Z ground 
coordinates were obtained by classical ground surveying methods 
and should be favored during the block adjustment. This is 
accomplished by increasing the size of the main diagonal elements 
of the normal equation matrix, which are the coefficients of 
the unknown dX, dY, and/or dZ correction terms. This serves to 
reduce the size of- the unknown dX, dY, and/or dZ correction 
terms when the normal equations are solved- . By .reducing the 
magnitude, of the corrections to the initial approximations, the 
least square adjustment is constrained in favor of these ini- 
tial values. Control stations that are not subjected to this 
type of weighting perform as pass points. These unweighted 
or withheld control stations provide' a means for evaluating 
the accuracy of the block adjustment solution. 

Presently, empirical values are used, for these weights at NOS 
instead of being based on the standard error of the observa- 
tions as required by rigorous statistical methods. Also, the 
present NOS block adjustment program multiplies the pertinent 
normal equation main diagonal terms by the control station 
weights instead of adding on- a number to increase their' size. 

The resulting accuracy of the block adjustment can be expressed 
as a photogrammetric RMS error and a geodetic RMS error. The 
photogrammetric RMS error is computed using the residual v^ and 
Vy plate observational discrepancies of all images contributing 
observation equations to the orientation solution. The geodetic 
RMS error reflects a comparison of the X, Y, Z results of the 
block adjustment computation with the known X,’y, Z coordinates 
for the control stations that were obtained by classical ground 
surveying methods. In most NOS operations, weights are selected 
that cause the photogrammetric RMS error and the geodetic RMS 
error to be about equal. This has the effect of providing for 
an equal distribution of the block adjustment errors between 
the photogrammetric observations and the geodetic field observations 



36 


The inherent errors in using nonmetric S-190B photography and 
office-identified photo control made it necessary to perform 
numerous block adjustment solutions involving different combina- 
tions of control and weights. The best results were achieved by 
using 14 weighted control stations distributed uniformly along 
the perimeter of the strip as shown in Figure 1. In NOS mapping 
operations using metric photography and field-identified photo 
control, eight weighted photo control stations would normally 
have been sufficient for a block adjustment of the 12 photos. 

As previously noted, the pass points were drill holes in the 
emulsion and not images of specific terrain objects. The stereo- 
comparator operator cannot' remove parallax exactly on the drilled 
pass point holes when observing their photocoordinates. This 
reduces the reliability of the pass point images during the photo- 
coordinate measurement process. The office-identified photo control 
stations, on the other hand, were prominent ground features 
providing sharp images on the overlapping photographs. Since 
drilling was not necessary for these stations, the comparator 
operator was able to remove parallax directly at the images before 
observing their photocoordinates. As a consequence of their 
higher reliability, a larger image quality weight was assigned 
to the control station images. 

Many combinations of image [quality and control station weights 
were applied to the 14 held photo control stations in an effort 
to optimize the accuracy of the block adjustment solution. In 
general all of these weight combinations yielded a horizontal 
position geodetic RMS error of approximately 15 meters for the 
15 withheld (unweighted) photo control points. The maximum 
horizontal position error on any withheld control station was 
less than 26 meters. 

Results of the Block Adjustment Solution : The weights used for 

each of the 14 held photo control stations in the final block 
adjustment of the SKYLAB data were: image quality weight = 6j 

control station weight for X and Y = 6; control station weight 
for ^ = 3. A smaller weight was used for Z because of the limited 
accuracy in the stereo height determination resulting from the 
low base-height of 0.10, even though the known elevations of the 
control stations were of a higher accuracy than the known horizontal 
positions for these stations (see Table 3). In fact, a higher 
Z weight was found to degrade the block adjustment accuracy. 

The pertinent output from the 600-photo block adjustment program 
for the SKYLAB study is shown in Table 7 . A summary of the 
residual errors remaining at each of the 29 control stations is 
given in Table 8. All of the ground coordinate data and the 
residual errors are in meters and expressed in the secant plane 
coordinate system. The results presented here are from a block 
adjustment solution in which the standard error of the ground 
control was assumed to be 4.1 meters. The residual errors given 
in Table 8 are also shown on Figure 1 in red. 
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PHOTO CONTROL 

STATION X 


Y 


Z 


A 288100 
A 288110 
A 288201 
A 288111 
A 288202 
A 288120 
A 288101 
A 290110 
A 290201 
A 290111 
A 291120 
A 291121 
A 292110 
A 292111 
A 293100 
A 293120 
A 293121 
A 293110 
A 294102 
A 295100 
A 296111 
A 296201 
A 296110 
A 297120 
A 297121 
A 298110 
A 299100 
A 299110 
A 299111 


11.738 
10.025 
10.991 
-7.517 
-10.718 
14.378 
14.961 
-2.030 
3.3.01 
1.143 
. -0.207 
' 12.644 
-1.150 
0.491 
-4.988 
-2.891 
8.579 
1-3,274 
1.034 
15.142 
1,938 
-1.008 
-4.219 
3.231 
1.667 
-6.028 
-16.286 
-3.290 
0.296 


.-12.483 

20.323 

7.727 

-5,612 

-2.895 

9.686 

19.031 

6.814 

-2.078 

-6.379 

6.891 

-6.815 

6.593 

-3.279 

12.935 

-0.096 

7.711 

9.672 

1.407 

-12.953 

0.900 

6.180 

-0.994 

7.052 

2.042 

-3.617 

-18.695 

3.970 

-1.523 


-30.343 
-233.104 
-196.716 
172.285 
134.556 
-205.982 
-210.726 
11. 736 
■40.588 
-86.241 
-78.019 
-25.212 
-38.472 
-5.3-78 
25,196 
-3.173 
-5 . 563 
-13.012 
2.332 
-70.711 
-7.411 
3.373 
29.125 
-137.767 
-135.449 
85.999 
-138.440 
6.269 
-26.832 


Note: A = weighted photo control station 

A - unweighted photo control station 


Table 8 . Residual errors in meters remaining at each of the 
computed positions for the 29 office-identified 
photo control stations after block adjustment 
solution, as expressed in the secant plane 
coordinate system. 


The block adjustment program is designed to terminate the iter- 
ative computation when the computed corrections to all of the 
angular camera parameters are each less than 0.00001 radians 
(two arc-seconds). Table 7 shows that five iterations of the 
block adjustment orientation solution were required to achieve this 
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condition fon the SKYLAB study. Usually, only one such pass 
through the solution is necessary in conventional NOS mapping 
projects employing metric photography and field-identified 
photo control. 


The residual errors at the control stations appear to be 
uniformly distributed throughout the test area, and there is no 
evidence that the least square solution was not able to absorb 
uncompensated _ systematic error, i.e. , no large 'isolated discrep- 
ancies exist in the solution. The horizontal- position geodetic 
RMS error for the 29 photo control stations was 12.218 meters 
and is equivalent to 12.915 microns at the SKYLAB photography 
scale of ^ 1 : 946 ,000 . The geodetic RMS error computed for only 
the ^15 withheld photo control stations was 15.068 meters. The 
maximum horizontal position error was 24.794 meters and occurred 
at withheld station No. 299100. No serious attempt was made 
to hold ^closely to the elevations of the control stations because 
of the inherent limited accuracy in the stereo height determina- 
tion. Consequently, several of the residual errors in Z exceeded 
-200 meters. ~ 

The photogrammetric RMS error was 12,996 microns at plate scale 
and was computed using the residual v and v plate discrepancies 
of ^ all the images created on the photography^by the 36 pass 
points and the 29 photo control stations. It should be noted 
that the photogrammetric RMS error is usually about eight microns 
on conventional NOS photogrammetric mapping projects. 

Inverse of the Secant Plane Coor dinate Transformation: After 

completion of the block adjustment solution, the adjusted 
computed secant plane ^ coordinates were transformed back into the 
original ground coordinate system (Geographic Positions and 
elevations based on sea level) by applying the secant plans trans- 
formation in its inverse mode. The results of this inverse 
computation are displayed in Table 9. 



TaBie 9, 

SBCANT PIANE COOROIMATES - OUTPUT OF BLOCK ADJUSTMENT Iff METERS 


PASS POINTS 

X 

y 

Z 

2S8310 

-175159..941 

-178654,007 

-3783.068 

2883Z0 

-224533.782 

-118960.309 

-4350.547 

288330 

-196938.725 

-154397.051 

• -3873.622 

289310 

-140968.327 

-150574i057 

-2564.325 

289320 

-186866.970 

-95485,277 

-2601. 182 

289330 

-163941.845 

-124266.567 

-2443.416 

290310 

-106548,364 

-122833.069 

-1039.850 

290320 

-152383*610 

-63563.072 

-1164.540 

290330 

-131685.875 

-96569.496 

-899,891 

291310 

-73506.450 

-91626.878 

, “629,. 069 

291320 

-121335.359 

"36081.66? 

-444.274 

291330 

-94988*940 

-67988.199- 

-315.295 

292310 

-36515,376 

-63241.457 

449.438 

292320 

-79693.386 

-3100.307 

'295.041 

292330 

-60177.549 

-36467.791 

522.021 

293310 

-557,116 

-33201.994 

880.547 

293320 

-46155.21,1 

23748.689 

755.439 

293330 

-22306.105 

-11480.969 

1059.436 

294310 

31898,162 

-4934,554 

690.663 

294320 

-3694,629 

43458.533 

304.189 

294330 

15111,681 

18897.434 

550.016 

295310 

67366,805 

22135.331 

82.479 

295320 

32090.883 

64818.444 

330,914 

295330 

49365.290 

49544.761 

285.874 

296310 

104040.087 

49038,235 

87.815 

296320 

59972,357 

107751,460 

-•348.662 

296330 

78054,828 

76365.232 

94.902 

297310 

136781.254 

86353.340 

“1159.746 

297320 

100539. .722 

133424.766 

-1563.397 

297330 

11976X..767 

104510-.673 

-T230.9-94 

298310 

176204.799 

105752,504 

-2617.263 

298320 

132263.976 

162905,948 

-3090.611 

298330 

154437.310 

140222.139 

-2944.173 

299310 

221333.928 

118339,172 

-4290.205 

299320 

162683,562 

191X83.725 

-4075.060 

299330 

188418.605 

168589.801 

-4134.662 

14 HELD 

CONTROL STATIONS 


288110 

-207159.931 

-116964.805 

-3657,839 

2881H 

-177004.709 

-175269.251 

-3729,028 

290110 

-145 821.193 

-64426.763 

-950.963 

290111 

-100286.078 

-14 1713. 5 38 

-1517.218 . 

292110 

-93657.861 

-20673.281 

295,369 

292111 

-27028.950 

-79654.702 

292.163 

293120 

-4299.309 

-52018.958 

701.227 

294102 

-23668.327 

37935.490 

804.347 

296110 

105353.664 

37357.931 

-104.927 

296111 

67450.678 

98991.736 

-238.805 

297120 

127581.797 

129937.544 

-1893.776 

298110 

156462.108 

98122.969 

-1757.118 

299110 

165977,778 

175878.308 

-3777,908 

299111 

213537.179 

124931.293 

-4024.813 

15 WITHHELD CONTROL STATIONS 


288100 - 

-201266. 9S4 

-146868.219 

-3905.180- 

288101 

-199181.963 

-121207.223 

-3451.964 

288120 

-200279.641 

-122471.264 

-3515.693 

288201 

-207333.116 

-117576.669 

-3661,657 

288202 

-178980.461 

-176476.017 

-3845,203 

290201 

-145925.244 

-64306.415 

-925.436 

291120 

-108298.113 

-25288.198 

32,891 

291121 

-106939.471 

-26688.740 

94.786 

293100 

-3040.499 

-51135.036 

747.140 

293110 

-28159.653 

-27103.189 

821.656 

293121 

-2734.791 

-51864.218 

701.808 

295100 

75565.486 

54087.760 

156.211 

296201 

67109.161 

98615.587 

-213.343 

297121 

127254.334 

130858.812 

-1904.261 

299100 

180311.924 

125630.652 

- 3126-697 



GEOGRAPHIC POSITIONS - OUTPUT OF BLOCK ADJUSTMENT 


PASS POINTS 



288310 

•34 

latitude: 

42 26.38609 

- ^ 

80 

longitude: 

39 43.87202 

Elevation ffeet) 
1096.07686 

I 

288320 

35 

14 

6.91221 


81 

13 

2.76273 

-279.18950 


288330 

34 

55 

18.68739 

- 

80 

54 

20.37638 

796.27163 

c 

289310 

, 34 

57 

57.28077 

“ 

80 

17 

37.42388 

-70.44538 

289320 

' 35 

27 

17.25727 

- 

80 

48 

31.79789 

176 .23997 


289330 

35 

11 

58.24487 

- 

80 

33 

1.74416 

259.67279 

( 

290310 

• 35 

13 

13,48098 

- 

79 

55 

13.16051 

783.04996 

290320 

35 

44 

54.64260 

- 

80 

26 

5.64609 

572.11904 


290330 

■ 35 

27 

14.74808 

- 

80 

12 

2.23300 

1293 .58307 

( 

291310 

35 

30 

16.84703 

- 

79 

33 

36.97538 

-1124.96987 

291320 

36 

0 

1.95660 

- 

80 

5 

44.95660 

44.75680 


291330 

35 

42 

57.32927 


79 

47 

59.30026 

-139.43936 


292310 

35 

45 

45. 64825 


79 

9 

13.58220 

233.10629 

292320 

36 

18 

7,55350 

- 

79 

38 

14.22851 

-13.95916 


292330 

36 

0 

10. 10432 

- 

79 

25 

2.77197 

370.14007 

■ c 

293310 

36 

2 

2.81021 

- 

78 

45 

22.25265 

557o'=^6828 

293320 

36 

32 

46,46501 

- 

79 

15 

55.66983 

555 .24755 


293330 

36 

13 

46.60473 

- 

78 

59 

53.15875 

1021.67941 

( 

294310 

36 

17 

18.00507 

- 

78 

23 

41.73789 

-82 .35886 

294320 

« 36 

43 

29.95754 


78 

47 

28.89599 

-1127.14234 


294330 

’ 36 

30 

12.68101 

- 

78 

34 

52.74101 

-660.65353 


295310 

^ 36 

31 

49.66979 

“ 

77 

59 

5 1.81921 

-1053.04192 

L 

295320 

! 36 

55 

.94303 

- 

78 

23 

23.53281 

-181 .68449 

» 

295330 

36 

46 

42.80926 

- 

78 

11 

49.19933 

-418. 66557 

( 

296310 

36 

46 

10.50551 

- 

77 

35 

4.84073 

1073 .23836 

296320 

37 

18 

8.87137 

- 

78 

4 

24.81194 

159.90532 


296330 

.37 

1 

5.88875 

“ 

77 

52 

22.45203 

765 .15397 

: c 

297310 

37 

6 

6.14785 

- 

77 

12 

40.16500 

309.91260 

297320 

37 

31 

49.55 76,3 

- 

77 

36 

44.86165 

-554.28613 


297330 

37 

16 

3.46029 

- 

77 

23 

58.88423 

-151.24184 

t 

; < 

298310 

37 

16 

12. 17141 

- 

76 

45 

46.80372 

-338 .95441 

298320 

37 

47 

32,17061 

- 

77 

14 

53.24438 

-1411.50456 


298330 

37 

35 

4.08519 

- 

77 

0 

4.26612 

-1072 .32185 

^ ( 

299310 

37 

22 

26.46691 

- 

76 

15 

1.24412 

-487,89597 

299320 

38 

2 

32.00796 


76 

53 

47,32371 

248 .66309 

1 

' r 

i 

299330 

37 

50 

1.81199 


76 

36 

33.11741 

277.17958 

14 HELD CONTROL STATIONS 






( 

288110 

35 

15 

25.53695. 

_ 

81 

1 

37.17181 

-56.53923 


288111 

34 

44 

15.06445 

- 

80 

40 

58.91790 

1130.99762 

( 

290110 

35 

44 

30.27598 


80 

21 

43.80370 

798.41836 

290111 

35 

3 

2,98861 

- 

79 

50 

57.40463 

173.57077 


292110 

36 

8 

32.95024 


79 

47 

26.25236 

716.77843 

c. 

292111 

35 

36 

54.17346 

- 

79 

2 

53.97097 

167.48599 

293120 

35 

51 

52,23900 

- 

78 

47 

51.36030 

387.59883 


294102 

36 

40 

29.64527 

- 

79 

0 

53. 15813 

538.16070 


296110 

36 

39 

51.15007 

“ 

77 

34 

17.49878 

251.29327 

( 

296111 

37 

13 

22.87241 

- 

77 

59 

24.02989 

297.80444 


297120 

37 

29 

44.66704 

- 

77 

18 

25.65847 

-290.11984 

■ ( 

298110 

■ 37 

12 

17.00557 


76 

59 

14.24742 

393.36137 

299110 

37 

54 

13.47135 

- 

76 

51 

44.87622 

51 .63104 

! ^ 
1 

i 

299111 

37 

26 

6.62123 

“ 

76 

20 

11.53555 

-75.09492 

15 WITHHELD 

CONTROL STATIONS 




! ^ 

288100 

34 

59 

19,78010 

— 

80 

57 

17.56339 

550.25723 


288101 

■ 35 

13 

13.88282 

- 

80 

56 

17.75518 

46,38242 

i 

288120 

35 

12 

32.05602 

- 

80 

57 

.06754 

29.57440 

! ^ 

288201 

35 

15 

5.55994 


81 

1 

43.45203 

-13.58050 


288202 

34 

43 

34.58196 

- 

80 

42 

15.70143 

1040.36233 

C 

290201 

35 

44 

34.13316 

- 

80 

21 

47.99968 

885.94229 

291120 

36 

5 

57.72387 

- 

79 

57 

9.5 7444 

669.99247 


291121 

36 

5 

12.83931 


79 

56 

14.54521 

816.69433 

c 

293100 

35 

52 

20.94618 


73 

A7 

1.19852 

512.32645 

293110 

36 

5 

19.20237 

- 

79 

3 

45,55598 

472.98178 


293121 

35 

51 

57.27979 

- 

78 

46 

49.00424 

382.53015 


295100 

36 

49 

3.99475 

- 

77 

54 

11.05976 

118.26105 

: 

296201 

37 

13 

10.74701 

- 

77 

59 

38.01136 

350.34280 


297121 

37 

30 

14.71046 

-* 

77 

18 

38.41834 

-283.95531 


299100 

37 

26 

54.19126 


76 

42 

42.61679 

-446.64844 


■Rnfl of O 
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Accuracy Analysis : In order to evaluate fully the accuracy 

potential of the analytic system, a final computer program is 
used to develop the inverse of the matrix of normal equations) 
the variances, and the standard' errors in X, Y, and Z at all of 
the points throughout the project area. 

It can be assumed that the camera parameters and the ground 
positions of the .pass points provided, by the final pass through 
the block adjustment solution V70uld not change significantly 
should an, additional pass.be made. Thus the same refined image 
coordinates, together with the final camera parameters and ground 
positions for the pass points, will yield essentially the same 
normal equations that occurred in the final block adjustment pass. 

This is the basis for the accuracy-analysis program. 

/■ 

The standard error E of the coordinates determined at a point in ' 
the block can be expressed as E = Qe where ^ is the variance of 
the point as derived from the inverse, and ^ is the standard 
error of unit weight for the problem and is considered to be 
essentially equal to the photo gramme trie RMS error determined 
in the block adji^stment solution. Both Q and e are relatively 
Independent and provide a means for the comparison of tests 
conducted under varying conditions. 

The variance' Q is affected by the geometry of the block, such as 
the number of photographs and the .number and distribution of 
horizontal and vertical control. Its value can be computed 
from simulated photographs before the pictures are actually taken 
and is unaffected by poor techniques. The standard error of unit 
weight e is a measure of the precision of the system and is 
affected by the camera, comparator, effectiveness of the correc- 
tions for systematic errors, overlap, premarking, operational 
techniques, etc. Its value is relatively constant for a given 
set of techniques and allows one to upgrade the system by 
improving, the techniques. Thus, for example, the nonmetric 
characteristics of the Earth Terrain Camera and the use of office- 
identified photo control resulted in a standard error of unit ' 
weight of nearly 13 microns for the SKYLAB study. , This is 
significantly larger than the eight microns or less that is found 
in NOS operations employing metric aerial cameras and field- 
identified photo control. 

Table 10 shows the horizontal standard errors in meters in the 
secant plane coordinate system for each of the 15 withheld 
(unweighted) photo control stations. The horizontal position RMS 
error computed from these values is 16.414 meters and substantiate 
the validity of the geodetic RMS errors found in the previous 
block adjustment solution. 



PHOTO 



CONTROL 



STATION . 

X meters 

Y meters^ 

288100 

• 9-,.'201 

9.529 

288201 

9. 961 

15.053 

288202 

10.164 

12.360 

288120 

9,003 

13.679 

288101 

8.916 

14.029 

290201 

7.271 

10.185 

291120 

13.554 

16.164 

291121 

13.158 

15.701- 

293100 

15.650 

11.110 

293121 

15.736 

11.260 

293110 

9.967 

9.494 

295i00 

9.377 

10.147 

296201 

7.949 

7.844 

297121 

10.949 

8.993 

299100 

9.084 

13.678 

RMSj^ = 

10.963 RMSy = 

12.2.16 


RMS = 16.414 meters 
YY 


Table 10. The standard errors in meters in the secant plane 
■ coordinate system for 'each of the IS withheld 
(unweighted) photo control stations. 



53 


FURTHER DISCUSSION OF SKYLAB ANALYTIC AEROTRI ANGULATION RESULTS 

In evaluating the results of the SKYLAB aerotriangulation 
study, it must be remembered that these results were achieved 
using a strip of photography instead of a block of overlapping 
strips of pictures. For the case of a strip, analytic compu- 
tations are usually terminated after strip adjustment because 
there is no evidence of a significant improvement in accuracy 
by continuing on through block adjustment. However, the strip 
adjustment of the SKYLAB photography appeared to show apparently 
systematic adjustment discrepancies that were assumed to be 
attributable to a failure to compensate completely for the 
systematic errors introduced by the nonmetric characteristics 
of the camera and/or the office-identification of photo control. 
For this reason, the block adjustment program was applied in 
an effort to optimize the accuracy of the aerotriangulation 
solution. 

The results obtained from the block adjustment were reasonably 
close to the values to be expected from the SKYLAB photography. 
Our experience indicates that the block adjustment of a strip 
of metric 1:946,000 scale photography using field-identified 
photo control would yield a geodetic RMS error of approximately 
14 meters. This figure must be modified to allow for the 
additional errors introduced by the nonmetric characteristics 
of the Earth Terrain Camera (focal plane shutter, camera lens 
distortion assumed to be negligible, imprecise location of 
the photograph principal point, etc.) and the office-identifi- 
cation of photo control. Assuming a maximum error of 20 meters 
introduced by the ETC and a maximum error of 15 meters for the 
office-identified photo control, the overall expected geodetic 
RMS error for the block adjustment of the SKYLAB strip of 
pictures increases to nearly 16 meters. As noted in -this report, 
the actual geodetic RMS error was 12.218 meters for all 29 
office-identified photo control stations and 15.068 meters for 
the 15 withheld or unweighted photo control stations. 

The National Standards of Map Accuracy require 90 percent of 
all map points to be correct to within 1/50 inch or 0.51 mm. 
for maps published at scales of 1:20,000 or smaller. Statis- 
tically, the SKYLAB results indicate that 90 percent of all 
29 office-identified photo control stations were held to 
within 20 meters, and 90 percent of the 15 withheld or unweighted 
stations were held to within 24.7 meters. It is evident, 
therefore, that if the positions of all the planimetric detail 
required to construct a map of the project area were developed 
digitally by analytic block adjustment methods, 90 percent of 
these planimetric points would also be correct to within less 
than 25 meters. Thus, the analytic aerotriangulation method 
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can _ be used in this manner with the 1:946,000 scale SKYLAB 
strip photography to construct a 1:50,000 scale map that will 
meet the National Standards, of Map Accuracy. 

The usual practice in mapping operations is to compile the 
planimetric details from stereoscopic models oriented to 
horizontal -photo control established principally by analytic 
aerotriangulation procedures. Experience has shown that 90 
percent of the horizontal photo control should be known to 
within 0.15 mm,, as measured on the manuscript. This is 
equivalent ^to 24.75 meters at a map scale of 1:165,000. Thus, 
stereocompilation techniques can be combined with the analytic 
aerotriangulation methods to construct a map at 1:150,000 to 
1:200,000 scale from the 1:964,000 SKYLAB strip photography. 
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THE SKYLAB S-190B EARTH TERRAIN CAMERA 


John D, McLaurin 
U.S. Geological Survey 
McLean, Va. 22101 

One of the major objectives of the SRylab manned space station to 
be launched in 1973 is to collect earth resources data. The Earth 
Resources Experiment Package (EREP) of Skylab includes instruments for 
collecting data in several regions of the electromagnetic spectrum, 
ranging from the visible to the microwave. The sensors include the 
S-190A multispectral camera, the S-191 infrared spectrometer, S-192 
multispectral scanner, the S-193 microwave system, and. the S-19A L- 
band radiometer. These systems have all been described elsewhere 
(NASA-MSC, 1971) . However, a new sensor, the S-190B Earth Terrain 
Camera (ETC), has recently been added to EREP. Because the ETC is not • 

as well known as the other sensors, I will describe its characteristics 

and indicate some potential applications of ETC photographs. 

The ETC was included in the EREP of the Skylab mission as an 

addendum to the S-190 experiment. It is designed to supply high- 

resolution photographs of areas within the field of view of the other 
EREP "sensors to aid in the interpretation of data gathered by them. 

In soma cases information from the ETC photographs will substitute 
for groynd truth and for photographs obtained from .aircraft underflights 
Furthermore, the resolution of the camera will permit certain investi- 
gations that would be impossible for -any of the other EREP sensors 
alone . 

The EROS (Earth Resources Observation Systems) Program of the 
Department of the Interior is interested in the ETC primarily because 
it has approximately the same photograph scale and ground resolution 
as the film-return satellite camera system which was recommend in 
1967 by the National Academy of Sciences for cartographic and photogram- 
metric applications . That proposed satellite would be in a near-polar 
orbit at an altitude of 200 km and would include a metric camera of 
300 mm focal length. In 1970 the Department of Interior proposed to 
NASA that a satellite of that type be flown, and we have a continued 
interest in it. 

The ETC is a modified version of the Lunar Topographic Camera 
carried on the Apollo 13 and 14 missions. It is being built by Actron 
Industries, Inc. (formerly Hycon) under contract to NASA, The body is 
• an extensively modified KA— 74 reconnaissance camera body with a focal- 
plane shutter and vacuum film flattening. The lens has a focal length 
of 460 mm, a fixed aperture of f/4, color correction, and maximum 
radial distortion of 10 pm. Forward image-motion compensation is 
provided by rocking the entire camera in its mount during the exposure. 





The frame format is 115 mm by 115 mm so that at the'Skylab altitude 
the format covers an area of 109 km by 109 km. Characteristics of 
the camera can be summarized as follows : 

Lens — 460 mm focal length, f/4 fixed aperture color corrected 

Lens distortion — Radial, ± 10 vira; tangential, ± 5 pm 

Shutter— Focal plane, bidirectional; 1/100, 1/140, 1/200 sec. 

Forward-motion compensation — By rocking camera, 0 to 25 mrad/sec. 

Film — 125 mm (5 in.), 2-mil. base; 450 frames/roll 

Format — 115 x 115 min 

Framing rate — 0 to 25 frames/min. 

Overlap — 15% Standard; 0 to 80% availabli 

'Ground coverage at nadir — 109 x 109 km. 

The camera has a control box with a switch for ' selection of manual 
or automatic operation. The forward-motion compensation system can be 
set to operate within the range of 0 to 25 mrad/sec, and the framing 
rate from 0 to 25 frames /min. Three shutter speeds are available, 

1/100, 1/140, and 1/200 sec. Figure 1 is a photograph of the camera. 

The ground resolution of the camera depends on the film used. The 
three films being considered are SO-242 high-resolution color, 3443 color 
infrared, and 3414 high-resolution black-and-white. To obtain estimated 
ground resolution for the different films, Actron has run computer . 
simulations that model the forward-motion compensation system, the 
attitude-error rates of the spacecraft, the shutter speed, the lens 
characteristics, and the film and filter characteristics. Table 1 
summarizes the simulations (the shutter speeds were subsequently 
changed). As can be seen, the expected ground resolution varies frbm 
10 to 39 meters per optical line pair. In addition to the films listed, 
it is possible that a color infrared film of higher resolution will be 
available for at least one of the Sky lab missions. • . 

The ETC will be mounted’ in the Scientific Airl.ock of the orbital 
•workshop of sky lab. The other EREP- sensors wilT be located in the 
Multiple Docking Adapter. Figure 2 presents a view of the Skylab 
space station with the major components indicated. The ETC will be 
boresighted to record the same ground areas that the other EREP' sensors 
are viewing. 

Skylab will be operated as four missions (fig. 3). The first 
mission, Skylab 1, will launch the unmanned space station. The next 
day the first three-man crew will be launched in an Apollo spacecraft 
as mission Skylab 2, which will reoulre the crew to occunv the space . 



TivBLE 1» Predicted ETC resolution, in meters on the ground 
per optical line pair. 


Case 

Film 

Shutter 

speed 

(sec) 

High- • 
contrast 
resdlut ion 
(1000:1) 

Low 

contrast 

resolution 

(2:1) 

1 

3443 

1/100' 

21 

39 

2 

3443 

1/200 

21 

38 

3 

3443 

1/500 

21 

38 

4 

3414 

1/100‘ 

8 

15 

5 

3414 

1/200 

6 

11 

6 

3414 

1/500 

5 

10 

7 

SO-242 

1/100 

12 

22 

8 

SO-242 

1/200 

11 

20 

9 

SO-242 

1/500 

11 

20 

Based on 

> 

1 computer simulations by Actron Industries 

, Inc. , July 1971, 

station 

for about -28 days. The crew 

will conduct 

several experiments. 

including the EREP series. About 2 months after the return of the 

Skylab 2 

! crew, the Skylab 3 crew will be launched 

to occupy the space 

station 

for as long as 

56 days, again conducting a 

variety of experi- 

ments . 

The final mission, Skylab 4, 

will start about 1 month after 

the return of the Skylab 3 crew and will also have 

a duration, of up to 


56 days. 

Skyiab will follow a 435 km circular orbit with an inclination of 50® 
.which will carry the station over any portion of the earth between 50°N 
and 50®S latitude. The normal attitude of the space station is referred 
to as the solar-inertial mode, required by the solar panels and the 
heating constraints. For the EREP experiments the space station must be 
maneuvered into the z-local— vertical mode. Because of thermal and other 
constraints, the mode can only be used for a limited niimber of orbits. 
Current plans call for approximately 60 z-local-vertical passes. 

The orbital and mission constraints will thus limit the number and 
location of EREP data-collection passes ^ it will not be possible to 
collect data for large contiguous areas, and repetitive coverage of an 
area will be limited. Nevertheless, considerable data of value to 
various earth resources investigations should be collected during the 




Figure 2, — The Skylab spacecraft. 

Bkylab program. For. the Skylab 2 mission, 4 rolls of film containing 
450 frames each will be available, and 6 rolls will be available for 
aach of the other 2 missions. As many as 7200’ ETC photographs may be 
obtained. ... 

The exact areas of the United States where ETC photographs will 
be acquired have not yet been selected. Generally, the final areas 
for photographic coverage will be determined during the mission. Investi- 
gators whose EREP proposals were accepted by NASA have been notified 
of the addition of the ETC to Skylab, and many of them have requested 
ETC photographs of their test areas. In addition, many Federal and 
State agencies have requested ETC photographs of specific areas. The 
requests are being coordinated with the Skylab mission planners in an 
effort to take photographs of as many areas of interest as orbits, 
weather, and other constraints will permit. ETC photographs will be 
available to the public at nominal cost through the EROS Data Center 
of the U.S. Geological Survey, at Sioux Falls, S. Dak. 

■ The design of the ETC will limit its applications. First, the 
ETC is not a metric camera in the photogrammetric sense. Because the 
image frame is a part of the removable film magazine and because of the 
use of a focal-plane shutter, the geometric quality of the photographs 
is limited. The principal point cannot be precisely located, and 
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Figure 3,-"7Sequence of Sky lab missions. 

-therefore analytical applications will be limited. The ETC has a 
limited field of view, 14°. When the camera is operated for 60% overlap 
the base-height ratio is only 0.15; thus, the use of the ETC for 
stereoscopic height determination will be especially limited. 

In spite of the limitations, the ETC represents a significant 
advance in camera systems for earth resources obseirvations from space. 
The ground resolution is considerably better than that of any camera 
previously used. A recent paper by Colvocoresses (1972) compares the 
image resolution of ERTS, Skylab, and Gemini/Apollo space photographs. 
Table 2, compiled from data in that paper, summarizes the ground 
resolution of the various systems. From the tabulated data, it is 
obvious that the ETC has ground resolution from 3 to 20 times better 
than the other space photographic systems. Moreover, the ETC fills 
the gap between the other space systems and high-altitude aircraft, 
cameras, which normally have ground resolution of i m or better. 

The ETC will also permit a comparison of multispectral and multi- 
spatial data collection. The S-190A multispectral camera will provide 
narrow-spectral-band photographs useful for multispectral interpretation 
The ETC, on the other hand, will provide photographs of a different 


A 
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TABLE 2. Comparison of Ground Resolution for Space Imaging Systems. 
Ground resolution given in terms of the photographic 
criterion of optical line pairs, in meters on the ground per 
line pair. 

System 

High-contrast 
(100:1 or 1000:1) 

Low-contrast 
(2:1 or 1.6:1) 

ERTS-A 

BBV, green band 

L26 

180 

RBV, red band 

126 

180 

RBV, infrared band 

156 

275 

MSS 

244 

316 

Sky lab 

S-190A 



Hlgh-resolution film 

22 

. 38 

Low-resolution film ' 

49 

99. 

S-190B (ETC) 



High-resolution film 

10 

15 

Low-resolutlon film 

20 

•38 

Gemini/Apollo 

High-resolution film 

50 

70 

.Low-resolution film 

80: 

125 

scale and resolution, which can be compared with- the 
and thus provide an evaluation of multi spatial data. 

S-190A photographs 


The S-190A and ETC photographs may also be used with aircraft 
photographs for multistage sampling, a technique that starts with the 
interpretation and classification of small-scale photographs of a large 
region. Interpretations are made on progressively larger scale 
photographs of smaller and smaller areas within the large region. The 
application of the technique in forestry has been described by 
Langley (1969) . 

Therefore, the principal applications of ETC photographs will be 
in experiments where high resolution is required. Each individual 
photograph will be a nearly orthographic view of the ground, with 
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rather low distortion within the frame. The high resolution will 
greatly benefit experiments in v/hich photointerpretation is important. 

An example of the kind of experiment planned for .the ETC 
photographs is photomapping at 1:250,000 and 1:100,000 scale. The 
image scale of the ETC photographs will be about 1:945,000. Doyle (1971) 
has proposed criteria for the resolution required for photomapping and 
the useful enlargement of the photographs: 


where 


R 10~^S„ 

g ,1^ 

R = required ground resolution (m/lp) 

'O 

Sjjj = map scale number. 


The suggested criterion for photographic enlargement is expressed as 


M = 12. 
^ 10 


where rp = original photo xesolution (Ip /mm) 

M = allowable enlargement from photograph scale. 

& 

According to criteria, the required ground resolution for 1:250,000 
and 1:100,000 scale photomaps is 25 m and 10 m. Assuming the use of 
3414 film, the approximate photograph resolution is 80 Ip/mm and the 
allowable enlargement is 8X'. Enlargement to 1:250,000 scale would 
require 3.8X and to 1:100,000 would require 9.5X. At 1:100,000 scale 
the image would still have a theoretical- resolution of 8 Ip/mm, which 
may be satisfactory in the practical sense. 

The U.S. National Map Accuracy Standards (NMAS)- require that '90% 
of the well defined points tested be no more than 0.5 mm from their 
correct position at map scale. For a 1:250, 000-scale map the tolerance 
converts to 125 m while for a 1.: 100, 000-scale map it is 50 m. A 
computer program prepared by DBA Systems, Inc. , for the Geological 
Survey has been used to determine how much relief can be tolerated 
before planimetric image displacement exceeds NMAS. The program 
includes the effects of earth curvature, atmospheric refraction, 
terrain relief, location of the image in the photograph format, and 
map-projection scale factor. The computer analysis indicates that 
about 500 of relief can be tolerated at the extremes of the usable 
photo format for the photograph to meet the standards for 1:250, 000- 
scale mapping. For mapping at 1:100,000 scale, only about 300 m of 
relief can be tolerated. The UTM was used as the map projection in 
the analysis. 


Conditions which could significantly affect the positional 
accuracy of the images, however, are the effects of the focal-plane 
shutter. and of errors in the forward motion compensation (FMC) system. 
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A NASA study (McDermltj 1971) considered the effects of spacecraft 
residual rates, FMC errors, earth rotation, shutter type, and spacecraft 
rigidity. The study concluded that, in the worst case, errors due to 
the sources considered v/ould amount to about 35 ym displacement between 
the leading and trailing edge of an image. That is, the dimension of a 
discrete image will be changed 35 ym in the direction of motion. Much 
of the error could be reduced by proper operation of the camera, but a 
random component of about 22 ym would probably remain, equivalent to 
21 m on the ground. Thus, there is some question whether the 1:100,000 
positional accuracy requirement can be met. Maps at the scale of 
1:250,000, however, appear to be well within the capability of the camera. 
The uses plans to conduct photomapping experiments at both scales in 
order to determine the usefulness of space photographs of ETC resolution. 

- Other mapping experiments planned for the ETC photographs include 
map revision and thematic mapping. Experiments will be conducted to 
determine the types of map revision information that can be derived from 
the photographs and applied to maps at scales of 1:24,000 and smaller. 

Thematic mapping consists of the preparation of maps depicting such 
data as vegetation distribution, surface-water distribution, snow cover, 
and the massed works of man. Thematic mapping experiments at scales of 
1:250,000 and 1:100,000 are planned, using color infrared photographs 
as the most suitable input for this kind of mapping. 

ETC photographs will also be used for land-use mapping; urban 
development studies; sediment loads and dynamics of San Francisco Bay; 
geological synthesis of the Colorado Plateau; study of hazards and 
tectonics in the Cascades volcanoes; marine geology of the Pacific 
Northwest; and geologic studies of areas in California, Oregon, Oklahoma, 
and a portion of the Great Plains. The experiments will depend largely 
on photointerpretation of the ETC photographs. Some will require normal 
color photographs while others will require color infrared photographs. 

One of the primary objectives of the experiments will be comparison of 
ETC photographs with aircraft photos, S-190A photos, and ERTS-A images. 
Thus, a better assessment can be made of the scale and the resolution 
that are optimum for each particular investigation. The results of the 
experiments will be of great value in defining future data-collection 
systems for earth resources. 

In conclusion, the Earth Terrain Camera provides an opportunity 
to acquire high-resolution space photographs for the study of mapping 
and for investigations of earth resources. The camera fills the gap 
between high-altitude aircraft photographs and other hitherto random 
space photographs. Although the camera has several limitations from 
the photogrammetric standpoint, it will supply high-resolution photographs 
of value to several disciplines. 
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Appendix B 



PROGRAM 


B6C0 


COC 663 0 FTN V3.D-32'* OPT=? D'./02/7'5 20.40.54 


e >5 

t 
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.15 


20 


25 


PROGRAM B6D0 tlNPUT.OUTPJr 

C BLOCK AOJUSTNENf'"ES5A-33 
_C_ MAXIMUM NUK0ER OF PHOTOS I 
C MAXIMUM NUH0ER OF PHOTOS [ 
_ OIHENSION At2G9.259».G(123 

IPtlS.lC) iE <3,563. TITueca). 

OIMENSION IN0tXlt6i30l,INDi 

COHHON/CDiJUF/LEN,H£xr,IFU 
_ COHKON/OOOUF/LtNG.NEX ,IFU 

CALL LTRI0(5LTAPtl,lUB,3, 
CALL LTRI0l5LTAPE2,1160,Qi 
CAUL LTRIO(5LTAPE3, 1160.3. 

CALL OPENHSdl, INDEX! ,60D, 

CALL "OPENH j < 12, INOEX2 .60S , 
□ ATA(IEF=lOOOOCOG011DOrOOJ . 
1 <IP=2A01200536000000C033 

OATA<XSTOP=lH J 
rzYXH = 0 

CALL PATCH _ „ 

IPdZYXH.EO.d) GO TO 999 

PRINT 98 

'98 FORHAT<»l JOB STEP ABORTED 

PRINT 99, 13UFF (II . 1=1. 132 

99 FORMAT ( IX, 8A10) 

lF<(jST.AND.ieFI .NE.DI S3 

CALL LTRIO ( I Q,4B XX XX, 
999 IZTXH = 1 


.TAP£1,TAPE2,TAPE3 .TAPEll.TAP 
AST ANDGEOUETi: SURVEY 1969 
S 0LOCK = 600 
S A STRIP = 20 

3 6) ,CAH(600,12I , iJ( 3.1801 ,C(3, 
3(6.91 

<2(6001 

3 T, IXUUFTOUFFdO 24T,I(3, JSf 
, IBUF,BUF(10 24), IB.ENDFLCt 
3 tKS) 

Q.KSI 

3,KSI 

3 I _ 

j , * “■ 

( ENDFLQ=0 .3 , 

3) 


- INPUT AREA »XI 
4 I 


El 21 

(K EL LERI 


860 ) ' 


B6D0S 

0600 

B600 


0600 ' 

- ’ 5 ' 

A600 

1 

0600 

8 

96CCS 

2 

0600 

9 

0600 

10 

060 0 

u 

8600 

12 

B600 

■ 13 

D600S 

3 

860 OS 

■ ~ 4 

0600 

17 

8600 

18 

060 02 

i 

8600 

19 

06130 

20 


TO 100J3_ 
"JSTl" 


30 


lOOO IF (£NOFL(3,EQ.O.) CALL LF2IO IIP, 1158, XX ,X X, KSTI 
ENOFLQ = 1. 

ITERAT=1 

LCNG =0 $ NEX=IFIR =1 

CALL OaUF (1) 

ENCODE (80,571,OUF(IOUFII 
571 FORMAT (19HU JOO STEF 


B600 
0600 
060 0 
8600 
Q60G 
0600 
0600 
0600 
0600 
8600 
0600 
8600 
060 0 


XSTOP 
■P l,Ail 


35 


INPUT PHASE 


"1(31 H 


40 


45 


50 


55 


KD = KC = LINE = ME 
IGO = IRE = KOISK = 1 
READ 531, TITLE 
IFITIILEdl .EQ.XSTOPICALL 
PRINT 569, TITLE 

PRINT 5 32 

IPHO = JOYCE=HO,RT=0 
LLL=259 

READ BLOCK CONSTANTS AND 
READ 524, LARRY, FL, HT712 ,HT 
IF(JERRY)1,1,2 

PRINT 567 

GO TO 3 
PRINT 568 
K=i 

PRINT 576,FL,HAX,RESID 
PRINT 577,WT712,HTi2P 
PRINT 573 


ER“ LOS '5" 0" 
DBUFC 01 


CONTROL HEIGHTS 
12P,HTCON,HAX,R£SIO,JERRY,INV 


ER 


21 

22 

23 

24 

25 

26 
27' 
20 

29 

30 

31 

32 

33 


0600 

34 

B690 

35 

0600 

36 

0600 

37 

0600 

38 

8600 

■ 39 

8600 

40 

8600 ■ ‘ 

41 “ 

0600 

42 

8600 

43 ■ 

0600 

44 

B600 ■ 

45 ‘ 

0600 

46 

B600 

47 

A600 

4 

■ 0600 ■ 

49 

0600 

50 

“ ■ 8600 

51 

860j 

52 

8600 

53 

B600 

54 

860 0 

55 

8600 

56 
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program B60 0 . 


CDC &&30 FTN V3.1)-324 0PI=2 _ _CA/02/75 _ 20 .GO . 5<*. P AGE 


%% 

^ Q 



Q 


.60 


65 


ei 

, ^ 


70 


75 


60, 


«5 


4 R EAP 55 8 , A ( 1 1 Klj AJ.2 ,K ) . A t3 j KJj msj 

IS=A(1,K) 

IF(JERRY.NE,0I GQ_T0.321 _. . .. 

' NEIU'HTCON 

JOYC£ = A <2t<l - ... 

HORT=A(3,K) 

GO TO 322 

321 NEIL=A(2,KI 

IFlNilL.EQ.O) NEIL=A13.K» 

IF(A(2, Kl .NE.0.0) J0YCE = Hr30N 

IF(A(3, K> .NE.O.OI HORT=HT:ON 

“322 PRINT 574, IS, NEIL, JOYCE. MDU 

JOYCE =H OR r = 0 

IF(ITEST) 5,5,6 

5 K=KU _ . . . 

IF(K-25C) 4,4,27 

__ 27 PRINT 5 70 „ 

STOP 

^ J=1 

"do 7 1=2, K 

IF(Atl,I) .GE.A(l,I-ll_)GO I3_7___ .. 

J=-l 

DO 8 L=l,3 _ _ 

SAYE=AIL,I> 

A (C,I ) = ATL,I-y ■ 

A (L.I-l 1 = SAVE 

a CONTINJE . 

7 CONTINUE 

IFU.LT.QIGO TO 6 

C ■ ' READ IN ALL GROUND COORDINATES FOR BLOC< 

LEN=0 $ NEXT=IFIRST=1 

10=240120053400000000003 
CALL OBUFCGl 
CALL OBUFU) 




95 


iOO 


1 o“5 




j=l _ _ 

9 MCON = J45 

NCOH = NCON-2 

CALL ORUF(l) 7~ 

DEC0DEl8D,527,6UFF(IXQUFn (G(I) , 1= J ,NC0H) , IT EST 
00 15 L=liK 

IFCGt Jl -A <i,L I) 11, 12. 10 

10 CONTINUE 

GO TO 11 . 

'iZ G( Jf4l = A 12.L) ‘ 

G(Jf3)=A<3.L) 

GO TO 13 

11 G(J+4) = G(J+S)=fl . 

13 IFIIIEST) 14,14,15 

14 J=J*-6 

GO TO 9 

15 NC0H=NC0N/6 

00 22 J=l,256,3 
00 22 1=1.207 
A(I,J)=0. 

22 CONTINUE _ 


B&flO 

57 


B6C0 

58 


6600 

59 


6600 

60 


8600 

61 


6600 

62 


O6C0 

63 


SqCO 

64 


6600 

65 


0600 

66 


8600 

67 


6600 ' 

68 


A6U0 

5 


B600 

69 


B60C 

70 


B604 

71 


B600 

72 


0600 

73 


B6DC 

74 


a&oo 

75 . 


3600 

76 


6600 

77 


e&oo 

78 


8600 

79 


6603 

80 


BbOa ~ 

81 


8603 

62 


B603 

83 


B600 

84 


B60 0 ‘ 

AS ■ 


B&OQ 

86 


8600 

■ 37 


B600 

88 


B6G3 "" 

89 


8600 

go 


B600 

91 


6600 

92 


B6uO 

93 


B60J 

94 


8600 “ 

95 — • - 


B600 

96 


0603 

97 


6603 

98 


■ B60Q 

99 


B600 

IGQ 


B600 " 

101 


8600 

12 2 


S600 

lu3 


860 3 

104 


■ B600 

105 


BGOO 

106 


B60D ■■ 

10 7 ■ 


B60O 

108 


6600 ■ 

109 


0600 

110 




pf^OGRAM 


060 0 


COC 6600 FTN V3.0-32A OPI=2 



U5._ 


17 


120 


16 

93 


68 


.125 


130 


101 

102 




135 


140 


14S 


150“ 


15J 


16D 


.LeN=0 t NEXT=IFIRST=i 

1(1=2401200 5 3560 00 00 00 000 
CALL OBIJF(O) 

CALL OCUFtll 

IF(NCON*2000)16tl6,17 , 

PRINT S29,NCON 

STOP . 

rPHO=IPHO*-l 
rF(IPHO-44»47,93f93 
NROM=HE=IPHO-43 

GROUP ALL NEW ODJECTS A'O 
DO 63 I=HE,fJROH 

jOYCE=Z15H-43»ni-l( 

■ N£IL=A< JOYCEtPU 
L=l' 

DO ICO K=1»MEIL 
DO 101 NU=liLINE 

IFTA( JOYCE »K»-G<6*NU-5) )131 flOZ.lfli 

CONTINUE . . 

H0ftT = NU-86» t UtU'-l J/86) 
M=(lMNU-l)/'86»’'9 - 8 , „ 

J = M 

106 iF(A<J,3»MORT-2)li08,lG8,236 
296 C(l,L»=A(J,3»KORT-2l 

C(2fU = A{ J,3»M0RT-1) 

C(3tL)=A{ J,3*M0RTT - 

_ J=JH . 

L=LFl 

IFU-M-0) 106il06»108 _ 

108 AtM,3»(10RT-2)=-l. 

_10Q CONTINUE 

‘ NA=3» (L-l> 

CAHU,12)=NA 

IF(NA.t0.ai GO TO 83 
CALL WRITHStll,C,NA,KOISKI 
'■ kOISK=KOISK + (NA+4 07 1/40 6 

83 CONTINUE _ 

GO TO (.47;95» ,IGO 

READ IN REFINED IMAGE COORDINATES 
"47 JOYCE=215HPHO-43*( (tPHO-11 /43I 
H£IL=0 
CALL OSUF(l) 

DECODE(80i525.BUFFtIXBUF))3(l»50l 
IF {Bt2t50» .ea.i.O I NCOH=I>H0 ‘ 
IF(Q13.50)“1. »318,317<318 
I=IPHO-NCOHfl 
IFa-20 1352.352,351 
PRINT 572 

STOP 

352' IFU.6T .LOS) LOS=I 
318 IFtLARRY)91,91,90 

90 FL=812,51J 
CAHtIPHO,21=FL 

91 IrtAGE=l 

71, r.Ai I nnuF n » 


THEIR IMAGES "OR PLATE AMD PUT IN 


FOR PLATE BEING RESECTED 


a<2 «5Q ) t U(3,50> t6(l t5l> «8 C2| 


'317 

351 


04/C2/75 20 

.40.54. PAGE 3 

a&oo 

111 

8600 

112 

0600 

113 

■ ‘ B600 

114 

B60(l 

115 

B600 

116 

B60G 

117 

B600 

118 

B60u 

119 

‘ B6D0 

120 ■ 

SK B6G0 

121 

B60(j 

122 ■ ■ 

B60D 

123 

B6o0 

124 

S6GC 

125 

B600 

126 

□ 6CC 

127 

B600 


D60G 

129 

B60Q " 

130 

B600 

131 

O60C ■" " 

132 ■ ■■■ 

8600S 

5 

■ B600 

138 

8600 

139 

B680" 

140 ■ 

B600 

141 

B6(!0 ’ 

142 

B600 

143 

B60D 

144 

B6O0 

145 

" ueoj ‘ " 

146 

660 0 

147 

0600 

148 

O600S 

6 

B600 ■ ■ ■ 

150 

6600 

191 

— B600 ' 

■152 ■ 

6600 

153 

B600 ■ 

154 

B6CC 

155 

B60Q 

156 

11 B6GC 

157 

8600 

158 

6603 

159 

6600 * "■ 

160 

86GO 

161 

8603 

162 

B&CO 

163 

!360C‘" 

164 

8600 

165 

0600 

166 ■ " ■ ■■ 

B60C 

167 

B6DQ 

168 ■" 

U600 

169 



. .PROGRAM B600 , 


CDC 6630 FTN i/StO-Sai* OPT=2 04/02/75 20.40.54. 


PAGE 


0© 


_gECOOE< 80 .52 6.0UFF(I XaU FH ( 0 ( J. I MAGE) < J=l< 3»_*._IJi?L 


170_ 


IFUTESn 1 . 8 , 18,21 
18 IMAGE =IHAGE *-1 

IFaHACE- 30 ) 36 , 36,20 
20 IS= 011 , 50 > 

PRINT 545,13 

STOP „ 

SET UP ORD^R LIST 


175 


_ioa_ 


185 




150 



205 


210 


215 


220 


LIST FOR PHOTO 
21 LAST=0 

IF«Dt2,5fl>>39,39,38 

38 NR0W=4 

GO TO 42 

39 IF(OU,50>)40,4Q,41 

41 KROW=l . . _ 

NHI=6 

_ last=i _ . ... 

GO TO 43 

40 HROH=l 

42 NHI=IHAGE 

43 00 23 K=NROW,NHI 

J=01i,K>».D00001 
FK=J»lfl00300 _ _ 

FK=B(1,K)-FK 

NU=i 

32 MORr = MU-86*((HU”l)/86) 

H= (It tNU-i)/66l *9 - 8 
IF lAlH, 3* rtORT -21) 25,24,25 

24 AtM,3*-HORr-21=BTl,K) 
A(M,3»MORT-l) =B(2,K1 

A(M,3*-H0RI1 _jB(3.K) 

LINE=UIH£H 
I=&*-MU-5 

IF{Gtl)-FKI26,49.26 
26 G(12C01)=C(I1 
C(12002l=G<Itil 

Ga2003) = C(It21 

Ca20i)4)=G(I»31 
G(120O5l=G(I*4) 
G(12006)=G(I+5> 

NU!=HU+i 

00 28 L=NU,NCON 

IF(Glb»L-5) -FK) 28^29,26 

"28‘ CONTINUE 
IS=FK 

■■ PRINT 526 , IS 
STOP 

25 GUI =G(6"U-5t 

G(I*-1> = G(6»L-4I 

G(I«^2> = G(6»L-3) 
G(lF3)=G(6»L-2) 

G{IF4> = G(6*L“1) 

GaF5> = G(6»LJ 
Gie^L •’5 1=0( 12 001) 

G(6*L-4 1 = 0( 1200 21 


- STO^H IMAGES BT OBJECT - SET UP OBJECT FILE 
SAVE S53RDINATES OF POTENTIAL RESECTION OBJECT! 


Q600 

170 

■ C600 ‘ 

1 

C60Q 

2 

C600 ■ 

3 

C603 

4 

0600 ■ 

5 

C600 

6 

' 0600 

'178 ■ " 

8600 

179 

B600 

180 " “ 

8600 

181 

0600 

182 

0600 

163 

“ 8600 

164 

0600 

185 

B600 

186 

8600 

18 7 

B600 ' 

188 ’ ■■ 

6600 

189 

B600 

150 • " 

0600 

151 

B600 ■■ ■ ■ 

152 " ‘ 

8600 

153 • 

B600 ■■ ■ 

154 

B60il 

155 

36 G V 

156 

8600 

157 

8600 

198 

B600 

159 

8600 

2C0 " ' 

8600 

2C1 

8600 

202 

8600 

203 

8600 

204 

8600 

205 

B60U 

2C6 

0603 

207 

8600' ” 

208 ■ ~ 

8600 

205 

8603 

210 

S6QQ 

211 

B600 ■ ■ ■ 

212 ' ~ 

8600 

213 

B600 

214 

B&CQ 

215 

8600 " 

£16 •- 

0600 

217 

0600 

218 ' 

8600 

219 

B6CG 

220 

BGOO 

221 

B60D 

222 

8600 

223 

0600 

224 

0600 

225 
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G{6'fL-3) = G{120C3l 

G(B*L-21=Gtl20Q4) 

• G(6»L-il=G(i2C«05l _ 

''G(6»L) =G(12006»' 

225 49 NEIL = NEILH . 

A(JOYCE,McIL)=GCI> 

If <K-9) 50 ,50,23 

^(1,K^6C)=G(I♦•11 
nt2,Kt'60) = 6tIi'2) 
a(3,Kv6fl) = G(If31 
GO TO 23 

IFTFK-G (6»NLr-5l ) 30, 31,30 

HU=NU«-1 

GO TO 32 “ 

L = H + 1 

IF{A<L,3*M0RT-2) >33,34,33 
L=Cfl 

IF(L-H-6) 35,35,37 

IS=Gt6»NU-5) 

PRINT 5 37, IS 
STOP 

A (L,3»MORT-2)=0tl,KI 
A (L,3«M0RT-1) -B12.K) 
A<L,3»MORT) =BC3,K) 

I=6»NU-5 

CO TO 49 
23 CONTINUE 

IFILAST146,46,44 

_ 44 IFtIMAG£-9)46,46,45 

'250' " ' 45 NROH=lO 

____________ LAST=C 

“go to 4f 

_ 46 AlJOYCE,31)=NEIL 

C 

_2S5 C PHOTO RESECTION PHASE 

c‘" 

C INITIAL APPROXimnONS 3 

■ “3'i5"n=0 


50 

2 30 



30 

235. _ 31 

35 

... 33 

37 

240 

34 


245 



CAH^A_PARAHEr£RS 


DO 51 J=l,2 

260“ C<J,1)=B<J,66> 

_ C(J,2)=0. 

C<J,3)=1. 

51 CONTINUE 

^ ' C" “ ADJUST APPROXIMATE AZIHJTH PARAMETERS FOR SWING AND 

265 _ 0(1, 1)=Q(2, 41-0(2, 5) 

0(1, 21 = 8(3, 4)-Bt3, 5) 

0(2,11=8(1,641-0(1,651 

0(2,21=8(2,641-0(2,651 

0<2,4)=Dt2,l)*D(2,l)f0(2,21*0(2.2> 

270 '0(1,31= (0(2,11*0(1,11 *0(2, 2 l»Oti, 21 I/O (2,4 1 

0(2, 31 = (0(2, 21*0(1, 11 -0(2, 11 *0(1, 21 1/0 (2, 41 
SCAUE=SORT (0(1,3) *0 (1 ,J1 O (2,31 *0(2,3 1 1 
C(3, 21 = 0(2, 31/SCALE 
C(3,11=FL/SCALE 

275 C(J, 31=0 (1,31/SCALE 


06CO 

226 

86Ct 

227 

8600 

22 0 

“'060 0 

229 

8600 

230 

8600 

231 

8600 

232 

8600 

233 

0600 

234 

8600 

235 

8600 

236 

B6U0 

237 

BbdO 

230 

“ 9600 

239 

860 0 

240 

8600 

241 

0600 

242 

060i) 

243 

6600 

244 

B60u 

245 

0600 

246 

■ 8600 

247 

8600 

240 

' 8600 

243 

6600 

250 

■ 8600 ■■ 

251 

8600 

252 

8600 

253 

8600 

254 

860 0 

255 

8600 

256 

" 8600 '■ 

257 ■■ 

86C0 

250 

S6CC 

259 

8600 

260 

060 3 

261 

6600 

262 

060 0 “ 

267 

060 0 

266 

3600 

269 ■ ' 

0600 

270 

8603 

271 ■ 

060 0 

272 

■ 0600 

273 

8600 

274 

B6CQ 

2/5 

8600 

276 

8600 

277 

0600 

270 

8600 ' " 

279 

0600 

280 

OoCO 

281 

BbOG 

£82 

J 

£63 " 

8600 

264 
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■C& 




52 

_ ..280 .. 


53 




285 

— 

54 

92 

290 

■- 





295 






3 00 







305 






310 





55 


C 

56 

315 




c 

57 

320 


58 

- “ 

— 

59 

6 G 

325 



- - - ' 


61 


.ORIENTATIQH.FACJORS_IN 3. ARRAY. 


IS=0(1. 501 
PRINT 535. IS 

PRINT 5 57, as, <C (J.n . J=l. 3 » .I=i*2I 
IRE=2 . 


62 


330 


K=G 

C<1.Kt!, 1=C(2,K+31 »C(3,l<f3l 
C tl,K*5>=-C(2.Kt3)*C(3,KT2) 

C(1,K + 6)=C<2,KT'21 
Ca,K»iei =-C<2,KfZ) *C l3,<i-3 J 
C(i.KUl)=CI2.K<-21»C<3,Kai 
C(l,KH21=C<2,Kf3J 
C(2,KU0) =Ca,K«-A) *C(l,!<t-2l 
C(2,Ktll>=C(l,K*^5)*C(i,KT2) 

C{2,K+i2)=C{l,K+2l »C(2,Kai 
C(3 ,Ktl0)=-C(l,K*’4) *-G [l,<»-3) 

C(3,KtLl> =-C(l.K*^5l »CC1,<^3 ) 

C{3,Ktl2J = -Ca,K»3l*C(2t<i-2> 

C t^,K^4) = G a ) K^3) »C 13, K*-2t + C t 2 , K«-12> *C ( 3 ,< +3 1 

C (2,K*-5)=C (1, Kf3 I '•C (3 ,K*3> - C I 2, Kf t2) *C ( 3 *-2 ) 

C(2,K<-6) = -Cll ,K^2I »C(2,Rf3 I 

CO.K+A l = C (1,K + 2I *C(3 ,Kf2> 1, K+'IO ) *C ( 1 ,<< »3 1 

C(3,Kt5) = C(l,Kt2)*‘CC3,K*3l «• C{ 1 , K*-lt I *C ( I , < *3 I 

0(3,KT5)-C(liKT3l*C(2,K+3» . 

00 55 1=7, M 

«=KH .... . 

C 11,M) = 0. 

C{2,Hl=-Ct3,M-3) 

c(3,rt)=c i2,M-3i ■ 

C ( I-5.K+13 I =C ( I-6,K»5I , . . , 

C(I-6,K+14>=-C(I-6,K»4| 

Ca-6,Ktl5» = l). 

CONTINUE 

GO 10(56,71) ,IGO 

' 'CLEAR NORMAL EQUATION' 0 ARRAY TO ZERO 
DO 57 1=1,6 . ...... 

00 57 J=I ,7 
0(1, J) = 0. 

CONTINUE 

COMPUTE P TERMS FOR PASS POINTS USED FOR RESECTION 
iF(B(2, 50)159,59, '58 

H=4 . . . 

NR0W=8 

GO TO 62 

IF(B(3,50) 161,61,63 

N=1 

NROH=5 
GO TO 62 
H=1 
NROW=9 

00 63 NU=H,NROM 
00 64 K=l,3 


! fl,*/02775 

HO,40,54« PA6E 6 


6600 

265 


“'6601) 

286 


8600 

287 


0600 

288 . ■ 


B600 

289 


8600 

290 


B60C 

291 


0600 

292 


B600 

293 


8600 

294 


B6G0 

295 


0600 

296 


B6GD 

297 


B600 

298 


B60Q 

299 


B600 

SCO 


0600 

301 


B60C 

302 


8600 

303 


8600 

304 


B6DO 

305 


36li0 

306 


B60G 

30 7 


B60O 

303 


BdCC 

309 


■ ■' 06OG 

310 


8600 

311 


060u 

312 


8600 

313 


B600 

314 


66(>Q 

jlS 


d&o u 

Slo 


8600 

317 


■■■ B600 

316 


B6D0 

319 


B60O 

■ 320 


B6Q,j 

321 


0600 

322 


B6C0 

323 


860J 

324 - • ■ 


U6Q0 

32 5 


8606 

326 


6600 

327 


' B60C " 

■ 32 8 


B60O 

329 


8600 

330 


8600 

331 


8600 

332 “ 


B6if0 

333 


■ 9600 

334 ■ ■■ 


B60u 

335 


8600 " 

336 “ ■ 


B6C0 

337 


O60C 

338 


86CC 

339 
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. C { K , 16 ) =8 t K , N,U*-&0.!_;C ( K , 1), 


335. 



64 COWriNUE 

'”dO 65 1=17*20 

DO 65 1=1,3 . , 

Ca*LI=C(l,K)'"Ctl*16) ♦ Cl2*K)»C(Z*t6) + 3 (3 ,K)*C(3,163 
. K=KM. 


340 


65 CONTINUE 
00 66 1=1,2 
OO 67 L=l,4 

P(I,L) = <OtI«-l,NU)*C<3,L + l6) +FL *C <1 , L+ 15 ) I /C (3, 17» 

67 CONTINUE 
00 66 L=5,7 


345 


350 


355 




360 


■ptI,L)= (-B(Ifl,NU»»C(L-4,5 t 
66 CONTINUE 

p{i,8i = -paTi) 

66 CONTINUE 

'C CONTRIBUTION TO NORMAL EQUATIONS 

00 63 1=1,6 

DO 63 J=I,7 
DO 63 K=i,2 

0 (I,J»=0(I,JI+P<K,ItlJ*P(<, Jtl) 

63 CONTINUE 

C FORHARO solution 

00 69 1=1.6 

soR = i./sqrt"io'<i,ij 1 

00 72 J=I,7 _ 

0(I,J)=0(I,J>*SQR 

72 CONTINUE 

IFtI-6)73,74,74 
73 IP1=IF1 


FL“*C(L-4,'Iv' 3> i*C(3,i)/Ct3,17r 


365 


370 


375 


380 


385 


00 65 L=tPl,6 
00 69 J=L,7 

D(L.J> = 0(L,JI-Da,L)»0tI,J> 

_69 CONTINUE 

■ OACK SOLUTION 

_74 0(6,7)=D(6,r»/0(6_,6) 

"■ 00 75 1=1,5 

NHI=6-I . . , 

NHIP1=NHI«-1 

00 76 J=NHIP1,6 . . _ . . 

0 CNHI ,7 > = OtNHI,7l-0 tj,7)*3 ( NHI, J3 

76 CONTINUE 

0 tNMI,7l=D (NHI,7I /O (NMI ,N>in 

75 CONTINUE . . 

00 77 1=4,6 
DCI,7)=OtI,7)*CU,ll 

77 CONTINUE 

ADD LEAST SQUARES RESULTS_ TO_CftHERA PARA Hi TERS _IN__C ARRAV__ 
“'■"do 76 J = i,3 

ClJ,l)=C(J,lM-DtJ*3,7» 

C(J,4»=CHJ,7) 

CU,5) = SQRT(i.-C(J,4)'fC(J,4)) 

C (J,6J=C( J,2»*C (J,5 l+C( J,3) *CtJ,4) 
C(J,7}=C(J,3)*CIJ,5)-CtJ,2)*C(J,4) 


8600 

B603 

B6C0 

B6CU 

660G 

B6O0 

8600 

' B600 ■ 
B6U0 
B600 
0600 
6600 
0600 
B60l 
6600 
6600 
0600 
R60I 
6600 
6600 
6600 
B6C9 
8600 
BdOO 

060(i 

■ 360d ' 
6600 
B6C0 
8600 
6600 
8600 

■ 8600 
6600 
G60o 
6600 
SpOO 
8603 
0600 
B600 
0600 
B600 
0600 
0600 
B6CG 
8600 
BoOC 
8600 
B&bO 
8600 
8600 
8600 
B60J 
8600 
860u 
B6C0 


340 

341 

342 

343 

344 

345 

346 
'347 

346 

349 

350 

351 

352 
"353 

35 4 

355 

356 

357 

358 

359 

360 

361 ■ 

362 

363 

364 

365 

366 

367 ■ 

36 6 

369 

370 
371" 

372 

373 

374 

375 

376 

377 
376 

37 9 

380 

381 

382 

383 

384 

385 

386 
337 

38 6 

389 

390 

391 

392 

393 

394 
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U Srf WU w 


M • V ” 


tj3 

-S) 





\ 

76 1 


..390 

C 

1 

79 1 



C 

70 1 


3 95 


i 

( 




1 

) 


4Q0 


60 < 




61 i 
1 

1 


405 


82 ; 



c 

95 : 


410 


85 1 



c 

6^ 1 



c 

c 

i 


415 


271 1 

1 




86 1 


420 

c 

c 

1 


* 


301 
234 1 




1 

105 1 

— 

425 

— 

1 

297 


430 


87 
S3 1 




89 1 




96 

- 

435 

, . 

94 

104 


440, 




OF CORRtiriONS FOR ORIENTATION PARAHtTERS 


C( J, 2 t=ct J,6I 

C<J,3I=C(J,7) 

CONTINUE 

TEST MAGNITUDE 

00 T9 1=1,3 , 

IF (AiJS( 0(1,7) >-.0 00 01 179, 7 9, 52 

CONTINUE _ _ „ _ „ _ 

STORE CAMERA PARAMETERS AS (SOMPUTED FROM PHOTO RESECTION 
CAH(IPHO,l)=Ba,50) 

00 90 J=l,3 

CAMCIPHO, J18)=C( J,3> 

CAHCIPHO, Jt5)=C(J,2> 

CAM(IPHO,Ji-2)=C (J,ll_ 

CONTIHUE 

IF(Qtl,51) 116,16,81 
NR0N=IPMQ 
HC=M5H 
IGO = 2 

LOS=2*UOS+3_____ 

IF(IRE-2) ',6,82,82 
STOP 

CLOCK SIZE 

IF(LINE-197a) S',, Bit, 85 
PRINT 5AA,LINE 

STOP 

PRINf ‘’543‘,IPHO,LINE 

BLOCK ADJUSTMENT PHASE 

00 86 1=1,259 ■“ 

00 65 J=I,259 

8CI,J)=C, 

CONTINUE 

RESET PLATE WITH NEW OBJECTS AND FORM C .ARRAY ROWS 1-15 FOR ALL 
PLATES ON WHICH THE NEW 03JECTS APPEAR _ 

INSECTsO ■■ 

K0ISK=1 

ME=1 

NA=CAH(HE,12) „ 

IFINA.EQ.O) GO TO 297 
CALL R£ADHS(11,U,NA,K0ISK) 

KOtSK=KOI3KKNAt407 1/40 0 

IF(HE-11 87,87,96 

IF(IPH0-43) 88,88,89 

UASr=IPH0 _ _ _ 

GO TO 94 ■ ■ ■ -- - - ^ . 

LASr=43 
GO TO 94 
NHI=LAST=HET42 
IF(NHI-IPH01104,lCI4,”l23‘ 

NMI = 1 

00 71 L=NHI,LAST 
K=20»<L-1J - 86G‘*((L-l»/43) 

DO 97 J=l,3 ' 

Cf J,K*1 ) = CA«(L, Jt2) 


B6B0 

395 

8600 

396 

B6 Ou 

397 

660J 

398 ' 

B600 

399 

B6i)u 

40 0 " 

BS'OO 

401 

B&GO 

4(i2 

B600 

40 3 

06bO 

40 4 

BoOtt 

40 5 

B600 

406 

B603 

40 7 

0600 " 

408 

B600 

409 

B60C 

410 

0600 

411 

B60U 

412 ' 

8600 

4l3 

S600 

‘414‘ 

8600 

415 

8600 ■ ~ 

416 ~ 

8600 

417 

B600 

418 ” 

B600 

419 

0600 

420 

8600 

421 

8600 

422 "■ 

B600 

423 

B600 

424 ” 

0600 

425 

B6fld 

426'“ 

B600 

427 

8600 

432 ■ ■ 

8600 

433 

B&diJ ' 

434 

B6G0 

435 

B600 

436 

0603 

437 

B600 

.438 

B6 00S 

7 

8600 

440 

B600 

441 

8600 

442 ~ 

B600 

443 

0600 

444 

0606 

445 

8600 ■■ 

446 

B500 

447 

8600 

448 ~ 

0600 

44 9 

B6DG 

450 

B600 

451 

0600 

452 

B6C0 

453 



Ci> 

o 



PROGRAM 
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i 



1 

. 97 1 

445 

— 

- 


450 






455 

C 

71 

123 

303 




460 


110 




465 


111 

115 



107 

' 


109 

470 


125 




475 



146 


480 


114 

126 


c 

112 

485 





116 

490 



495 


117 


C ) = CAM(t, Jf 51 

C<J,K»3 ) = CAH<L» J^6) 


If < INSECT .€Q. 0> GO TO 92 
Ca,K+A)sC(2,Kf3) »C(3,K+3t 
C(lfK+5) = -Ct2 tK^SI-'CCStK*-’ ( 
Ca,Kt6l = C(Zf K«-2) 
C <^> K ^‘*>= C < 1 » K ^■ 3 » * Cl 3 tKf 2 l 
C(2,K*5) = Ctl,K<-3l *Cn,K<-3l 
C<2,Kt-6 ) = -C(l ,K^2I •C(2,K«-3I 
C (3.K+A l = C( 1, K«-2l •C(3tK»2l 
C(3,Kf5> = C(l,K*-2l »Ct3,K»3l 
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20,40.54 


8600 

454 


~ ' ■ 0600 

455 


0600 

456 

■ ■ ■ 

■■ ' B6Q3 

457 


8600 

45 8 


B600 

’ 459 


8600 

460 

rcTlTK'+ 2 } * C(' 2 , Kv“2?*CU7K^3) 

B600 

461 ■ 

- C<l,K + 2) *0(2, K»21 »C<3,K<-2) 

B600 

462 


B6Q3 

46 3 

- C<2,K + 2)»C(3, Kf3»»C«1,KF3» 

B6C0 

464 

♦ C<2,K + 2) *C(3, Kt 2) »C (1,K+31 

" ■ 3600 

465 


B60C 

466 


PAGE 


71 CONTINUE 

Ntw OUJECTS being PROCESSiO ON PLATE 


IFtINSECT.EQ.O. AMO. NA. £0.53 CO TO 169 

IfdNSECT.EQ.il GO^JO 36p 

liO J = B(1,NI».D 00 0IH 
FK=J 

00 111 L=1,IPH0 


IF(FK-CAH(L,1I) 111.115,111 
CONTINUE 

K=20'‘(L-il -_860»( CL-^/!t3_l_ 
£F(LARRT1 109, KToill)? ■ 


FK=y!i,Nt-FK 
IF<H-1) 125,125,126 


00 1A6 1=1,3 
E d,5fl) =0 . 

DO 1L6 J=I,3 
£d,Jl = 0. 

CONTINUE 
HORT=0 

00 114 NU=r,LIN£ 
IF(SQR-GC6*NU-5) > 114,126,114 


COMPUTE C ARftAY COLS lo THROUGH 20 FOR IHAGE 


D(1,HORT)=NCON=L 
00 115 1=1,3 
M = 6»NU-5*-I 

Cd,K+l6l=G(HI-CtI,K+U 


J=K«-4 

DO 117 L=17i'Z0" 
K=K«-L 

00 117 1=1,3 


C(I,KI = C(i,J>»C(l,K*-16l 
J=J»i 


* : (2,J1 »C (2,K+l6l *■ C(3,J»*C«3,Kd6) 


0600 
B6G0 
0600S 
0600 
8600 
0603 
O&QO 
8600 
8600 ■ 
8600 
B600 
B 6 G 3 
8600 “ 
06CC 
B609 
8600 
0600 
86 G 0 
■B6C0 ■ 
B6Q0 
B6GJ 
8600 
8600 
D 6 C 0 
B6C0 
0600 
B60Q 
06QC 
6600 
S60.I 
8600 ■ 
8601) 
B6C0 
D&OO 
B60G 
B60I) 
B600 ■ 
B600 
B6G0 
B60 9 
B600 
6603 


46 7 

46 8 
8 

473 

474 

47 5 

476 ■ 

477 

478 

479 
490 
481 
432 '■ 

483 

484 

48 5 
,.86 
467 

488 ■ 

489 

490 

491 
492" 
493 
49-, 

495 

496 '■ 

497 

498 ' 

499 
5£!C“ 
501 
SC 2 
50 3 
50 4 
505 
50 6 
507 
50 8 
50 9 

510 

511 
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PAGE 


10 


-§i 




500 , 


119 


COHPLtTE P COEFFICIENTS DF OBSERVATION ejUAJIONS FOR IMACE 
Vo 118 I=lt2 

00 119 U=ii4 . _ . .... 

H=KH. 

P(I,U = (Q (H-1 ,N>*C( 3.H + 15I t-FL’C (r,H«-l6) >/: (3 .R*!?) 

CONTINUE 

H=KH ■ 

DO 118 L=5,7 



505 


pii,u = (-tj(iniN)*ca 
'llfi CONTINUE 

ARRANGE AUGMENTED COEFFIIIENT MATRIX 
■ 00 120 1=3,4 

00 121 L=l,3. 


4,<^5 »-FL*CtL-4,Mt3» J* CAH(1.5>/Ct3,K*17> 


510, 


515 


_52 0. 


525 


121 

122 

'120 


129 
135 
137 

130 
,.131 

,132 

133 


STATION OBSERV. EQUATIONS FOR TARGET QUALITY 


530 


535 


.134 

..136. 

136 

139 


140 

141 


540 


142 

143 


545 


550 


144 


P(I,L>=-P[I-2,L+4) 

CONTINUE 
00 122 L=4,9 
P(I,LI=P(r-2,L-2) 

CONTINUE 

PI 1,10) =-P!I-2,lJ 

CONTINUE 

WEIGHTING CONTROL 
H=6»NU-1 

IF(G(HJ 1130,129.130 
IFIGIH*!) ) 135,136,135 

IF<JERRY1131,131,13 7 

WEIGHT=G(H + 1) 

GO TO 133 

IFUERRY) 131,131,132 

WEIGHT = HTCON _ „ 

GO TO 133 “ 

H£IGHT = G[M) 

DO 134 1=3,4 
00 134 J=l,10 
Pa,J)=P(I,J)»HEIGHT 
CONTINUE 

WEIGHTING EQUATIONS FOR LOCATION 

RA0IUS=SQRnBt2,N)»»2 3 _ . 

'iFtRADIUS-. 07)143, 143 ,’138 

IFIRAOIUS-. 12) 139 ,139 ,140 . 

HEIGHT=WT7i2 

GO TO 141 . . . , 

HEIGHT=WT12P 

00 142 1=3,4 

"OO 142 J=l,lfl" 

P<I.J1=P(I,J)*HEIGHT __ 

CONTINUE 

COMPUTE NORMAL EQUATIONS- FOR IMAGE AND STORE IN E AND A ARRAYS 
00 144 1=5,13 

00 144 J=l,10 

PtI,J>=0. 

CONTINUE 
00 145 1=1,9 
DO 145 IP1=3,4 
DO 145 M=I,10 

p:if4,M)=p<i+4,H) ♦p(iPt,ri opaPitHi 


OF IMAGE ON PLATE IRESOLUTION) 


B6CQ 

512 

~ 9608 

513 

8603 

514 

8600 ’ 

5l5 

B6C3 

516 

B608 

517 

B600 

518 

9600 

519 

B600 

520 

■ B603 

521 

B60O 

522 

8600 ■ 

523 ' 

B60J 

524 

■' 9600 

525 

B6QQ 

526 

0600 

527 

B6G0 

528 

B&CO 

529 

9600 

530 

B6Cu 

531 

ceuB 

7 

D60C ” 

533 

06CD 

534 

B600 

535 

060Q 

536 

B600 ~ ■ - 

537 ■ 

B60D 

538 

B60u 

539 

B603 

54C 

Q&OO 

541 

8600 

542 

■ ■ 6600 

543 ' 

B6G0 

544 

9600 

545 

9600 

546 

8600 ' 

547 

0600. 

548 

B600 

549 

0600 

550 

B600 

531 

0603 

552 

B60 0 

553 

0600 

554 

■ ‘ 8600 

555 " 

8600 

556 

8600 

557 

S60C 

556 

B60l) 

559 

B6QC 

560 

“ 9600 

561 ■ 

9600 

562 

B600 

563 

9600 

5 64 

0600 

565 

9600 

566 
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555 


560 


CONTINUE 

00 147 I=ii3 
00 147 J=Ii3 
£U,JJ = E(IiJ1+P(H-4,J) 
147 CONTINUE 

E <1 ,581 =E<1,58) *P<5, 101 

£<2,581=E(2,S8)+P(6,10I_ 

£t3,58t=E(3,S8l+P17,10) 

DO 148 1=1,3 

N£IL = 6’‘HORT-3 

DO 148 J=4,9 

N£H. = HSILtl 


5 65 


570 


• 




E(r,NEIL>=P(I+ 

4., J)_ 


148 

CONTINUE 




IF(IPHO-431ZOO 

,20Q, 

201 

201 

IF(ME+42-IPH0) 

202,203,203 

2C2 

I=6»NCON-6-6*( 

HE-ll 



GO TO 204 



.20 3 

I = 6*NC0N-6-6» ( 

IPHO.- 

431 


2C0 

204 


575 

1 


580 

149 ( 

585 

151 1 

'c 

113 


153 

154 


590 

155 


■ 152 
157 

595 

159 

156 

C 

C 

15_6_ 

600 


605 

'128 


GO TO 204 
I=6*NCON-6 
DO 149 K=8,13 
N=K-4 
J = I 

r=in 


AII,LLL1 = A(I,UL) l-PCKtlD) 
00 149 L=H,9 
J = JH 

A(I,JI = Aa,JI+P(KrLl 
CONTINUE 

N=NH 


IF(N-NA/3I 110,110,151 
KC=1 

APPLY POSITION HEIGHTS 3- CONTROL STATIONS 

IFtCt6»NU-llll52,lS2,153 _ 

IFtJERRY) 154,154,155 ’ 

£(1,1)=E(1,11*G(6»NU-1I 

E(2,21=e (2,2 l*G<5*NU-ir 
GO TO 152 

E(1,1)=£(1,11*HTC0N 
S(2,21=E(2,21*HTCON 
IF(G(6»NU1 1156, 156,157 

IF(JERRY) 159,159,158 

E(3,3) = E(3,31»CC6*tlUl 
GO TO 156 

E(3,3I=£(3,31»HTC0N 

FORrfARO SOLUTION OF NOR-IAL EQUATIONS FOR NEW OBJECT 
FORHARO SOLUTION OF OBJECT ROWS IN £ ARRAT 

NEIL=6*«ORT+3 _ 

DO 127 1=1,3 ' 

SQR=1./SQRT(F(I,I>) 

00 128 J=I,MErL 

E(I,J1=E(I,JI*SQR 

CONTINUE 

E(I,5B1=E(I,56I*SQR 


04/02/76 

28.40 . 

8600 

567 

66Ca 

568 

B600 

569 

8600 

570 

B600 

571 

0600 

572 

B&OO 

573 

B60O 

” 574 

D50Q 

575 

8600 

576 

B6DO 

577 

B600 ■ 

578 

B6GQ 

579 

■ B600 ■ ' 

' 580 

0600 

581 

B60G 

582 

6600 

583 

0603 

58-, 

0600 

585 

'■ B600 “ ' 

■ 586 

B6QO 

567 

■■ 860Q 

'■ 58ft 

B60G 

589 

B60C 

690 

060J 

591 

“ B600 

* 592 

B60C 

593 

B600 

594 

8600 

535 

' B50G 

596 

8600 

597 

■ 8600 

59ft 

B6G0 

599 

8603 

630 

860U 

601 

8603 

602 

8600 

603 

"■ 3600 “ 

■■■ 694 

8608 

6C5 

8603 ™ 

636 

6600 

6C7 

Q60C 

60 ft 

8600 

60 9 

“ B60C 

■ ' 610 

8600 

611 

sees 

612 

0600 

613 

8600 

614 

8600 

615 

** 8600 

■ 616 

B6QG 

617 

86C0 

61ft 

8600 

619 

B600 

620 

8600 

621 


PAGE 
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# 










. , „ . 610. _ 


615 


620 


625 


630 


635 




645 


650 


655 


C 



IFtI-3) 160,161,161 

B60C 

622 

160 

IPl = Ifl 

- ■ 3600 ■ ' 

623 


00 127 H=IP1,3 

8603 

624 


DO 162 J=M,NEIL 

0600 

625 


E(M,J)-EtM,J)-£CI,H)»E(I,J) 

0600 

626 

162 

CONTINUE 

0600 

627 


E (H,58) =E(«,S6) -£ (I ,H) »£([, Sd) 

0600 

628 

127 

CONTINUE 

3600 

629 


EFFECT OF OBJECT RONS O') PERTINENT CAMERA ROWS 

0600 

630 

161 

DO 163 J=1,H0RT 

0600 

631 


IPl =0(1, J) 

B600 

632 


IF (IPHO -43) 20 5,20 5,20 6 

B6C0 

633 

206 

IF (HE*-*, 2- IPHO) 207 ,2 08,206 

B60u 

634 

207 

LASr=6*IPl -6*(HE-1) 

Q6C0 

635 


GO TO 209 

0600 

636 

208 

LAST=6»IPl -&»«IPHO-43) 

B600 

637 


GO TO 209 

B600 

638 

205 

LAST=6» IPl 

" ■ “ 8600 

639 

209 

H=LAST-5 

860 0 

640 


00 163 NROW=i,3 

3600 ' “ 

641 


NEIL=5* J-3 

3630 

642 


DO 163 I=H,LAST 

■ ' B60u 

643 


164 

TVs' 

■'■ 167 

211 
^ 212 . 

213 

210 

214 


166 

163 


HHI^J _ 

NCON=MEIL 

NEH,= NEILt-l 

A (I,LLL ) = A(I7i.uC) “E(NR6w.MiILJ»E(NROH758> 

00 164 JOYCE=I.LAST 
NCOH=NCONfl 

A a, JOYC£l = A( I» JOYCE) -EtNRDH,NeiL)*EtNROW, NOON! 
CONTINUE 

I F < J- HO RT ) 1 65 , 1 63_._1 b3 . 

NMI=NMH-1 ‘ ' 

IF(NMI-MORT 1167, 167, 153 
IS*D(l,NHn 

IF(IPHO-43»210,210,211 

IFlMEt42-iPH0)212,213,213~' ■ 

N0=6»IS-6»(HE-1) 

GO TO 214 ■■ 

NU=6»IS-6" (IPHO-43) _ „ „ 

GO TO 214 
NU=6'‘IS 

IHAGEsNU-5 " 

00 156 JOYC£=IHAG£,NU 

HCON=NCONn ■ “ " 

A (I ,JOYC£) = A( 1, JOYCE) -E (NRO H , N£I L) * E( NROH, NOON) 

CONTINUE 

GO TO 165 

CONTINUE 

IFIKC)12S,125,169 

FORHARO SOLUTION ‘0F“C‘AM:R’A ROHS'TOR Th£"RESET“ PL’ATE 


3600 

644 

B600 

645 

9600 

646 

■ 0600 

6*7 

0600 

648 

3600 

6*9 

3600 

650 

B60G 

651 

3600 

652 

B6CC 

633 

B600 

654 

B600 

655 

86S0 

656 

3600 

- 6S7 

B600 

658 

~ B60O 

659 

8600 

660 

0633 

‘ ■ 661 

3600 

662 

86CG 

663 

3600 

664 

B60'0 

■ 665 

B600 

666 

' 8600 

" 667 

3600 

668 

0600 

669 

0600 

670 

" 8600 

671 


169 

K = 0 

B6G0 

672 


IFTIPHO-43) 170,170,171 ■ 

0600 

673 

170 

N=6»ME 

B60D 

674 


GO TO 174 

9600 

675 

171 

If (HE +4 2- IP HO >172,173 ,173 

B6C0 

676 
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172. N=6 

K=1 

GO TO i7k 

173 H=6*HE-6» (IPHO-A3I 

665 17i< «=N-5 

HCOH=H-l+6»LOS 

IF(NCON .GT.256) NCON=258 

DO 175 I=M, N 
SQR=l./SQi?T(A (I, in 
670 DO 176 J=l,NCON 

_ A{I,J)=A(I,J) »SQRl 

176 CONTINUE 

A<I,LUU = A(I,LLLn;.SO^ 

” ■ A<I,I)=SQR 

^ 675,.„„. IF(I-NI 177, 180, 180 

177 IPl=Ii-l 

00 175 L=IP1,N 

00 193 J=L,NC0N 

A<L,J) = AtL,JI-A(I,U»-ACI,JI 

660 193 CONTINUE 

ATLiULl ) = Aa,LLU-A(I,U*AtI,LLL» 

175 CONTINUE 

„ 160 IF(KI 178, 178, 181 

I=H,N 


B600 

677 

S603 

678 

8600 

679 

3600 

660 

B6CQ 

661 

BbOO 

662 

8600 

683 

fl&OO "■ " 

684 " ■ 

3600 

685 

B600 

666 

6603 

687 

B600 ' " 

688 ■ 

b'600 

68 9 

0600 ■ 

690 

0,600 

691 

B605 

692 

B605 

693 

8600 

694 

B600 

895 

8600 ' ■ 

696' ■ 

B600 

697 

B600 ■ ■ 

696 

6600 

699 

BpOC 

7G0 


685 


690 

C 


695 


Veo' 


705 


710' 


DO 186 J=I,NC_0N 

AtJ,I) = At I, j) 

186 CONTINUE 

A(LLL,I)=.AtI,LLL) 

320 CONTINUE _ 

CALL HRITHS(12,A,155A,ME»‘ - 

EFFECT OF RESET PLATE ON <3HS OF FOL LOHIN G CAM ERAS 
178' IF (KE-I PHO) 18 27179, 17 9 
182 NEIL = NH 

00 183 I=H,N 

IFtl-NJ 18l,,185,185 

185 IF <K) 16<,,18A,305 

305 DO 196 L=NEIL,NCON 

■ '■ DO 197 J=L,NCON 

ATU-6,J-6I = A(L, J) -A (L ,I1»A C J,I1 
197 CONTINUE ’ ' 

A<L-6,LLU = A<L,LLLI -A(L,I) *A(LLL,I) 

196 CONTINUE 

J0YCE=NC0N-5 

DO 298 L=1,NC0N 

00 250 J=JOYCe,NCON _ 

■ AIL, J)=C. 

298 CONTINUE 

00 299 L=JOYCE,NCON 

AIL,LLL)=0. 

299 CONTINUE 
GO TO 183' 

ia<* DO i9<4 l=neil)NCon 

DO 195 J=L,NC0N 

AIL,J) = ail, J l -A (I ,L) * All, J} 

195 CONTINUE 


OcCO 7G1 

'B6C3 '■ ■ 702 

8600 703 

B6B0 704 

8600 705 

P600S ■ ■ 9 " 

0600 707 

8600 "" 7E8 

B60C 7G9 

B600 710 

B600 711 

86033 ■ 10 ■ 

i)600 717 

B600 ' • 718 —■ 

8600 719 

* B60J 720 

B6C0 721 

B60J 722 

8600 723 

B60u 724 

8600 725 

“ " B600 726" 

0600 727 

B6C0 728 

B6D0 729 

8600 730 ■■ 

B600 731 

8603 732 

B600 733 

B600 734 

8605 735 


715 
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194 

183 




IIP 

K 


725 


725 




179 

198 


730 


J.99 

215 

225 


735 


186 

189 


740 


745 


221 

217 


750 


218 


219 


755 


760 


222 

"216 

220 


765 


223 

307 


770 


A{L,LLL) = Aa,LLU-A I.LLU _ . 

cowriNue 

CONTINUE 

H£=ME»1 

KC=0 .... 

GO TO 105 

SACK SOLUTION OF CA HERA_lOHS_ANO AOO RESULTS JO CAM ARRAY 

‘SCALE=0AM (1,5J 

IF(rPH0-43»198,199,199 
K=JOVCE=6*IPHO 
LAST=IPHO 
GO TO 215 

K=JOYCE=256 

LASr = 43 ~ ‘ " ■ 

SAVE=1./A (JOYCE, JOYCE) 

PRINT 556 , HERAT, A( JOYCE, -LU , SAVE 

ME-1 _ 

DO 216 NU=HE,LASI 

NROM=IPHO»l-NU 

IFIDU. GT. LOS) J0YCE= JOYCE-5 

IF(JOYCE-6»LOS)188,189,183 

J0YC£=6»L0S 

00 216 1=1,6 _ 

IF(NU,EQ.1.ANO.I.EQ.1)CO ro 217 

H=K 

K=K-1 

DO 221 J=M, JOYCE 

A(K,LLL» = AIK,LLL) -A{J,LLL» » AIK, J1 

continue • _ 

A tK,LLL ) = A(K,LLL) »A (K,K) - . - . 

IF (I -3) 218,218,219 

J=6-I 

CAHINROH, JI=CAH(NROH, J1 4-A ( < , LLU *SC AL E 

GO TO 216 
J=12-I 

FK=SQRT(1.-A1K,LLU»A(K,L.I.)) 

SAVE=CANtNROM,J>»F< f C AH( N ROM , J*3 » »A (K,LL 
CAHCNROH, J4-31=CAK(NR0H, Jf3) »FK -'CAM(NR0W, 

CAKCNROW, J)=SAVE 
IF(NU.EQ.1.ANO.I.£Q.4)GO TD 222 
IF(ABS(A{K,LLU )-ABS(G( 12 JO 1))> 216,216,222 
NHI=NR0W 

G <12001 )=A(K,LLL) 

CONTINUE 
HE=LAST=LASTfi 
IF(HE-IPHO)220,220,227 
NROH=IPHO ■H-ME 
CALL REA0HS(12, A,1554 ,NRD<() 

DO 223 J=7,256 

‘l=265-J 

A(I,LLL)=AtI-6,LLU 
CONTINUE 
DO 224, 1=1,6 
DO 224 J=I,2S9 
A(I, J) = A( J,I) 


U 

J)'*AfK7LLr>' 


B60C 

736 

8600 ■■ ■ 

737 '■ 

B60J 

738 

B600 

739 

8640 

740 

B600 

741 ■ ■■ 

8603 

742 

8600 

743 

B6G0 

744 

6600 

745 

8600 

746 

8600 

1-^7 

8600 

748 

8603 ' “ 

749 '* 

8600 

750 

8600 

751 

G600 

75 2 

B600 

753 

8600 

754 

8600 

755 

8600 

756 

6600 

757 

8600 

758 

86 0 0 

759 

8600 

760 

8600 ■■ *' 

761 

8600 

762 

8600 

763 

8600 

76 4 

8603 ' ■ 

765 " ■ ■■ 

8600 

76 6 

B6C0 

767 

8600 

768 

0600 

769 

8600 

770 

8600 

771 

8603 

772 

8600 ■" * ■ 

773 

8600 

77 4 

8600 

77 5 

6600 

776 

?b00 

777 

8600 

770 

06 00 

779 

860C 

780 

B6C0 ■ 

781 

B6DQ 

78 2 

B60'OS 

11 

B60C 

784 

8600 ■ 

785 

8600 

786 

8600 

787 

660 0 

792 

860j 

793 

8600 

794 
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775 


780 


785 


790 


795 


8G0 


805 


810 


.6 15 


820 


224 CONTINUE , 

K=7 

GO TO 225 
227 IS=CAM(NHI,1) 

PRINT S33,ITERATtG(120011»IS , . 

IF (ASS «G (12001) 1 -.000 01»2JO »230 t228 

_228 ITERAI-=ITERATH . 

1F(IT£RAT-HAX1 229’.'229.230 

229 IRE=1 . 

GO TO 231 

230 1RE=2 ' . . 

C 

C_ CLOCK ORIENTATION SOLUTION COflPLETEO 

■ C “ 

PRINT 541 

DO 232 I=1,IPH0 

IS=CAH(I,1> 

PRINT 561,lS,tCAHa,J),J=3,5) 

_ 232 CONTINUE 

PRINT 539 

DO 233 I=1(1PH0 

IS=CAM(I.1) 

PRINT 557 ,IS,ICAM(I,J),J=5,8) , , _ 

PRINT 557,15, CCAMa,J>,J=3, 11) 

233 CONTINUE . 

“"print 547 " 

PRINT 534 
PRINT 548 
C 

C INTERSECTION OF O0JECTS U3EO IN BLOCK ORIENTATION SOLUTION 

C 

Z3l”lNSECT=:l 

00 329 1=220,225 
00 329 J=l,20a 

A(I,J)=0. „ . ^ 

32'9 CONTINUE 

KC=3 

IT£ST=0 “ 

CO TO 234 , 

360 IF(NA.EQ.Q) GO TO 268 

C SET UP IMAGES OF OSjECT TO BE CONPJTEO DV INTERSECTION 
355 IHAG£=1 

361 A(1,IMAGE)=0(1,N) _ 

"■ A(2,IMAGE)=8(2.N) 

_ A(3,I«AGE>=8(3,N) 

270 J=A(i, IMAGE)*. 000001 
FK=J 

00 235 L=l,IPHO 

IF(FK-CAK(L.in235,236,235 

235 CONTINUE 

236 A(4,IHAGE)=L 

A(5,IKACE)=20*(L-1) - 063M (L-D/43) 
rF(LARRV)358.35S,357 
“357' A(6, IMAGE )=CAN(L,21 
CO TO 359 


0600 

8600' ' 

8600 

8600 

B6Q0 

8603 

BoGJ 

8500 " 

8603 

860J 

B6ca 

8600 ' 

B6G0 

8660 

8600 

()663 

8600 

A600 

B60C 

"■ 8600 ■' 
(36ti0 
860 0 
^603 
8603 
3600 
36QG 
8600 
8605 
8600 
B60u 
86C0 
8600 
8600 
(360 0 
3600 
8600 
Q60C 

"8600 
8603 
8600 
8603 
6600 
8600 
8600 "■ 
B60G 
B600 
B600 

' B60D 
B60J 
B6D0 " 
8600 
8600 
B600 
B60C 
8600 


795 

796 

797 

798 

799 
SCO 
ecl 
oG2 
833 

634 

635 
80 6 
60 7 
60 6 
8C9 
810 
611 

6 

813 
" 614 

615 

616 

617 

618 

819 

820 
821 
822 
823 
82 4 
625 
82 6 

627 

628 

629 

630 

631 

832 

833 

834 

835 

836 

837 
836 
839 
84 0 
841 
&hZ 
843 

■ 8h4 

845 

846 

847 

848 

849 


825 
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358 

* 

359 

830 „ ,. . 

238 

835 “ 

246 

237 

840 

287 
2 68 

241 
C 

242 

645 


. 

124 

850 




289 

855 


8 60 

239 

323 

363 

865' 

240 

870 


A (6,IMAGE)=FL 

FK=J»10 OOOC-0 

A (7,IHAGE)=Atl,IHAC£) -FK 
IFdOTHER.EQ.U GO TO 313 
IF(A{7,1) -A (7, IMAGE »> 237t’I0.237 
M=N+1 

IHAGE = I KAGE + 1 . 

IF(N'NA/3)361,361.2A6 

KC=1 

IrtAGE=IHAGE-l 
IS=A<T, II 

P Hit 1> =P<11»31 =0 . 

00 207 NU=ltLlNE _ 


875 


630 


2AE 
.3 62 

311 


2<.g 


IFI At7. ll-G(6*NU-5> )287,2i9 ,287 
CONTINUE 

GO TO <241,242) ,IR£ 

IF(G<6»NLI-1I. NttC.O, AMD. G(6»NU) .N£, 0.01 GO TO 267 
P COEFFICIENTS, CONSTANT TERMS, CONTRIBUTI OSS TO NORMAL E(3l)Al 

DO 124 1=1,3.. 

00 124 J=I,4 

E<I,J)=0. . . 

CONTINUE 

DO 362 l'=l, IMAGE 
K=A(5 ,M) 

IFdOTHER.Ea.l) GO^TO_323 

’DO 269 1=1,3 
J=6*MU-5i-I 

C (I 6) =C(J)-C <I,K«-11 

CONTINUE . . .... . . 

J = Kf4 

DO 239 1=1,3 

C<I,K«-17I=C(1, JMC(1,K+16I 4C(2, J) *C(2,Ki'15)i'C(3,J>*C<5,K*16) 

J = JH 
CONTINUE 

IFtlOTHER.EO.OIGO TO 363 

IF<A(10,HI .ER.O.OIGO TO 362 

DO 246 J=1 ,3 .... - . . .. ..., 

P (1,J) = <A (2,M) »C(J,K*6» +A(5 ,Hl*C(J,K«-4H/0 (J ,K417) 
P(2,J)=(At3,M)*C(J,K«-6M-AC6,H)*C<J,Kl-S>)/3 (3,K«-17) 

CONTINUE 
L=A(4,H) 

P tl,4l = P<l,ll»CAH(L,3l+P(l,2) •CAH(L,4 ItPU,3 )*CAH(L,5) 

P(2,4l = Pt2 ,1) *CAK(L,3) ♦•PIE , 2) *CAmU,4H-P<2 ,T 1_^CAHIL,5)^ 

"do 245 1 = 1,3 " - 

00 245 J=I,4 . 

DO 245 L=i,2 

E(I,J) = E<I, JM-P(L,II »PIL,JI 
CONTINUE 

CONTINUE 

SOLUTION OF NORMAL eQUATIDNS' FOR' XYZ 
00 248 1=1,3 
SOR = l./SftRTIEIItin 
00 249 J=I,4 
E(I,J1 = £H, J)»SCIR 
CONTINUC 



8600 

85 0 


8606 

851 


B&Ou 

852 


8600 

853 


8603 

854 


8600 

8S5 


8600 

856 


B60J 

857 


8603 

858 


Q6CQ 

859 


8603 

660 


8603 

861 


86.00 

862 


86uu 

863 


0600 

664 


8600 

865 


86G3 

866 

3HS 

8600 

867 


8606 

868 


96l0 “ 

869 


BbOG 

870 


8600 

871 


8600 

872 


8638 

■ 873 


8603 

874 


86,00 “ 

” 875 


8600 

076 

■ 

8600 

877 


6600 

878 


e>30 

879 


B60.J 

880 


8600 

881 


8600 

862 


8600 

883 


B600S 

12 


660QS 

13 


86035 

14 


■ 8500 

886 


8608 

887 


360G 

686 


86Uu 

869 

“ 

8608 

890 


86U0 

891 


8600 ■ 

692 


8600 

693 


8600 ■ 

894 


8600 

895 


8630 

896 


8600 

597 


' 8608 

898 


86C0 

899 


8603 

9C0 


8600 

9C1 


6600 

902 


8603 

90 3 
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20.A0.5A- 



eas 


ago 




895 


9 BB 


9CS 


910 


91S 


920 


925 


930 


935 


e(i,n = sQR 
IFa-31 250,251,251 

250 IPl^IK 
00 248 M=IP1,3 
00 248 J=M,‘* 

_248 COMTINUE 

251 t(3,4l=E(3,4)*E(3,3r 
00 252 1=1,2 
NMI=3-I 
NMIP1=NMI«-1 
00 253 J=NM1P1,3 

E(Mm ,4 ^ = E(NKI,4| -E (J,4)»E( NHI, J) 

"esY continue 

E(NHI,4> = E(NHI,4I»E«NNI,N'1I I 

252 CONTINUE 

IFaOTHER.EQ.l.ANO.IGO.EO.l I GO 70 274 
IFaOTHER.hQ.l.AHD.IGO.Ea. 21 GO TO 276 

test to TERMINATE OOJECT INTERSECTION SOLJTION 

SAVE=AE}S(Ell,4) -G (6-»NU-l) I ' < -A 1 6 , IMAGE) /C ( 3, K»17> > 

IF( SAVE-. QOOOOIJ 257 ,257 ,255 

'257 SAV£=AQStE(2,4)-G C6»NU-3J I » l-A(8, IHAGEI/Ct 3, K<-17> ) ■ 

■ IFISAVE-. 0000011259, 259, 253 
258 NEIL=2 

GO TO 273 

259' NEIL=1 

273 GO TO (243,264) ,IRE 

OLOCK ORIENTATION SOLUTION NOT COMPLETED - , UPDATE ONLY NONOONTROL 
XY2 FOR ITERATION OF THE 3.0CK ORIENTATION 'SOLUTION 

243 IF (G<6'‘NU-i> ) 260 ,244,260 

244 G(6*-NU-4)=Gtl,4) 

'g(6*NU-31 =£(2»4) ‘ 

260 IF(G(6»NU) >262,261, 262 

261 G(5*NU-2)=E(3,4) 

262 GO TO (267,242) , NEIL _ 

BLOCK ORIENTATION SOLUTION COMPLETED 

264 IF(A(7,1) .£Q.G(6‘LINE-5)) ITEST=l 

IF(KD) 280,280,265' ' - - 

280 IF(G(6*NU-1) .NE.O .0 .AND.3(5‘>NU) .NE.fl.OI GO TO 282 _ _ 

IF(G(6»NU-1) .NE.O .0 ) GOT! 261 

IF(G(6»NU) .CO.fl.OJ GO TO 338 
G(120D3)=G(6»NU-2) 

PRINT 560.1S.G(1 20B3 ) 

'go to 2 84 

PRINT 562 _ _ _ 

GO TO 2 84 - . - . . . 

C(12C01 )=G(6*NU-4) 

G(12D02)=G(6»NU-3J 

PRINT 554 , IS, G(1200i) ,6 (123 02) 

GO TO 2 84 ’■ • - ■ 

PRINT 5 54,I5,G(6»NU-4) , G( 5 NU-3 ) t Gt6* NU-Z) 

CALL OOUFU) 

£NCODE(80 ,561,BUF (IBUF) ) I S , G ( 6» NU-4 ) , 6 <6»N U- 3) ,G(6*NU-2) ,ITEST 
KO=l 

GO TO (346,265) ,NEIL 


30 8 
281 

282 


284 

285 


B60G 
B 6 C 3 
B 6 02 
B6C0 
B60a 
C6GS 
B60a 
B60()" 
B6C0 
B60j 
6600 
B680 
B60C 
8600 
B6C: 
8600 
B6u0 
8600 
B50G 
B6C0 
B6CG 
B600 
Q6D0 
B&OC 
0600 
S60u 
8600 
0600 
B6GJ 
3600 
B600 
06 Ov 
B 6 CG 
8600 
0600 
' B6C0 
0600 
360C 
B600 
86 C u 
8600 
(3600 
8600 
■0600 
B600 
360 0 
8600 
0600 
8600 
8630 
8600 
0600 
B&oa 
B600 
B600 


90 4 
905 
9 G 6 
907 
90 8 
9C9 
910 
9X1 

912 

913 

914 

915 

916 

917 

918 

919 
92 0 

921 

922 
92 3 
924 
92 5 
926 
92 7 
928 

92 9 ■ 

93 0 
931 
93 2 
93 3 
93 4 

935 

936 

937 
93 8 

93 9 

940 

941 

942 

943 ■ 

944 

945 

946 

947 

94 8 

949 

950 
551 

952 

953 

954 

95 5 
95 6 
957 
95 8 
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PAGE 


o 

w 

Y 









.. 940 . . 

266 

945 

290 1 

^50 . 


♦ 

294 

1 

955 

292 


291 
. ZG3 

960 




1 

965 


970 

286 



C 

267 

975 

268 


269 

980 

272 

. 


985 


990 . . 

348 


G<6*NU-4)=E(1,41._ 

G(6*NU-3>=E<2»4) 

C(6»NU-2)=E(3,4) 

GO TO 242 
00 285 J=l, IMAGE 
K=A(5, J) 

P(13 1 11 = lA (2, J) +A (&,J) *ca , K + 17)/C (3tK*17J l*10COOCO. 

Ptl3,2l = (A(3,J» tA (o,JI»C«’, K417 )/C{3,K»17l I *10 0000 0 . 

IPl=A(liJl 

IF«J-l)291.290,291 

PRINT 5 53 ,IS.£a,4»fE(2.4).EC3,4I.IPl,P<13.11,P(13,2) 
IFlGt6»NtJ-ll .NE.O.C. AND. Gt5*NU) .N£. 0.01 GO TO 2G3 

CALL OnUF ill . 

IF(G(5*NU“1) .NE.S.ai GO TD 292 
IF(G(6*NIJ1 .NE.O.Ol GO TO 294 

ENCODE I 80 ,961 iCJUF (lauFI 1 rS,E(l,4J,EI2,41,E{3,4),ITEST 
GO TO 263 

ENCODE <80 ,E61,0UF<rBLIFI I I S , £ (1 , 43 ,E (2 ,4 I , G t 1 200 31 , 1 TES T 
GO TO 263 

ENCODE ( 8C , 561 (DUFdOJn 1 I S ,G (1 200 11 , C <129 02 1 , £ C3 , 4) , ITEST 
GO TO 263 

PRINT S36,IP1,P<13,11 ,P<13, 2) 

P<12,1)=P(12.11*P<13,11*P<13,11 *P<13,2I*P<1S,21 . . 

P<12,31-P<12.3»T2, 

P<11,11 =P<ll,l>TP<13,JL.l *iLl-3ai TP (13, 21 »Pt IS, 21 

P(U, 3) =P<li, 31+2. 

L=A(4,JI 

IF<L.L£.2Q0) NROW=220 
IF<L,GT . 200 .AND.L.Lt. 40 01 'JR0H=222 
IF(L.GT.4C01 NR0H=224 
K=L-2D0 *< (L-1 l/'2B01 

A< NROW.K) =A<NR0H,K1 tP < 1 3 . L I'^P'd 3 , 1 1*P <13,'2 r*P(13, 21 
A<NROWfl,K)=A<NROm-l,K) *2. 

CONTINUE 

P<U,41=GQRT(P<U,11/P(11, 3 II 
PRINT 57D,P(11,4> 

K0=O - , 

RECYCLE FOP, NEXT OBJECT NEXT PLATE 
IF<KC.£O.01 GO TO 355 
ME=MEH 

KC=0 . . 

IF (ME-IPNOl 10 5,10 5,2o9 

GO TO <271,2721 , IRE _ 

P<12,4l=SORT<P(12,ll/P(i!,311 
PRINT 538,P<12,41 
PRINT 578 
DO 348 L=l.IPHO 
IS=CAH<L,1) 

lF<L.LE.20fll NROH=220 _ 

IF(L.GT.20Q.ANO.L.LE.4001 NROH=222 
IF<L.GT.40S1 NR0W=224 
K=L-200 » < (L-1 1/200) 

A(NROH, KI=SQRT (A (NROH,K)/A ( NROHfl.KI 1 
print 557,1S,A(NR0W,K) 

CONTINUE 


BGOO 

B&oa 

B60U 

B6DQ 

B&QO 

0600 

B600 ■ 

B60D 

S6QU 

3600 

BbO 0 


18 


B60Q 

959 

B600 

960 ' 

9680 

961 

0600 

96 2 

860o 

963 

B60G 

964 

G&QO 

965 

6600 

966 ' " ■" 

3600 

967 

B600 

968 

3600 

96 9 

B6J9 

970 

8605 

971 

8600 

972 

B&OS 

973 

B&GO 

97 4 

6600 

97 5 

0608 

976 

6600 

977 

B6QO 

978 

B600 

97 9 

3600 

980 

B60l) 

961 

0600 

38 2 

B60a 

98 3 

B6CC 

'984 ■ 

B600 

98 5 

B685 

98 6 

3603 

98 7 

B600 

98 8 "■ 

B60C 

98 9 


99 0 
991 
99 2 
993 
99 4 
99 5 
• 59 5 
99 7 

998 

999 

io:o 


B680 

lOCl 


B680 ■ 

10C2 


B60a 

1053 


B600 

lCi.4 


6SC0 

1005 


B6GJ 

100 6 


3683 

10C7 


B600 

"IOC 8“ 


B600 

1009 


B6Q0 " 

laiu 


3608 

iOll 


B6Gu 

1012 


8600 

1013 
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o, 

« % 

TO 


<o 


995 . .. 

C 

. c 

" " 1000 

. c 

316 


313 

.. 1005 

314 

I ' 

. . 315 

IDifl 

- - „ 

1015 

1 

1 

- 

i 

i 

328 - 

1020 1 

1 


1 

1 

1025 

1 

( 

i 

1 

1 

1030 

327 1 


C 1 


1 

E7«t \ 

< 

1C35 1 


C ( 


\ 

j 

1040 

1 

) 


. 326 ( 


1045 


INTERSECriON OF OBJECTS NOT USED IN THE BLOCK ORIENTATION SOLUTION 

PRINT 530 
PRINT 534 
PRINT 548 

IOThER=IKAGE=ICO=HORT=1 
JOyCE=LINE=0 

READ IN IMAGES OF THE OBJiCT 
call OOUF a ) 

DECODE t 80 ,526,BUFF( IXdUFH ( A ( J, IMAGEJ ♦J=l, 3) 

IFIAtl, IMAGE). EQ. 0,01 GO T3 315 
CO TO 270 

IFUC7, 1) -A (7, IMAGE) ) 315,314. 315” " ‘ 

IHAG£ = IHAG£tl 
GO TO 316 
IMAGE=I HAGE-1 
IFUMAGE.Ea.O) GO TO 350 

IS=A(7, 1) 

PRINT 5 62 

DO 327 H=l, IMAGE 
A<i0,M)=l. 

K=AT5,H) 

L=A(4,h) 

00 328 J=l,3 _ 

CU.KFl ) = CAM(L, Jt2)' 

C( J.Kf2)=CAM(L,Jf5) 

C<J,Kl-3 ) = CAM(L, J+8) 

CONTINUE 

C(l,K^4) = CI2.K«-3) »C(3‘,K+3) 

C(l,K*5) = -C)2,K«-3) *C(3,Ki-;) 

Ca.K»6J = C(2,K+2) '■ 

C(2,K<-4) = C(1,KF3)*C(3,KF2) f 
C(2,Kt5)=C(l.Kf3)»C(3,K*3) - 
C(2,Kf6)=-C(l.K<-2) •Ct2,K*-3) 

C{3,K«-4) = C(l,Kt-2)*C(3,K+2) - 
C(3,Kv5) = C(l,Kf2)»C(3,K + 3) » 

CI3)K*6) = C(1,K*'3)*C(2,K»3) " 

C13,K<-17) =1, 

CONTINUE ■ ■ ~ 

COMPUTE L.S. SOLUTION OF XTZ USING A3B EQUAL 1 
GO TO 242 ...... 

CU2001 )=e(l,4) 

G<120C2 ) = £(2,4) ' 

G(12003)=£(3,4) 

COMPUTE L.S. SOLUTION OF <i Z USING COHPUTtO A30 ' ’ - - ~ 

IGO=2 

DO 325 H=l, IMAGE 
K=A{5,H) 

DO 32 6 1 = 1,3 ' ■ " 

C(r,K+16) = £(I,4)-Ca,K+l) 

CONTINUE 

J=K»-4 

DO 325 1=1,3 

C(I,Ktl7)=C (1, J) »>C(l,K+lo) tC(2, J)*Ct2,K+l9)>C<3, J)*Ct3,K+15) 


C«1,K*2) *C(2, K»2)*Ct3,K»^3) 
C(1,K+2)*C(2, K»2)*Ct3,K + 2> 

C(2,K»2) *C( 3, Kt-3) •CCl,K<-5> 
C(2,K«-2) »C(3, Kf2) *Cll, K<-3) 


Q60U 
B60O 
0600 
0600 
8600 
0600 
0600 
B 6 DG 
6600 
860u 
0600 
0600 
0600 
0600 
oeoj 
0600 
B6G3 
B 6 GG 
B6GC 
BdOC 
86G0 
B600 
0600 
3600 
8600 
8603 
0608 
0600 
0600 
G6G0 
0600 
8603 
8600 
8600 
0600 
8600 
8609 
B 6 C 0 
8600 
0600 
8600 
8600 
B60e 
8600 
8600 
8606 
6600 
8600 
8600 
B6DC ■' 
8600 
8600 
8600 
0600 
8600 


1014 

1015 

1016 
1017 
lOlS 

1019 

1020 
1021 
1022 

1023 

1024 
1C25 
1C26 

1027 

1028 

1029 

1030 

1031 

1032 
1C33 
1034 
1055 
1C36 

1037 

1038 

1039 

1040 

1041 

1042 

1043 
1C44 
1C45 
1046 
1347 
1048 
1C49 

1050 

1051 — 
1852 
1D53 ■“ 
1054 

less " 

1G56 

1C57 

1058 
1G59 "" 
1C6Q 
1061 "" 
1062 

1063 

1064 

1065 ■■■ 
1C66 

1067 

1068 
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ZO 



325 

„ 1050 

C 

276 


330 

1055 

C 

. _ 331 


.347 

.. 1060 

- ■ 

- 

332 

C 

1065 


277 

1070 

278 


319 


279 

1075 

275 

C 

226 

1060 

333 

335 

334 

1085 

C 


336 

1090 

337 

— 

338 

339 


1095 


1100 


3A1 

342 


B600 

coNriMue 8600 

IF(IOrHER.EQ.O) GO TO 266 8600 

GO TO <2A2,347> ,MORT 8600 

TEST TO TERMINATE OBJECT INTERSECTION SOLUTION B60J 

SAVEsAOSIEdiD-GtlZOOD* * < -A ( 8 , IHAG E) /C 1 3 , K +171 ) B60C 

IFISAVE-. 000001)330, 330, 2rL _ 8600 

SAtfE=ABStE[2,G)-G<12002)J » (-A t 6 , IMAGE ) /C ( 3, K Hyi 3 ' 6600 

IF (SAVE-. 0 00 0 0 1) 3 31, 331, aft B60 0 

TERMINATION TEST GOOD. CUPUTE PLATE RESIOJALS ’ 8600 

MORT=Z 8600 

GO TO 346 8603 

DO 332 H^l, IMAGE _ _ _ 8600 

IF(A(1C,M) .EO.0.0) GO TO 33Z ' ' ' ' - 

K=At5,H) 8600 

A{8.HJ = tA(Z,H)*-A(6.HJ*C<l,<H7)/C(3,K*17n»1000 00 0. 860 0 

A«9,H)= £A (3,M) tA(6,H)*C(2, <H7> /C(3,K«-17) ) »10 OCflflO. 850 0 

CONTINUE 8600 

DETERMINE MAXINUH_PLArE R-SIDUAL AND TEST AGAINST RESIDUAL LIMIT D6CC 

N = d Bo U V 

H=J=1 3600 

DO Z75 1=8,9 ” 

MsHU BSuli 

iFrM-rMACEjzra.E/atzrB ' O6oo 

IF(AaS(A(N,J>)-AOS(A(I,H>) >319, 277 ,27 7 B6C0 

’N=I - -- g-gj 

J=M B600 

GO TO 277 ■ ■ B60G 

H=0 _ _ 860 Q 

CONTINUE . 

IF(ABS( A(N, J) ) -REST D) 165.15 6,22 6 360 U 

QAO IMAGE ' “ ‘ *' ' ' - ~ jjggj 

JOYCE=J(3YCEd 8600 

IF(J0YC£-i)333,333,334 , ’ B60d 

IF(IHAGE-3)334,334,335 3600 

A£8,J)=A(9,J»=A(10,J)=0. ' “ - -- 

MORT=l 8600 

"go to 274 ■ ” B60C 

A(10,J)=2. _ 8606 

OUTPUT XYZ AND PLATE RESIOJALS 6600 

IP1=AC1,1) BcQO 

IFtAdO, 11 -1. )336, 337, 338 B600 

PRINT 549, IS,£(1,4» ,E(2,4) ,£13,4) ,IPt B6Cu 

GO TD 339 '■ ■ - - - B600" 

PRINT 553,IS,E(1,4),£(2,4) ,E(3,AI ,IP1,A(B,H ,A(9,1) _ 8600 

GO TO 339 ~ ' ”■ 8600 

PRINT 550 , IS,E d»4) ,E<2,4) , E<3, 4) ,IP1, A(8,l) , A(9,l) 3600 

CALL OBUFd) B600 

eNGOOE(aD,561,BUF(IBUn) IS, Ed, 41 ,£(2.4) ,£(3,4) 8600 

DO 340 H=2, IMAGE “■ " ’ ' "8600 

IPi=Atl,M) B6D0 

IF(AdO ,M) -1. ) 341 ,342 ,343 8600 

PRINT 551, IPi B6C0 

GO TD 340 86Cu 

PRI!JT 636,IP1,A(B,M) ,A(9,'1) B6CC 
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1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1078 

1079 

1080 
1081 
1082 

1083 

1084 
1035 
1086 

1087 

1088 
1C89 

1090 

1091 

1092 

1093 

1094 

1095 

1096 

1097 

1098 

1099 

1100 
llCl 
1102 
11Q3 
11C4 

1105 

1106 
110 7 

1138 

1139 
1110 
1111 
1112 ‘ 

1113 

1114 

1115 

1116 

1117 

1118 

1119 

1120 
1121 
1122 
1123 
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31.3 
3 40 
3 44 


1105. C 


345 


1110 


1115 


1120 


1125_ 


350 


3 65 
3 66 


1130 


GO TO 340 

’print 555,'iPl,At8THUA't9,'t 
CONTINUE 
LINc=I-IN£vi 

RECTCLE FUR NEXT OTHER OBJ 
IF(A(1, EHAGE+iI.EQ.O.D) 3D 

00 345 1=1,7 

A {I,1) = A(I,IHAGE+1> 

CONTINUE 

JOYCE=0 

IGO=IHAa£=MORT=l 
GO TO 314 

BLOCK ADJUSTMENT “completed" 

PRINT 552, LINE 
IF(INVER.Ea.O) GO TO 366 
00 365 I=l,IPHO 

IS=CA«(I,1} 

CALL OBJFdl 

£NCOOE<80,561,SUF(rBUFn I 
CALL OeUF(l) 

ENCODEC80 ,557,BUF(I6UF)) I 
CALL OBUF(l) 

ENCODE (50, 551, BUF(IBUF)) I 

CONTINUE 

PRINT 542 

CALL OBUF(O) 

GO TO lOGG 

FORMAT STATEMENTS 


EOT 

TO 350 


3 ,ICAM(I,J) ,J=3,5» , INVER 
S , tCAKa, Jl ,J=6, 5t 
3 , (CAM(I,J) ,J=9, lU 


1135 


1140 


1145 


1150 


524 F0RMAT(lX,ll,lX,ei4.7,2(J< 
'525 FORHAT(7X,F3.0,3(2X,Fi.0), 
526 FORMAT(1X,F9.0,2<2X,E14,7I 
”527 FORMAT(4X,F6.0,3<1X,F12,31 
528 FORMATtFSaH GROUND COORDH 
'529 F0RMAT(44H PRDVISIOi(AL GRD 
ICTS. EXCEEDS 2D00 LIMIT-- 

530 FOP.MAT(lHi,28X,56H»»'f H 
ION ***) 

531 FORHATCaAiO) | 

532 FORHAT(/SOH ORIENTATION ®A 
"533 F0RHAT(14H PROGRAM PASS U' 

lECTION OF E14.7, IIH FOR = 
■534 FORMAT<3ax,37H (SECANT PLA 

535 F0RHATI/7H PLATE 13, 36H V 

536 FORMAn48X,I16,F16.3,F19.J 

537 FORMAT{/8K OBJECT 16, 31H 

'538 FORHAT(//733H RMS RESIDUA: 

539 FORNAT(/52H PLATE 

541 F0RMAT(/44H PLATE 

542 F0RMAT(/54H BLOCK A 

■ 543"F0RNAT( 716H BLOCK CONTAINS 

544 FORMATl/arH BLOCK AOJUSTM- 


3X,£14.7> 

,37X,I1> 

,30X,I1) 

ATES MISSING FOR 16, SH-STOPF) 
JND COORDINATES ENTERED FOR 14, 
STOP/} 


raheter correction 
46H PRODUCES A MA 
-ATE 13) 

NE COORDINATES * M 
EEOS OVER 5 RESECT 
J 

APPEARS ON OVER 9 
■ FOR ENTIRE BLOCK 
OMEGA P 

XO YD 

DJUSTHENT 
I3,17H PHOTOGRAPH 
NT CONTAINS 14, 41 


LIMIT IS 0 .C u001/> 


ICRONS)//) 

ION PASSES--STOP) 

PHOIOS--STOP) 

=■ £14,7, 6K 'rflCROT 
HI KAPF 

70/) 

C 0 N P L E T 



B6CQ 

1124 


B66C 

1125 


B6C0 

1126 


0600 

1127 


B6G0 

1128 


B600 

1129 


3600 

1130 


BoOO 

1131 


Q&CG 

1132 


3600 

1133 


36 Gu 

1134 


B6Q0 

1135 


36G0 

1136 


8^,0 0 

1137 


8600 

1138 


3600 

1139 


A60S 

7 


A6Q0 

8 


A6DG 

9 


A60D ‘ ■ 

10 


A6Q0 

11 


A6GD 

12 


A6Q0 

13 


A603 

14 


A6Q0 

15 


A60C 

16 


A6Q0 

17 


3600 

1141 


360 0 

11..2 


36'QO 

1143 


360 0 

1144 

' " 

8600' " - 

1145 

> 

A603 

18 


8600 

1147 


8600 

1148 


860 0 

1149 


860 0 

1150 

)BJE 

86u0 

1151 


3600 

1152 

,UTI 

860 u ■ 

1153 


B6DJ 

1154 


0600 

1155 


3600 

1156 

;ORR 

"860 0 

1157 


860 U 

1158 


0600 ■ 

1159 


3600 

1160 


860 0 

1151 


860 C 

1162 

u 

8604 

1163 

i/) 

8600 

1164 


8600 

1166 

■ 0) 

8660 

1167 

;/) 

8600 

1168 

:os 

3600 

1169 


1155 



. PROGRAH 


86C0 


GDC 6&30 FTN V3.0-32^ OPT=2 Q'*/02/75 20.<»at54 


ILIMIT OF 1978 — STOP) ' •' ^600 

S'*? FORHAr(/ 7 H PLATE I 3 t EgH" EXCEEDS 'SQ-IHAGE LIMIT — STOP) ' ” ’ "■ C 600 

, „ 546 F 0 RHATI 24 H » DISK tPARITT) ERROR ») 6600 

547 FORHAniHliSOXfSEH*** INTERSECTION OF 03 JECTS USED IN SOLUTION B 60 Q 

1160 . , I-***) 66 oa 

548 FORHAU* OBJECT X SiDUNO Y GROUND Z GROUND 6606 

1 IMAGE VX PLATE RESIO VY PLATE RiSIO *) B 6 G 0 

549 FORHAT<lX,ir,F 14 . 3 , 2 F 13 . 3 , [ 16 , 28 H OVER LIMIT — IMAGE DISCARDED)"" 0600 

, . 550 F 0 RNATtiX,I 7 ,FI 4 . 3 , 2 F 13 . 3 ,ri 6 , 2 H • , FI 4 . 3 , - 19 . 3 ) 66110 

1165 551 F 0 RMAT( 4 BX,I 16 , 28 H OVER LMIT— IMAGE OISCAROEDl 8600 

552 FORHAT</'/ 40 H TOTAL NUMOER OF INTERSECTED OBJECTS IS 14 ) 8600 

553 FORHAniX,I 7 ,F 14 . 3 , 2 F 13 . 3 , I 16 ,F 16 . 3 ,F 19,31 8600 

554 FOR«AT(/ 7 XlX,I 7 ilHT, 3 F 13 . 3 ) B 6 ufl 

, 655 F 0 RHAT( 48 X,I 16 , 2 H • , F 14 . 3 , ? 19 , J ) ' - ‘ g^OO 

1170 . 556 F 0 RMAT(/ 14 H PROGRAM PASS II, 19 H LAST PLATE DZO IS E 14 . 7 , 12 H Oltfl BfaOO 

IDEO BY E 14 . 7 ) B 6 C(j 

557 F 0 RMAT( 1 X,I 9 , 2 X,E 14 . 7 , 2 X,£ 14 . 7 , 2 X,E 14 . 7 ) 8600 

558 FORKAT( 4 X,F 6 . 0 , 2 X,F 2 ,B, 2 X,- 2 .G, 61 X,I 1 ) ' 8600 

^ 560 FORMArt///lX,I 7 ,lMT,F 39 , 3 l B 60 o 

1175 561 F 0 RHAT( 1 X,I 9 , 1 X,F 12 . 5 . 1 X,- 12 . 3 , 1 X,F 12 . 3 , 33 X,'I 1 ) ' ‘ A 6 C 0 

.. 562 FORHAMIH //I ggag 

567 F 0 RMATI 723 H POSITION HEIGHT VARIES/) ' ' C 600 

568 F 0 RMAT</ 25 H POSITION HEIGHT CONSTANT/) C 600 

569 FORHAT<lHl, 6 Aia) ’ ■ - - 

1160 1 570 format (/ 40 H OVER 250 HEIGHTEO OBJECTS ENTE RE D— STOP) 3600 

572 FORMATHBIH NUMBER OF PHOTDS IN ORDERED STRI= EXCEEDS 'EO — 57 qP) ' B 6 C 0 
__ .. „ _ 573 F 0 RmAT< 54 H control EGJATION XY Z <Q = N 0 CONSTRAINT)) B 6 GC 

574 FORMAT< 2 X,I 6 ,aX,I 2 , 8 X,I 2 , 5 X,I 2 ) 8630 

- 575 FORHAT< 60 X, 13 H RMS 03 JECT= E 14 . 7 , 8 H MICRONS) 3 {, 0 ,j 

1165 576 F 0 RMAT( 15 H FL CONSTANT = E 14 . 7 , 24 M MAXIMUM ITERATIONS = I 2 , 33 H B 600 

1 RESIDUAL LIMIT (MICRONS) = £ 14 . 7 ) 0600 

577 F 0 RHAT(. 28 H RESOLUTION HEIGHTS HT 7 l 2 '= £147 771 GH"’HT 12 P '= ElVV?/) 8600 

... , 578 F 0 RHAT<iHl, 26 H PLATE 2 HS (MICRONS) /) 8600 

SNO 0600 


1170 

8 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
U82 

1183 

1184 

1185 

1186 

1187 

1188 
1169 

19 

1191 

9 

10 

1194 

1195 

1196 

1197 
1193 

1199 

1200 
1201 
1202 

1203 

1204 



SUBROUTIKE DBUf 


CDC 6630 FTN V3.0-32A 0PT=2 0<./Q2/75 


2O.'.0.S4 


-Jd 


SUBROUTINE DBUF CNORECS) _ „ _ „ 

C ' ~ 

C THIS PROGRAM IS SET UP FOR A DOUBLE BUFFER AREA ALTHOUGH IT IS USING 
C KAIT MODE INSTRUCTIONS. TO 3E CHANGED IN THE FUTURE. 


_S C. 



COMMON /COOUF/ LENGT, HEX! , I FI RS T , I XBUF ♦ 3U- F( 1 0 2 4) , IQ, KS 
DATA (ICOUNT=C) 

DATA II T 3=24012005360 00 DO 0 0 00 00) , ( I £F = 10 00 0 3 0 0 C 00 00 0 6 ) , 

1 tiPR=2ooooacocooooooDB) 

IF (NORECS) 70,10,70 

CALL LTRIO ( I D, 1113 tOUFF ! I F IRST ) , BUFF ( IFIRST ♦ 5 1 1) ,KS I 

IF IKS .LT. 0) GO TO 76 _ 

IF ( (KS.ANO.IPR) .NE. 0) ^RINT 100, IQ 
IF ( (KS .AND.IEF) .N£. Q) 30 TO 60 
IFIRST = HOO( IFIRSTf512,13 24) 

ICOUNT = 0 
RETURN 

ICOUMT = ICOUNT F 1 

IF (ICOUNT .UT. 2) GO TO 10 ” " 


PRINT 64 


CALL LTRIO ( I T3 ,115B, A , 6, J S ) 

_ ___ STOP 

7C IF (M0D(LEN0r,512)) 79,74, 70 

25 74 call LTRIO ( I Q, 11 IB, BUFF ( I F IRST ) ,BUFF ( IFIRST ♦ 511 KS I 

IF (KS .LT. C) GO TO 76 ■ ■ 


IF ( (KS .AND.IPR) .NE. 0) ’RINT 100, IQ 
IFIRSr = MOOaFIRSTF5l2,n24) 

_ 79 LENGT = H0D(LtNGrF8»H0RED3, 1024) 

30 IX9UF = NEXT 

NEXT = LENGT F 1 

■return 

76 PRINT 77, KS 
STOP 

35 100 FORMAT (7» TROUBLE IN INPUT TAPE * 020 /) 

'64 FORMAT(»l JOB TERMINATED— END OF DATA*) 

^77_F0RHAT(»i JOB ABORTED-- STATUS WORD ♦ 020 ) 

ENO 


DQUF 

2 

‘ OBUF" 

' " 3 

DBUF 

4 

OBOF 

5 

DBUF 

6 

DBUF 

■ 7 

DBUF 

a 

DBUF 

9 

OBUF 

10 

OSUF 

11 

OBUF 

12 

OBUF 

13 

OBUF 

14 

DBUF 

" ' 15 

OBUF 

16 

OBUFl 

1 

OBUF 

18 

OBUF 

19 

DdUF 

20 

OBUF 

21 

DBUF 

22 

OBUF 

23 

OBUF 

24 

duUfi 

2 

OBUF 

26 

OBUF 

27 

OBUF 

26 

OBUFl 

3 

oeurt 

4 

OBUF 

31 

DBUF 

32 

■ OBUF 

33 

DBUF 

34 

DBUF 

35 

DBUF 

36 

OBUF 

37 

OBUF 

38 

” OBUF" 

39 



SUBROUTINE OBUF 


COG 663 0 FTN VStO-SZ** OPTs’ Uit/I}2/7S 20.<*a.5<t 


PAGE 





„.S.UOROUTINE„OOUF^(NORECS) OBUF 2 


5 

C 

, c 
c 

. c 
c 

THIS PROGRAM IS S2T UP FOR A 30UBLE BUFFER AREA ALTMOUGH 
WAIT MODE INSTRUCTIONS. TO Bi CHANGED IN THE FUTURE. 

COMMON /009UF /UENGT , NEXT. I FIRS Tt IXBUF .8U" F ( 1 0 2 A) , I P, 

IT IS USING 
ENOFLQ 

OQUF 

OBUF 

OBUF 

OBUF 

OBUF 

OBUF 

3 

4 

5 

6 

7 

8 




IF (NORECS) 10,40,1.0 


OBUF 

■ "■ 9 



IQ 

UENGT = LENGT ♦ NORECS 


OBUF 

10 

10 



IF «LENGT .GT.faA) GO TO 20 


oauF 

11 




IX3UF 5 NEXT 


OBUF 

12 




GO TO 30 


OBUF 

13 



?n 

CALL LTRIO (IP,U2Q,BUFF(If IRSn,8UFF(IFIRSr*511) ,KS) 


OBUF 





IF IKS .Lr. 0) GO TO 80 


OBUF 

IS 

15 



IX9UF = IFIRSr = MOD(IFIR3rt312,i02Al 


08UF1 

1 




LENGT = NORECS 


OBUF 

17 



30 

NEXT = IXBUF f NORECS*8 


OBUF 

18 




RETURN 


OBUF 

19 



An 

ItlOEX = IFIRST + (LENGT-NORECSI »8 - i 


OBUF 

20 

20 



call LTRIQ aP,U2B,aUFF(IFIRSn ,BUFF(IN0EXI ,KS1 

“ ‘ ^ ' 

OBUF 

21 




IF (<S .LT. 01 GO TO 80 


OBUF 

22 




ENDFLQ = 0- 


OBUF 

23 




RETURN 


OBUF 

2A 



so 

PRINT tCC,KS 


OBUF 

25 

25 



STOP 


OBUF 

26 



IOC 

FORMAT <♦! JOS ABORTED-- STATUS WORD ‘.DEn 


OBUF 

27 




ENO 


OBUF 

26 







